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ABSTRACT 


The  ’modified  dipole'  has  its  origin  in  the  consideration  of  the  ge¬ 
neral  properties  of  a  satellite  antenna  which  bears  great  resemblance  to 
a  dipole  modified  to  incorporate  at  the  center  a  conductir,  ;  volume  which 
is  used  to  radiate  electromagnetic  waves  and  to  house  a  power  supply  and 
radio  frequency  generators,  etc.  The  object  of  this  research  is  to  pursue 
a  theoretical  and  experimental  exploration  of  the  effects  induced  by  the 
presence  of  the  conducting  volume  on  the  antenna  performance,  i.  e.  ,  in¬ 
put  characteristics,  current  distribution  along  the  surfaces  of  the  entire 
radiating  structure  and  radiation  properties. 

In  Volume  I  a  mathematical  model  consisting  of  a  perfectly  con¬ 
ducting  sphere  from  which  project  the  ends  of  a  thin  biconical  antenna  is 
chosen  to  simulate  the  actual  sphere- cente red  thin  dipole.  The  conical 
antenna  is  driven  at  its  junction  with  the  sphere  by  a  rotationally  symmetric 
electric  field  maintained  across  the  gap  by  a  biconical  transmission  line 
excited  by  the  TEM  mode.  The  attractive  features  of  this  model  include 
the  fact  that  it  has  surfaces  that  permit  a  simple  specification  of  boundary 
conditions  and,  hence,  a  rigorous  formulation  for  the  electromagnetic 
fields  and  a  shape  such  that  its  properties  should  come  reasonably  close 
to  those  of  a  modified  cylindrical  antenna  as  the  cone  angle  becomes  quite 
small. 

The  measurements  of  both  input  admittances  and  current 
distributions  on  modified  dipoles  (with  either  conical  or  cylindrical 
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antenna  projecting  from  the  sphere)  are  also  presented  in  Volume  1. 
Comparisons  were  also  made  between  modified  conical  and  cylindrical 
antennas  with  the  same  sphere  radii  and  antenna  heights.  The  radius  of 
the  cylindrical  antenna  is  the  same  as  the  smaller  end  of  the  cone.  The 
fact  that  the  admittance  curves  for  modified  cylindrical  and  conical 
antennas  involve  only  slight  shifts  suggests  that  by  introducing  an 
equivalent  antenna  length  that  is  a  little  longer  than  the  actual  physical 
length  of  the  conical  antenna  a  good  approximation  is  obtained  for  the' 
cylindrical  antenna. 

An  infinite  set  of  algebraic  equations  was  solved  numerically  in 
Volume  II  for  small  cone  angles.  Comparisons  were'  made  between  the 
modified  conical  antenna  and  its  limiting  biconical  antenna  which  provides 
both  an  extrapolatory  numerical  check  for  the  modified  conical  antenna 
with  shrinking  central  sphere  and  an  understanding  of  the  underlying 
physical  phenomena.  Theoretical  and  experimental  results  are  in  very 
good  agreement. 


Accepted  for  the  Air  Force 

Joseph  R.  Waterman,  Lt.  Col.  USAF 

Chief,  Lincoln  !  .ahomtorv  Project  Office 


1.  INTRODUCTION 


The  formal  solution  for  the  b,  1  s  and  with  them  the  determination  of 

k 

the  driving  point  admittance,  the  current  distribution  and  the  far-field 
pattern  of  the  modified  conical  antenna  with  half  cone  angle  6  ^  and 
central  sphere  of  radius  b  is  contained  in  the  infinite  set  of  linear 
equations,  (1-57)  and  (1-64)  which  appear  in  Volume  I  of  this  report.  No 
general,  exact  solution  of  this  set  of  equations  is  available  for  cones  with 
arbitrary  angles  0^.  Hence,  a  general  formula  for  the  input  admittance 
of  the  modified  dipole  cannot  be  provided.  Thus,  although  the  analysis  of 
the  problem  is  not  restricted  in  its  formulation,  mathematical  limitations 
make  an  exact  solution  unavailable.  However,  if  the  half  cone  angle  0^ 
is  sufficiently  small  (i.  e.  ,  0  ^  <  4°),  which  is  of  special  interest  in  the 
present  problem,  a  numerical  solution  of  (1-57)  and  (1-64)  can  be 
carried  out.  In  this  problem  an  exact  result  can  be  obtained  only  by  an 
infinite  sequence  of  steps.  Since  the  process  necessarily  has  to  be 
truncated  after  a  certain  finite  number  of  steps,  truncation  errors  are 
unavoidable.  There  are,  however,  three  sources  of  truncation  errors  in 
the  approach  taken  to  this  problem.  First, there  is  the  error  introduced 
when  the  infinite  matrix  is  truncated  to  one  of  finite  dimensions.  Secondly, 
error  is  introduced  by  truncating  the  infinite  sum  (over  v)  in  the 
calculation  of  each  matrix  element.  Thirdly,  the  double  series  re¬ 
presenting  the  input  admittance,  current  distribution  and  near  fields  is 
twice  truncated;  once  over  v  and  once  over  k.  In  addition,  because  of 
the  formidable  amount  of  computation  required  to  solve  the  truncated 
infinite  matrix,  small  round-off  errors  may  cause  errors  in  the  computed 
solution  out  of  all  proportion  to  their  size.  Above  all,  the  truncated 
system  of  linear  equations  has  to  be  solvable  (i.  e.  ,  the  coefficient 
matrix  must  be  non-singular)  before  problems  of  accuracy  and  con¬ 
vergence  of  the  computed  solution  can  be  considered. 

An  exact  error  analysis  relating  different  sources  of  errors  to  the 
final  results  is  very  difficult,  if  not  impossible,  when  many  equations  are 
involved.  It  cannot  be  made  here.  Hence,  after  all  possible  precautions 
have  been  taken  to  prevent  error  from  growing  out  of  proportion  and  to 
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preserve  the  accuracy  of  the  final  solution,  the  computed  results  are  to 
be  justified  by  comparison  with  expe ri mental  measurements. 

A  brief  review  of  the  major  steps  taken  in  the  course  of  the 
numerical  calculation  and  their  associated  theories  follows: 


(A)  Existence  of  the  Solution 

Once  the  infinite  set  of  linear  equations  is  truncated,  a  set  of  n 
simultaneous  linear  equations  in  n  unknown  is  at  hand.  It  is 

n 

X  a.  .  x  .  b.  (2-  la) 

j=i  1J  J  1 

Equation  (2- la)  is  conveniently  written  in  matrix  form  as  follows: 

AX  b  (2-lb) 

T 

where  A  la..l  is  the  nxn  matrix  of  coefficients,  X  (x .  ...  x  )  and 
T  L  ij  1  ’In 

b  (b^  .  .  .  b  )  with  T  denoting  the  transpose.  Let  A^  denote  the 
n  x  ( n-b  1 )  matrix  which  has  the  column  vector  b  appended  as  an  (n-tl)st 
column  to  A.  The  rank  of  any  matrix  A  is  r(A).  Then  the  basic  theory 
on  the  existence  of  a  solution  of  (2-1)  states  that  the  system  of  linear 
equations  has  a  solution  if  and  only  if 

r(  A)  -  r(Ab)  (2-2) 

This  familiar  theorem  may  be  proved  by  using  Gaussian 
elimination.  This  will  not  be  pursued  here  [  1  J  —  [  2  ]. 

The  theoretical  proof  of  the  existence  of  a  unique  solution  for  the 
simultaneous  linear  equations  does  not  in  any  way  guarantee  the  correct¬ 
ness  of  a  computed  solution  of  (2-la,  b).  There  is  no  hard-and-fast  rule 
in  carrying  out  a  numerical  solution  for  simultaneous  linear  equations. 
Different  approaches  need  be  employed  for  different  matrices  with 
different  characteristics;  e.  g.  ,  direct  methods  are  usually  used  for 
matrices  which  are  filled  but  not  large,  while  iteration  methods  are  use¬ 
ful  for  matrices  that  are  sparse  and  filled  etc. 
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(B)  111-  Condition 


In  numerical  calculations,  a  class  of  matrices  of  coefficients  is 
encountered  which  is  called  ill-conditioned  such  that  no  matter  how 
accurately  it  is  calculated,  the  solution  may  still  be  grossly  in  error.  If 
the  matrix  A  is  normalized  so  that  the  largest  term  in  magnitude  has  the 
order  of  magnitude  of  unity  and  is  such  that  A  ^  contains  some  very 
large  elements,  the  matrix  and  therefore,  the  system  of  equations  is 
ill-conditioned.  Conversely,  if  the  largest  element  in  magnitude  of  A 
has  the  order  of  magnitude  of  unity,  the  matrix  may  be  said  to  be  well- 
conditioned.  This  can  easily  be  seen  by  writing  (2- lb)  in  the  following 
fo  rm: 


X  A  lb 


(2-3) 


Suppose  that  A  ^  has  very  large  elements,  one  of  which  is  a^  A_/|a|. 
I  A  I  denotes  the  determinant  of  A  and  A.  .  is  the  cofactor  of  a.  ..  The 
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ij 


1J 


assumption  that  a^  is  large  means  that  A_  must  be  large  relative  to  |A|. 

Since  one  of  the  terms  in  the  expansion  of  |  A  |  about  the  ith  row  and  jth 

column  is  a.  .A..,  a  small  error  in  a.,  may  cause  a  large  relative  error 
iJ  iJ  iJ  .  -1 

in  |  A  I  and  therefore  a  large  relative  error  in  a^..  .  This  in  turn  can 

cause  a  large  relative  error  in  X.  Similarly,  a  small  change  in  an 

element  of  b  could  cause  a  large  change  in  X. 


(C)  Truncation  Error 

The  errors  caused  by  truncating  the  infinite  sum  by  finding  an 
approximation  (or  an  upper  bound)  to  the  remainder  in  each  series  after 
truncation  will  be  examined.  Having  found  these,  an  estimate  could  be 
found  of  the  relative  error  in  each  matrix  element  by  comparing  the 
remainder  with  the  corresponding  truncated  sum.  Finally,  the  remainder 
estimate  is  added  to  the  corresponding  matrix  element,  the  system  is 
solved  and  its  solution  compared  with  that  of  the  unaltered  matrix.  This 
comparison  yields  an  estimate  of  the  error  one  might  expect  in  the  field 
pattern,  current  distribution  and  input  admittance  etc. 
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(D)  Intermediate  Round-offs 

One  way  of  checking  the  round-off  error  in  the  course  of  com¬ 
putation  is  by  substituting  the  calculated  values  into  the  left-hand  members 
of  the  original  equations.  The  presence  of  the  deviations  between  the 
resultant  members  and  original  right-hand  members  serves  to  indicate 
the  presence  of  errors  due  to  intermediate  round-offs.  Again,  the 
relationship  between  the  magnitudes  of  these  deviations  and  the  magnitude 
of  the  errors  in  the  solution  column  is  not  a  simple  one.  Use  will  be 
made  of  this  simply  as  a  check  if  the  round-off  errors  are  too  far  out  of 
proportion  to  afford  a  reasonable  solution. 

(E)  Convergence  of  the  Solution 

Before  considering  the  convergence  of  the  solution  of  the  system  of 
linear  equations,  error  checks  were  made  for  several  cases  with  various 
combinations  of  the  extreme  sizes  of  antennas  and  inner  spheres.  The 
maximum  size  of  matrices  has  been  studied  so  that  the  different  kinds  of 
errors  mentioned  above  can  be  kept  within  one  per  cent.  It  has  been 
found  that  for  most  cases  a  40  x  40  matrix  is  still  well  within  the 
acceptable  range  and  the  solutions  show  rapid  convergence  when  the 
matrix  size  exceeds  15  x  15. 

The  convergence  of  the  solution  of  the  system  of  linear  equations 
will  be  dealt  with  empirically  by  solving  the  same  problem  with  matrices 
of  increasing  size  and  examining  the  successive  solutions  for  conver¬ 
gence. 


Z.  ROOTS  OF  THE  CHARACTERISTIC  EQUATION  L  (0n)  0 

v  0 

The  application  of  the  boundary  condition  given  by  equation  (1-35) 
in  Volume  I  together  with  the  property  of  anti-symmetry  of  E^  with 
respect  to  the  equatorial  plane  (i.  e.  ,  0  t\  f  Z  -  plane)  leads  to  the 

following  characteristic  equation  [see  equation  (1-30)  in  Volume  1]: 
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W  Lv(»-fl0>  0 

which  is  independent  of  the  presence  of  the  conducting  sphere  at  the 
center.  Once  the  half  cone  angle  0 ^  is  fixed,  an  infinite  sequence  of 
discrete  values  of  v  is  obtained  from  the  above  equation.  The  present 
analysis  will  be  restricted  to  small  half  cone  angles  (i.  e.  ,  0^  <  4). 

In  the  neighborhood  of  the  point  9  7 r,  Legendre’s  functions  can  be 

expressed  in  the  following  way  [3]  : 

sin  vjt  ?  9 

Pv(-cos0o)  =  — - -  F(-v,  v+  1;  1;  sin  — )  •  H(v,  0Q) 


00  (-v)*  •  (- v4a- I)(v4l)#  •  (v4a) 


9 


4  2 
a  1 


<j)  (v,  a)(sin -^)2a  (2-4) 


( CL !  >  ‘ 


where  9  ij  -  9  is  in  the  neighborhood  of  zero 


H(v,  9)  2  sin  —  +  l//(v4l)  4  4  2\ 


(2- 5a) 


and 


T(c)  oo  r(a4n)T(b4n)  n 


F(a.,  h;  c;  z)  r(a)r  (b)  t(c+n) 

n-0 


n! 


(2-5b) 


which  is  the  Gauss  hypergeometric  series  whose  circle  of  convergence  is 
the  unit  circle  |z|  1. 


v+1 


+  4-  +  - 

v4r  -  v 


1 


- v4  r - 1 


-  2  (y  +  7  +  *  •  •  +  7) 


( 2- 5c) 


and  l//(z)  is  the  logarithmic  derivative  of  the  Gamma  function;  i.  e.  , 


l//(z) 


I  '(z) 

r  (z) 


( 2- 5d) 


or 


-  V  +  (Z-  1 ) 


Cl) 

v 


1 

(nl  1  )(n+  z) 


(2-5e) 


n  0 


where  y  denotes  Euler’s  or  Mascheroni’s  constant  equal  to  0.  57722. 

The  expression  (2-4)  will  be  single -valued,  on  account  of  the 
cross-cut  of  the  plane  of  \i  cos  6  along  the  real  axis  from  -1  to  -oo. 

Using  the  reflection  formula  for  the  ^//-function,  namely, 

\J;(z+\)  t//(-z)  -  77  cotTrz  (2-6) 

(2-4)  is  rearranged  in  the  following  form; 

E^(cos0q)  =  sin  V77  G(v,  0  (2-7a) 


sin  V77  [G  ^  (v,  0  q)  +  G^(v,0q)] 

( 2 - 7b) 

where 

Gi(v’V 

2  % 

F(-v,  v+1;  1;  sin^  ~)  •  F(v,  6 Q) 

(2-8a) 

F(v,0q) 

C(0n)  I  -  if/  (v+l)  -  tan 

(2-8b) 

In  ( 2-8b) 

C(B0) 

-  “  [ln(sin  —)  +  \] 

( 2- 8c ) 

G2(V,  0q) 

oo  (v-a+1  )•  •  v(v+l)*  •  (v+a)  <t)(v,^ri) 

2  (-dq - - — j -  - - 

,  a!  -  a!  77 

a  1 

(2-8d) 

Thus,  finding  the  roots  of  the  characteristic  equation  is  now 
equivalent  to  setting  G(v?0q)  0.  Before  actually  carrying  out  the 

numerical  calculation  of  the  roots  of  (2-7a),  it  is  instructive  first  to 
survey  the  equation  qualitatively.  Some  general  properties  of  these 
roots  regardless  of  the  cone  angle  6 ^  (as  long  as  it  is  small)  can  be 
derived  by  looking  at  (2-7a,b)  more  carefully.  The  vanishing  of  (2-7a) 
can  also  be  accomplished  by  finding  the  intersections  of  two  functions, 
i.  e.  ,  -G i  (v,  0q)  and  G^(v,  6^)  with  v  as  a  variable. 
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It  is  convenient  to  look  into  the  restrictions  upon  v  and  0^  such  that 
the  numerical  solution  of  G(v,  0^)  0  is  feasible  and  a  sufficient  number  of 

roots  is  obtained  to  provide  accuracy  in  the  calculation  of  the  matrix 
elements,  the  input  admittance  and  the  current  distribution.  It  has  been 
found  that  both  the  functions  Gj  and  G  y  are  well  behaved  (i.  e.  ,  the  series 
converges  quite  rapidly)  and  the  above-mentioned  requirements  were  met 
over  the  range  of  interest  so  long  as  v  is  not  too  large  in  magnitude  and 
the  cone  angle  0^  is  sufficiently  small  so  that  the  product  of  these  two 
quantities  does  not  exceed  the  order  of  magnitude  of  unity  [i.  e.  , 

0o 

v  sin  ^  0(1)]. 


The  function  -G^(v,0^)  behaves  essentially  like  tan  with  its 

values  continuously  adjusted  by  first  adding  to  it  a  slowly-decreasing  quantity 
2 

[i.  e.  ,  -C(0  )  -  —  \p  (v+ 1 )  ]  and  then  multiplying  it  by  a  slowly- va  rying 


0 
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2  °0  2  % 

function  [i.  e.  ,  F( - v,  v+  1 ;  1 ;  sin  ***  ( 1  -v(v+l )  sin  +  •••)].  Therefore, 

-Gj(v,0q)  is  a  continuous,  monotonically  increasing  function  within  each 

equally- spaced  segment  which  has  a  total  length  of  two  and  the  odd  integers 

as  its  end  points.  G^(v,  0^)  becomes  singular  in  the  neighborhood  of  the 

end  points  and  changes  sign  whenever  it  crosses  the  end  points.  The  ^ 

2  2  U  0 

derivative  of  G9(v,0  )  with  respect  to  v  is  proportional  to  (4-3v  sin  -r- )• 

0  L  U  0Q  C 

(v  sin  which  is  positive  as  long  as  v  sin  — y-  ^  1.  So  G^(v,  0^)  is  also  a 

monotonically  increasing  continuous  function  over  the  range  of  interest. 

Accordingly,  the  functions  -Gj(v,  0  )  andG^(v,0^)  intersect  only  once  in 

each  segment.  Therefore,  the  form  of  the  roots  is  an  odd  integer  plus  a 

small  number  which  increases  with  v  and  0 q. 

The  simplification  usually  employed  [4]  to  obtain  an  approximation 
of  the  roots  of  (2-7)  is  actually  to  take  the  leading  term  of  G^(v,  0^) 

[i.  e.  ,  F(v,  0  )  ]  and  retain  only  the  two  terms  in  F[i.  e.  ,  tan  (v7t/2)  and 

d  0  Q 

ln(0g/2)]  under  the  as  sumption  that  ln(sin  — y-)  approaches  infinity  for 
sufficiently  small  angles  0^.  Therefore,  the  approximate  roots  are 


v 

app. 


k  + 


ln(  sin  — — ) 


(2-9) 
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This  approxi mation  is  in  error  due  to  the  fact  that  even  for  half  cone 

°0 

angle  0^  as  small  as  0.  1  ,  |  ln(sin  |  (  7.  044  in  this  case)  is  still  of  the 
same  order  of  magnitude  as  the  two  neglected  terms,  i.  e.  ,  (v-fl)  and  y. 

Furthermore,  ln(sin— ^-)  becomes  even  smaller  than  i//(vfl)  as  one  proceeds 
to  find  larger  v's.  In  this  range,  the  higher-order  terms  of  G  become 
comparable  in  magnitude  and  can  no  longer  be  neglected. 

In  Fig.  2-1,  the  results  from  the  approximate  formula  (2-9)  are 
compared  with  those  calculated  from  the  exact  formula  (2- 7a)  for  various 
cone  angles.  Generally  speaking,  (2-9)  provides  quite  a  poor 
approximation  of  the  accurate  formula  except  for  very  thin  cones  for 
which  the  approximate  formula  gives  the  correct  trend  but  the  results  are 
shifted  from  the  accurate  ones  by  a  constant. 

The  numerically  calculated  v’s  were  checked  by  substituting  them 
in  a  subroutine  which  generated  Legendre  functions  of  fractional  orders. 

At  0  the  generated  values  of  L v($q)  (which  are  equal  to  zero 

\  —  5 

theoretically)  are  of  the  order  of  10  which  indicates  the  accuracy  of  the 
solution. 

3.  COMPUTATION  OF  FRACTIONAL  ORDER  LEGENDRE  FUNCTION 

The  computational  difficulties  with  the  fractional  order  Legendre 
functions  lie  in  the  following  facts: 

(a)  The  Legendre  function,  defined  in  the  interval  (-1,  1),  can  be 
expressed  by  the  hypergeometric  series  with  various  elements,  i.  e.  , 

2  ^  0 

(cos  0 )  F(-v,  v+1;  1;  sin  (2-10) 

The  above  expansion  can  be  employed  for  the  calculation  of  values  of 

the  function  P  (cos  0)  only  in  a  limited  region  of  variation  of  the  argument 

0  v 

sin—  and  order  v.  With  an  increase  in  v  and  0,  the  convergence  of  the 
series  grows  worse,  which  makes  it  of  little  use  for  calculation  with 
large  values  of  v  and  0. 


FIG.  2-1  ROOTS  OF  CHARACTERISTIC  EQUATION  L„(0o)=O 
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(b)  The  process  in  generating  higher-order  functions  by  a  re¬ 
currence  relation  is  unstable,  since  the  generation  is  carried  out  with 
rounded  values  which  propagate  in  the  recurrence  process  of  increasing 
order  and  hence  cause  loss  of  significant  figures. 


The  Legendre  function  P^Jcos  0)  assumes  various  representations 
suitable  for  numerical  calculation  in  different  ranges  of  v  and  0  with  the 
aid  of  definite  and  contour  integrals.  In  this  work  use  is  made  of  re¬ 
presentation  of  P  (0)  in  the  form  of  Mehler-Fock,  namely, 
v 


P  (cos  0) 

V 


_2  cos(  Vi  1  /Z)t  _ 

71  q  Z(cost  -  cost)) 


(2-11) 


which  provides  (under  different  conditions  for  v  and  d)  rapidly  converging 
series  for  numerical  calculation. 

Equation  (2-11)  can  be  expanded  in  positive  powers  of  (v+l/2)  by 
substituting  in  it  in  place  of  cos(v+l/2)t  its  series  expansion.  After 
changes  of  variables  and  rearrangement  of  terms,  one  obtains 


P  (cos  0 )  E  B  (v+  1/2) 
v  n  m 

m  0 


2m 


(2-12) 


where 


Bni  t-1* 


m  2 


(2m+l)  1  / Z77  [sin  ^(psinx)]^m 


71  (2m)!  0  ( 1  -p^  sin^  x) 


dx  (2-13) 


and 


p  sin 


(0/2) 


(2-14) 


In  the  last  integral,  the  function  under  the  integral  sign  becomes 
infinite  only  when  p  1,  that  is  for  0  77.  Consequently,  outside  the 

neighborhood  of  the  point  0  77,  it  is  possible  to  calculate  this  integral  by 

the  usual  numerical  methods. 

One  limitation  in  using  (2-12)  for  computing  Legendre  functions  is 
its  slowness  in  convergence  as  v  becomes  larger. 
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An  expansion  of  Legendre’s  function  of  (2-11)  can  be  expressed  in  a 
series  of  Bessel’s  function  [5],  namely, 

oo 

P  (cos  0)  r  2  a  (0)J  .  [(v+l/2)0]  (v+l/2)"m  (2-15) 

m  0 


where  the  a  's  are  elementary  functions,  regular  in  0  ^9  ^  7T,  and 


a  (0)  -  — : — zr 

0  sin  0 


(2- 1 6a) 


a  (9)  (2m- 1)!!  an(0)  0nty  (9) 

m  0  x  m 


(2—1 6b) 


\ 1/  (8)  is  a  combination  of  functions  of  <p  which  is  obtained  from 

m  2  2m 
2(cos  t  -  cos  8)  in  terms  of  a  series  of  (0  -  t  ),  i.  e.  , 

2(cost-cos  9)  -  (9Z  -  tZ)  ~JL  (i  +  2  <02-t2)v-1<|>  (0)]  (2-17) 

m  2 


where 


4>  (0)  =  —  1 

rm 


m!  2 


J  1  (9) 

m-  2 


,  Jm-1)  fl(m-l) 


J,  (0) 


(2-18) 


Taking  the  inverse  of  the  square  root  of  (2—17),  one  obtains 

oo 


1 


1 


9 


7  7  S  (92  -  t2)n>  (0)  (2-19) 

/2(cost-cos0  \  8  -t  sin  0  m-  0 


Therefore,  the  combination  of  (2-17)  and  (2-19)  gives 


^0 

*1 

^2 

^3 


1 


cj)  3<j> 


2  f  8 


§>A  3 

-  X  +  4  ‘h  h 


2  r2 
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^4  2  ^5  +  8  (<b  +  Z<f)2<f)3)  "  16  ^3^2  f  128  ^2 

^5  “  {  *6  f  !  (<M5  +  *3+4>  -  T6  <3^a  +  ^\) 

.35,3.  ,  63,5 

+  32  ^2^3  +  250  ^2 


By  using  the  recurrence  relation  for  Bessel’s  functions  of  the  first 
kind,  namely. 


Jn-l(z,  +  Jn+l(z) 


2n  t  /  . 

—  J  (z) 
z  n 


(2-21) 


and  substituting  (2-21)  into  (2-15),  we  finally  obtain  an  expression  for 

P  (cos  0); 


P  (cos  0) 

V 


J  [(v+l/2)fl][a  -  - - - 2 

(v+l/2) 


4a 


(v+l/2 )  9 


+  •*] 


+  J1[(vtl/2)6]  - 

(v+l/2) 


(a  0-  2a?)  8a 

— ? - zL.  +  - __ — _  -|-  .  . 

(v+1/2)30  (v+1/2)502 

(2-22) 


At  first  glance,  (2-22)  does  not  assume  any  advantage  over  either 
of  the  more  general  expressions  (2-10)  and  (2-11).  But  the  asymptotic 
expansion  for  large  values  of  v  is  a  convenient  representation  of  the 
function  P  (cos  0)  from  the  computational  point  of  view.  It  gives  a  good 
approximation  to  the  functions  themselves  and  can  also  be  used  for  the 
computation  of  their  derivatives  and  integrals  with  respect  to  the 
variable  v.  Furthermore,  the  coefficients  of  the  representation  are 
associated  with  functions  of  spherical  Bessel's  functions  and  integral 
order  Bessel's  functions  of  the  first  kind  and  J^,  which  are  axailable 
as  built-in  functions  on  the  IBM  360/65  machine  in  use  at  the  Harvard 
Computing  Center.  Hence,  (2-22)  is  both  a  convenient  and  accurate 
formula  for  computation. 
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In  summary,  it  is  necessary  to  use  a  suitable  representation  for 
Legendre’s  functions  of  fractional  order  in  different  regions  in  the  v-0 
plane  for  purposes  of  computation  with  desired  accuracy.  Therefore,  as 
shown  in  Fig.  2-2,  the  v-0  plane  is  divided  into  four  regions.  In  region  I, 
the  hypergeometric  expansion  (2-10)  is  generally  suitable  as  long  as  the 
angle  is  small  and  the  order  is  not  large;  also  the  product  of  v  and 
sin  (0/Z)  must  be  less  than  one.  In  region  II,  (2-4)  is  a  suitable  expansion 
in  generating  Legendre's  function  with  the  same  restrictions  upon  v  and 
0  as  in  region  I.  In  region  III,  use  is  made  of  the  expansion  (2-12)  as 
long  as  the  order  is  not  large  (v  <  10).  Finally,  in  region  IV,  (2-22) 
provides  a  good  tool  which  is  also  valid  in  the  neighborhood  of  9  -  0  and  7 j 
for  large  values  of  v.  Therefore,  in  computing  L^(0)  P^(0)  -  PV(7T  -0), 

formulae  for  different  regions  must  be  used  simultaneously,  e.  g.  ,  for 
small  0,  formulae  for  region  I  and  II  must  be  used. 


The  quantity  Q  .in  the  second  set  of  linear  equations  is  generally 
r  ic 

a  small  quantity  for  small  half  cone  angle  0q.  For  k/r,  can  easily 

be  expressed  in  a  close  form,  i.  e.  , 

2,  [k(k+l)P;(^)Pk(^)-r(r+l)Pr(lx0)Pk(lx0)] 
rk  "  "  [k(k+l)  -  r(r  +  l)  ] 


(2-23) 


and  for  k  r 


Q 


rk 


where  H(k,  0 


2Hl-^)P^0)Pk(^0)  +  2k(k+l)(l-^)H(k,0o)] 

2% 

)  is  a  polynomial  of  sin  -y  of  order  k,  namely, 


(2-24) 


HM0) 


OO 

2 

m  0 


C  ,  9 

m  ,  .  2  0,m 

7m  (sin  "2” 


where 


m 


m 

2 

n  0 


a  k  A  k 

A  A 
n  m-n 


(2-25) 


(2- 26a) 


FIG. 


it/2. 


2-2  COMPUTATION  OF  FRACTIONAL  ORDER 
LEGENDRE  FUNCTION 
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and 


m 


.  m  T  (kbmi  1 ) _ 

{~i}  T  (mb  1 )  T  (k-mbl) 


(2-26b) 


Hence,  A  "  0  for  m  >  k. 

’  m 

The  generation  of  Bessel’s  functions  of  fractional  order  is  not 
significantly  different  from  that  for  functions  of  integral  order.  It  will 
not  be  discussed  here. 


4.  THE  BICONICAL  ANTENNA 

As  a  preliminary  to  the  later  discussions  of  the  modified  dipole, 
the  problem  of  the  thin,  center-fed  biconical  dipole  will  be  treated  first. 
Based  upon  equations  (1-90),  (1-97)  and  (1-99)  given  in  Volume  I,  a 
numerical  solution  for  the  input  impedance  and  current  distribution  with 
various  cone  angles  and  antenna  heights  can  be  obtained. 

Most  of  the  considerations  here  will  provide  an  opportunity  to 
develop  the  technique  which  will  be  employed  later  for  the  more  com¬ 
plicated  case  of  the  modified  dipole.  Furthermore,  they  will  also  serve 
as  a  numerical  extr apolatory  check  for  the  modified  dipole  with  a 
shrinking  central  sphere. 

Three  sets  of  linear  algebraic  equations  describing  the  same 
biconical  structure  will  be  solved  numerically  for  various  antenna  lengths 
and  cone  angles  with  proper  error  estimate  for  the  elements  of  the 
coefficient  matrices  and  the  terms  with  infinite  summations. 

The  general  terms  of  the  coefficient  matrix  in  (1-90)  include  the  sum 

(Wq)  9v 

[k(k+  1 )  -  v(v+l)  ][r(r+l )  -  v(v+l)]  9(iQ 


g(r,  k,  v) 


( 2v+ 1 ) 
v(v+l) 


(2-27) 
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To  discuss  this  case  for  large  v,  some  approximate  formulae  are 
introduced  for  the  Legendre's  functions  which  have  been  given  by 
Schelkunoff  [6]  provided  0  ^  is  not  near  0  or  77- , 

P  (cos  0}  as  cos(q  -  Ttf  4)  (2-  28a) 

v  (|psinfl)I/2 

and 

cos(q  77  -  q^n  -7t/4) 

P  ( -  cos  0 )  ^  ■  (2-28b) 

V  /  7T  •  q  \  -A  /  ^ 

(“P  sin  a) 

where 

p  [v(v+l)]^^  and  q  -  [(v+l/Z)^  +  l/4]^^ 


Thus? 


s  i  n  ( q  77  /  2  -  jr/4)  sin(7r/2  -  0) 

P  sin())1//Z 


(2-29) 


Using  the  fact  that  L^(d^)  0,  it  follows  that  for  v  very  large 

^  (2v+l)  3v  4  sin  d 

(1“^i0)  v(v+l)  ^  (77  -2^  (2-30) 

The  factor  1  /[k(k*t- 1 )  -  v(v+l )  ]  [r  (r+ 1 )  -  v(v+  1 )  ]  tends  to  be  quite 

large  and  sensitive  to  the  value  of  v  when  v  is  near  k  or  r.  Let  N  be  the 

matrix  size,  then  the  largest  value  that  k  or  r  may  have  is  (2N+  1).  Let 

v  be  the  value  of  v  at  which  the  series  (2-27)  is  terminated.  An 
m 

examination  of  (2-27)  shows  that  the  most  significant  values  occur  when 

v-k  and/or  v- r  are  small  as  v  passes  through  the  entire  range  of  values 

obtained  from  L  ((L)  0.  Therefore,  it  is  necessary  to  require  that  v 

v  0  ’  j  m 

be  greater  than  (2k+l)  and  (2r  +  1)  in  all  computations. 

Let  R(r,k,  v)  denote  the  remainder  such  that 
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oo 

R(r,  k,  v^)  E  g(r,  k,  v) 

V  V 

m 


F rom  (2-27)  and  (2-30),  one  obtains  for 
(i)  off-diagonal  elements,  i.  e.  ,  k  /  r 

A1  (~n - 7— r)  +  A9  (~-  -  — t__) 

1  v+k  v-k-1  2  v+r  v-r-1 


1  (2k+l)(k+r+l)(r-k) 

and 

A  _ l _ 

2  (2r+ l)(k  +  r+l)(k-r) 

From  (2-30)?  one  finds  that  the  separation  of  zeros  of  L^(0 
approximately 


4  sin^ 

s(r-k'v)  ~  <„-2  «0) 


wh  ere 


A  v 


2tt 

(77  ~  2  fig) 


then 


R(r,k,vm)  ^ 


77-2  0 


0 


77 


oo 

g(r,  k,  v)  dv 

V 

m 


2  sin^  0 , 


v  -k-1  v-r-1 

Ai ln  Hrnr*  +  Az ln 


m 


m 


(ii)  For  diagonal  terms,  i.  e.  ,  k  r 

4  sin^  0  1 

g(r,  k,  v)  ^  - -  - 

-20  Q  (2k+ 1 ) 


1  1 


v-k-1  v+k 


2 


(2-31) 

(2-32) 

( 2- 33a) 

(2-33b) 
D>  is 
(2-34) 

(2-35) 

(2-36) 
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thus  for  large  v. 


R(k.  k,  vm) 


Z  s  i  n  2  0 , 


~ 

Z 

k+T 


(Zkfl)‘ 


1  1 

v  +  k  v  - k -  1 
m  m 


v  +1 

i  ,  m  v 

lnlv  -k-l> 
m 


(2-37) 


For  those  terms  with  infinite  summation  over  v  in  (1-97)  and  (1-99), 
the  approximate  formulae  are  introduced  for  large  v  and  v  >>  x,  i.  e.  , 
Meissel's  formula  [7] 


Jv(x) 


(v+l/Z) 

/  X  exp 

_ 


(v+l/Z)2  -  x2  -  (v+l/Z) 


1  -  (  -  )Z 


1/4 


1  + ( 1-x  /v  ) 


2/  2.1/2 


V 


(v-x  »  V 


-1/3, 


(2-38) 


Actually,  since  v  is  much  larger  than  x,  a  further  approximation 
can  be  made,  namely, 


[(v+l/Z)2  -  x2]1  2  - 


X 

Z^ 


1  - 


X 


4  v 


[!+  f  1  -  (f)2]  -  2Vexp(-^) 


thus . 


Jv(x) 


(vf  1  )  ,  X 

l/2  x  exPl- 


4(v+l  Tl\ 


z(w  1  ]  r  (v+3/z) 


(2-39) 
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Similarly, 

J’  (x)  —  J  (x) 

V  X  V 


(2-40) 


Again,  let  gj(r,k,v)  and  g^(r,k,v)  denote  the  remainders  for  the 
first  and  second  set  of  linear  equations,  respectively;  one  finds  for: 

(i)  the  off-diagonal  elements,  i.e.,  k/r 


gj(r,  k,  v) 


and 


R2(r,  k,  v) 


9v 

2v+l  J'.(|30h)  Bu0 

v(v+l)  Jv(Pf)h)  [k(k+ 1)  -  v(v+l )  ][r(r+ 1 )  -  v(v+ 1 )  ] 

(2-41) 


gj(r,  k,  v) 


Jv(f30h)  2 

jMP0h) 


(2-42) 


For  v  large  and  v  >"  P(^h,  one  gets 


•  2  A 

sm  u 


RT<r>k>  vm} 


0 


77  (3Qh  (k2-rZ) 


ln( 


2  2 
v  -  r 
m 

2  v2 
v  -k 

m 


•  v  +k 

)  +  £ln(-2L_) 

V  -k 

m 


.  v  +  r 

-  ln( — — - ) 

r  v  -  r 

m 


(2-43) 


Simi  larly , 


R9(r,k,  v  ) 
Z  m 


2  P(^h  sin  0 q 


A 


v  +  k 

ln(— —  )  + 


1  k- 1  Mv  +  r  k  +  2 
m 


ln( 


v  -k-  1 

-S2 _ ) 

v  +■  1  1 

m 


+  A 


v  tr 

3  n  ( — — r )  i 


2  r- 1  \  fl  r+  2 

m 


\  -k-  I 

ln(~v"  n  1 
m 


(2-44) 


(ii)  the  diagonal  terms,  i.e.,  k  r 
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Rl<k'k'  vm> 


sin2  A  1  v  -k  1  -  v  / t‘ 

- 5  — j  ln(— — — )  +  E— 

”Poh  -  2r  vm+k 


Z  2 
v  -  r 
m 


(2-45) 


and 


RJk,k,v„  ) 


m 


2PQh  sinZ  0Q 
77 


-1 


(2k+  1 )  (k— 1 )  v  +k 

m 


(2-46) 


( 2k+ 1 )  (k+2)  (v  -k-  1 ) 


Having  found  the  error  estimate  due  to  the  truncation  of  the  infinite 

sum  in  the  calculation  of  each  matrix  element,  it  has  been  possible  to  get 

the  relative  error  introduced  in  each  element  by  comparing  the  error 

estimates  with  the  corresponding  truncated  matrix  elements  and  hence 

keep  the  relative  error  in  each  element  of  the  matrix  less  than  a 

3 

specified  value  (one  part  in  10  in  our  calculation).  Furthermore,  by 
adding  the  error  estimates  to  each  element  of  the  matrix  and  solving  the 
resulting  equations,  it  is  possible  to  determine  how  sensitive  the  solution 
is  to  the  degree  of  uncertainty  in  the  matrix  elements. 

Before  solving  (1-90),  it  is  of  interest  to  consider  an  approximate 
solution  derived  originally  by  Schelkunoff  [6]  by  assuming  a  sinusoidal 
current  distribution  along  the  thin  biconical  antenna  and  later  from  a 
different  point  of  view  by  Tai  [4]  by  neglecting  all  the  off-diagonal 
elements  of  the  coefficient  matrix  of  (1-90)  under  the  assumption  that  the 
half  cone  angle  0 ^  is  small  so  that  the  diagonal  elements  are  dominant. 

The  fact  that  the  same  formula  was  obtained  by  both  for  d^  in  (1-90) 
signifies  nothing  more  than  that  the  approximations  involved  in  the  two 
approaches  are  essentially  the  same.  In  other  words,  a  sinusoidal 
current  distribution  in  itself  implies  sufficiently  severe  restrictions  on  the 
magnitude  of  the  cone  angle.  Ab  in  the  thin  cylindrical  antenna  [4],  a 
sinusoidal  current  distribution  can  only  be  accomplished  with  a  line  of 
ideal  generators  distributed  along  the  length  of  the  antenna  which,  in 
practice,  does  not  exist  and  thus  the  assumed  distribution  is  acceptable 
only  for  determining  certain  properties  of  very  thin  antennas  over  a  very 
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limited  range  of  length. 

The  approximate  formula  for  the  input  impedance  of  a  thin  bi  conical 
antenna  is  [4  ], 


ZC  sinB^h  -  iZ  cos  B^h 
y  ^  me _ r0  J  c _  0  ^ 

C  Z  sin(3  h  -  j  Z°  cosp  h 
c  r0  J  me  r0 


(2-47) 


whe  re 


me 


RC  4  j  Xe 

me  me 


(2-48) 


R  30  I  2Cin  2  Bnh  4  (Si  4  (3nh  -  2  Si  2  (3 nh )  sin  2  (3nh 

me  L  0  0  0  0 


4  ( 2Cin  2  (3^h  -  Cin 4  Pgh)  cos  2  (3^h  ] 


(2-49a) 


xv 


me 


30  [2Si2(3^h  4  (ln4  -  Cin4(3^h)  sin2(3^h 
-  Si  4(3^h*  cos  2f3^h] 


(2-49b) 


where  Si  x 


x  x 

sin  u  _j  •  ,1  .  ,  i  .  ,  1-cos  u  j 

- —  du  is  the  integral  sine  and  Gin  x  - au, 


0 


J0 


Generally  speaking,  the  coefficient  matrix  of  (1-90)  has  symmetric 
elements  with  its  dominant  part  lying  on  and  near  the  main  diagonal.  For 
each  row  in  the  upper  triangular  matrix  of  coefficients,  the  diagonal 
element  is  always  the  dominant  quantity  with  the  other  matrix  elements 
rapidly  decreasing  in  value  as  the  distance  from  the  diagonal  is  increased. 
The  size  of  the  diagonal  elements,  however,  also  decreases  in  the 
direction  of  larger  row  numbers,  but  with  a  much  smaller  rate  than  that 
of  the  off-diagonal  ones. 

It  is  instructive  to  look  into  the  elements  of  the  coefficient  matrix 
of  (1-90)  for  the  case  of  a  quarter  wave  conical  dipole.  Note  that,  while 
all  the  off-diagonal  elements  are  decreasing  with  shrinking  cone  angles, 
the  diagonal  elements  continue  to  increase  substantially.  Therefore,  in 
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the  limit,  this  justifies  the  approximation  made  by  Tai  in  deriving  (2-47) 
for  extremely  small  cone  angles  (9^  <  0.  01).  Whereas,  on  the  other 
hand,  when  the  angle  becomes  larger,  the  opposite  trend  occurs  and  the 
dominance  in  value  seems  to  spread  from  the  diagonal  elements  to  the 
ones  near  them.  Hence,  the  neglect  of  off-diagonal  terms  will  cause 
great  errors  in  the  final  solution.  A  typical  example  of  a  5  x  5  matrix, 
normalized  by  its  maximum  element  for  the  cone  angles  1  and  0.  1  ,  is 
shown  below 


1.  000 

0.  040 

0.  020 

0.  010 

0.  007 

0.  Z30 

0.  070 

0.  005 

0.  003 

[C] 

0.  096 

0.  003 

0.  002  ! 

V‘° 

0.  056 

0.  00  1  31 

- 

0.  034 

1.  000 

0.  010 

0.  005 

0.  003 

0.  002 

0.  160 

0.  018 

0.  001 

0.  0008 

[C]  - 

1 

0.  083 

0.  0006 

0.  0004 

V  °- 1 

1 

i 

0.  050 

0.  0003 

0.  0320 


A  concentration  of  the  dominant  values  is  readily  observed  from 
the  above  examples.  But  even  for  angles  as  small  as  0.  1°,  element  C^j 
which  is  small  compared  to  is  still  of  comparable  magnitude  with 

respect  to  element  C^.  Hence,  the  neglect  of  C  ^  will  cause  at  least 
10%  error  in  the  final  solution  for  by  For  large  cone  angles,  elements 
in  the  first  row,  and  hence  first  column,  become  quite  comparable  to  the 
corresponding  diagonal  elements,  which  implies  that  a  greater  error  is 
introduced  by  neglecting  the  off-diagonal  terms. 

The  results  presented  in  Fig.  2-3  are  the  calculated  input  impedance 
from  (1-90)  for  the  biconical  antenna  with  cone  angles  0  1.  5°,  1.  1° 

and  0.  1  .  The  impedance  variation  is  given  as  a  function  of  the  electric 
length  of  the  antenna  (i.  e.  ,  p^h)  with  particular  points  where  the  electric 
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length  attains  the  values  \/4  and  \/Z  specially  noted.  Larger  cone  angles 
tend  to  flatten  the  impedance  curve  and  shorten  both  the  resonant  and 
anti  resonant  lengths  of  the  antenna.  A  discussion  of  the  characteristics 
of  the  curves  for  the  impedance  1  j  X^  may  be  centered  about 

certain  critical  points.  Of  these,  the  most  important  are  the  values  of 
P^h  for  which  the  reactance  vanishes  and  which  involve  some  interesting 
points  which  have  counterparts  in  the  cylindrical  antenna. 

Input  resonance  is  characterized  by 


X, 


R0  (Ro’res. 


and  X, 


(X 


07  res. 


(2-50) 


The  condition  of  input  resonance  occurs  first  when  p^h  is  somewhat 
smaller  than  77 / Z .  It  is  seen  that  the  resonant  length  of  a  biconical 
antenna  moves  toward  77/2  as  l/0 ^  (and  hence  Z^)  increases.  It  is 
significant  to  note  that  when  p^h  77/2  the  input  impedance  given  by  the 
approximate  formula  (2-47)  is 

ZQ  73.  13  4  j  153.  7  ohms  (2-51) 

Thus  the  input  impedance  of  the  thin  biconical  antenna  as  given  by 
(2-47)  is  a  constant  independent  of  the  cone  angle  0  ^  when  h  is  a  quarter 
wave  in  length  and  0^  is  sufficiently  small. 

For  the  cylindrical  antenna,  such  constancy  obtains  in  the  modified 
zeroth  order  approximation  when  for  p^h  =  77/2,  i.  e.  , 


( Z  ) 

J(r  cylindrical 


73.  13  4  j  42.  5  ohms 


(2.  52) 


It  was  observed  by  King  [1]  that  the  curve  of  Z^  at  p^h  77 / 2 
approaches  the  value  of  73.  13  4  j  42.  5  ohms  as  2  ln(2h/a)  increases. 

The  same  phenomenon  is  observed  for  the  biconical  antenna  as  the 
corresponding  characteristic  impedance  Z^  increases.  This  is  shown  in 
Fig.  2-4  in  which  the  impedance  curves  tend  to  converge  to  the  value  of 
73.  13  4  j  153.  7  ohms  as  the  cone  angle  becomes  vanishingly  small.  This 
suggests  that  the  formula  (2-47)  may  be  no  better  than  a  modified  zeroth- 
order  approximation. 


FIG.  2-4  INPUT  IMPEDANCE  OF  X/4  THIN  CONICAL 
ANTENNA  VS.  CONE  ANGLE  90 
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Input  anti  resonance  is  characterized  by 


XA  0,  Rn  (Rn)  .  and  p.h  (P  h)  (2-53) 

0  0  0  anti  res.  0  0  anbres. 


It  occurs  at  values  of  p^h  that  are  less  than  77  by  an  amount  that 

increases  with  decreasing  in  a  manner  resembling  the  occurrence  of 

resonance  near  77/2.  The  reactance  has  a  maximum  for  p^h  slightly 

less  than  (P^h)  .  .  The  resistance  has  its  maximum  very  near  to 

r0  ant]  res. 

anti  -  resonance.  The  fact  that  the  ratio  jX  .  I  / 1  X  1  characterizing 

the  behavior  of  the  reactance  near  anti  resonance  is  always  greater  than 

one  for  the  cylindrical  antenna  is  not  necessarily  true  of  the  biconical 

antenna.  The  ratio  IX  .  1/  X  |  for  the  biconical  antenna  is  quite 

1  min'  1  max' 

different:  in  fact,  ,X  .  |  / 1  X  I  is  greater  than  one  only  for  small  cone 

angles  (0  ^  <  1  )  and  becomes  less  than  one  with  increasing  0^.  This 
implies  that  the  capacitive  lobe  of  the  reactance  curve  increases  with  0 
and  becomes  larger  than  the  associated  inductive  lobe.  In  other  words, 
the  capacitive  loading  effect  of  the  upper  and  lower  conical  surfaces  which 
act  more  or  less  as  a  capacitor  with  electric  lines  between  them  to 
the  point  generator  are  much  larger  in  the  biconical  case  than  in  the 
cylindrical  one. 

The  difference  in  the  behavior  of  the  reactance  curves  might  be 
ascribed  to  the  quite  different  driving  condition  in  these  two  types  of 
antenna;  for  the  biconical  antenna,  there  exists  a  point  generator  at  the 
center  while  for  the  cylindrical  dipole,  the  excitation  is  actually  over  a 
circular  surface  with  vanishing  separation  between  them.  This  suggestion 
is  later  justified  in  the  experimental  comparison  between  modified 
conical  and  cylindrical  dipoles  in  which  both  are  driven  by  a  circular 
radial  electric  field  at  some  distance  away  from  the  origin.  Both  real 
and  imaginary  parts  of  the  impedance  curves  are  more  comparable  than 
observed  here. 

A  comparison  of  the  input  impedance  of  the  thin  biconical  antenna 
as  computed  from  (1-90)  with  the  corresponding  approximate  results 
using  (2-47)  is  shown  in  Fig.  2-5.  It  is  seen  that  the  approximate  results 
lead  to  lower  resistance  near  anti  -  re  sonance  and  resonance  and  larger 


Ohms  Ohms 


(a)  e0sCU° 


(b)  9  --11° 
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resonant  and  anti- resonant  length  than  the  computed  results.  The 
reactance  behaves  in  an  analogous  manner.  In  genera],  the  approximate 
formula  leads  to  results  comparable  to  those  with  (1-90)  only  at  very  small 
cone  angles.  This  suggests  that  for  cone  angles  greater  than  0.  1°,  (1-90) 
must  be  used  for  studying  the  behavior  at  resonance  and  anti  -  resonance. 

In  Fig.  2-6,  is  shown  the  variation  of  the  cylindrical  antenna 
impedance  as  a  function  of  antenna  length  with  various  values  of  the 
parameter  SI  and  a  comparison  with  a  biconical  antenna  with  half  angle 
0 q  1.  1°.  The  comparison  is  interesting  since  by  judicious  adjustment 
of  the  radius  of  the  cylindrical  antenna,  a  measure  of  correspondence  can 
be  achieved.  It  is  obvious  from  Fig.  2-6  that  for  the  cylindrical  antenna 
with  radius  equal  to  or  greater  than  the  maximum  radius  of  the  biconical 
antenna  a^,  the  impedance  behavior  is  quite  different  from  that  of  the 
biconical  antenna,  whereas  the  impedance  curves  tend  to  be  comparable 
as  the  radius  of  cylindrical  antenna  remains  within  the  range  of  l/4  or 
l/5  of  a^.  Again  the  resistance  curves  seem  to  match  each  other  better 
than  the  reactance  curves.  One  way  of  measuring  the  comparability  is 
by  defining  the  average  of  the  antenna  as  follows  [6],  i.  e.  , 

T_ 

?  n  ?  R 

<  ft  >  ~  ln(— )  dR  (2-54) 

ave.  h  q  r  v 

where  r  is  the  radius  of  the  section  at  radial  distance  R  from  the  driving 
point  at  the  origin.  Accordingly,  the  average  ft  for  the  cylindrical 
antenna  with  radius  a  is 


<  ft  > 


ave. 


2  [In  (^)  -  1] 


(2-55) 


and  that  for  the  conical  antenna  with  angle  6 ^  is 


<  ft  > 


ave. 


2  ln(2  cot  q) 


(2-56) 
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It  is  seen  in  Fig.  2-6  that  the  impedances  of  the  bi  conical  and 
cylindrical  antennas  are  comparable  when  both  have  approximately  the 
same  >(for  curve  (d)  <f2>  12.  5  and  <$~2  '  conical  This 

means  that  for  purpose  of  comparison,  one  must  choose  a  cylindrical 
antenna  with  a  radius  at  least  smaller  than  half  the  maximum  radius  of 
the  cone  antenna.  No  attempt  has  been  made  here  to  pursue  further  a 
theoretical  study  of  the  relation  between  these  two  radiating  structures. 

The  general  properties  provide  a  qualitative  guide  of  the  extent  to  which 
these  two  radiating  antennas  behave  reasonably  alike,  so  that  information 
obtained  from  one  can  be  generally  applied  to  predict  the  behavior  of  the 
other.  This  is  essential  later  in  interpretating  the  expe ri mental 
comparison  between  modified  cylindrical  and  conical  antennas  in  which 
both  the  cone  angle  and  antenna  radius  are  fixed  by  the  measuring  set-up. 

Numerical  values  for  (1-97)  and  (1-99)  for  a  corresponding  biconical 

structure  are  also  obtained  along  much  the  same  lines  as  for  (1-90).  All 

three  sets  of  equations  give  the  same  results  with  differences  less  than  1 

per  cent  among  them.  For  each  equation,  matrices  of  order  N-2,  N-l, 

N  etc.  are  solved  separately  and  the  solution  inspected  for  convergence. 

It  was  found  if  v  >  max.  (2k-f*l,  2r-fl)  that  the  solutions  converge  both 
m 

rapidly  and  monotonically .  However,  if  v  <  max.  (2k+l,  2r+l),  the 
solutions  tend  to  oscillate  widely  and  show  little  tendency  to  converge. 

The  choice  of  the  size  of  the  matrices  depends  upon  the  accuracy 
required  for  the  solution  on  the  one  hand  and  the  computation  time  on  the 
other.  As  a  compromise,  a  less  elegant  method  was  adopted.  That  is, 
the  extreme  cases  are  examined  with  great  care  (i.  e.  ,  the  longest  antenna). 
In  addition  to  checking  among  themselves  for  convergence  in  the  course  of 
a  numerical  solution,  solutions  from  a  truncated  matrix  were  continously 
checked  with  those  of  the  modified  matrix  in  which  each  element  was 
corrected  by  adding  the  error  estimates.  To  proceed  one  step  further, 
solutions  from  a  matrix  of  size  N  (which  is  the  number  when  solutions 
show  convergence  of  1  per  cent)  are  compared  with  solutions  from  a 
matrix  of  size  (N+10).  It  is  required  that  their  difference  lie  within  1  per 
cent.  Although  it  is  not  known  exactly  how  this  truncated  solution 
converges  to  the  true  solution,  at  least  a  self-consistent  one  is  obtained 
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in  this  manner.  The  number  N  is  then  used  for  antennas  of  shorter 
lengths  without  repeating  the  same  test  for  convergence. 

As  a  double  check  for  convergence,  obs e rvati ons  were  made' 
particularly  at  (3()h  1.  6  and  0^  1.  1°.  Numerical  results  from  (1-90) 

and  (1-99)  for  a  series  of  matrix  sizes  are  shown  in  Fig.  2-7.  It  is 
observed  that  values  of  (1-90)  tend  to  converge  toward  the  final  value  in 
a  monotoni cally  increasing  way  with  an  increasing  size  of  the  matrices 
while  the  solutions  from  (1-99)  also  approach  the  same  value  in  much  the 
same  but  monotoni  cally  decreasing  way. 

The  mutual  checking  among  these  three  sets  of  linear  equations 
suggests  that  the  procedures  used  in  the  derivation  of  (1-97)  and  (1-99) 
for  the  biconi  cal  antenna  and  hence,  in  quite  a  similar  manner  for  (1-57) 
and  (1-64)  for  the  modified  dipole  are  justified. 


5.  THE  MODIFIED  CONICAL  DIPOLE 


Based  upon  (1-57),  (1-64),  (1-73)  and  (1-76),  a  numerical  solution 
for  the  current  distributions  along  both  antenna  and  sphere  with  various 
sizes  can  be  obtained.  Also  the  input  conductance  and  susceptance  of  the 
antenna  are  readily  known  as  the  real  and  imaginary  parts  of  the  total 
current  at  the  driving  point.  In  the  process  of  the  numerical  computation, 
much  of  the  analysis  and  many  of  the  results  from  the  preceding  section 
can  readily  be  utilized  with  some  necessary  modifications. 

The  general  terms  of  the  coefficient  matrix  of  (1-57)  and  (1-64)  are 
the  following  forms: 


Kj (r,  k,  v) 


(Zv+1)  MM(3Qa)  1 

v(v+l)  M v ( P q a )  [k(k+ 1)  -  v(v+ 1 )] 


_1 _  3v 

[r(  r-t  1 )  -  v(v  +  1 )  ]  3|i^ 


(2-57) 


and 
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FIG.  2-7  CONVERGENCE  OF  NUMERICAL  SOLUTION  FOR 
X/4  BICONICAL  ANTENNA  (0O=1.1°) 


-33- 


M  (p  a)  2 

g  (r,k,v)  g  (r,k,v)  — - -  (2-58) 

_  M;(PQa) 

where  subscripts  1  and  2  denote  the  first  and  second  sets  of  linear 
equations  respectively. 

The  above  equations  differ  from  those  for  the  co r  res pondi  ng 
bi  conical  antenna  only  in  the  ratio  MM(3^a)/M ^  ((3^a) . 

By  applying  the  approximate  formulae  for  the  Bessel’s  function  of 
the  second  kind  for  large  v,  i.  e.  y 


Nv(x) 


(7T) 


1/2 


2V  r  (v  +  1/2)  x“v  exp  4(v^l/2) 


(2-59a) 


and 


N'  (x)  ^ 

v 


V 

X 


N  (x) 
v 


(2- 59b) 


and  by  using  the  fact  that  v  »  x,  one  obtains 


M'(x)  (v+1)  P(x) 

Mv(x)  ^  x  Jv.(x) 


(2-60) 


The  fact  that  h/P  (x)/M^  (x)  and  J^(x)/J  (  x)  have  the  same  asymptotic 

behavior  for  large  v  implies  that  all  of  the  formulae  for  error  estimates 

in  the  preceding  section  [(2-35)  -  (2-46)]  can  be  used  directly  for  the 

modified  conical  antenna  with  (3^h  replaced  by  p^a.  Furthermore,  the 

same  error  bound  obtained  for  the  two  cases  suggests  that  the  coefficient 

matrix  associated  with  the  modified  dipole  should  have  essentially  the 

same  convergence  properties  as  that  for  the  biconical  one;  especially  for 

a  small  central  sphere  in  which  D  ~  and  D’  ~  are  small  relative  to  D  t 
r  r2  r2  rl 

and  D  ^  [(1-57)  and  (1-64)]. 

Additional  terms  (i.  e.  ,  D  ~  and  DT  »)  which  characterize  the 

’  r2  r2 

presence  of  the  central  sphere  for  the  modified  dipole  are  also  present  in 
the  form  of  an  infinite  summation  in  the  driving  terms  of  (1-57)  and  (1-64). 
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For  a  small  central  sphere,  these  terms  are  of  the  order  (P{)b)  [where 
v  is  the  first  root  of  the  characteristic  equation  I jv(^q)  0)  and  vanish  in 

the  limiting  case  of  the  biconical  antenna.  On  the  other  hand,  ?  and 
D*  ^  become  more  important  as  the  central  sphere  becomes  larger. 


The  error  estimates  for  D  ^  and  D!  0  will  be  examined  first.  Thus, 

rZ  r  Z 


D 


rZ 


-JV0  2 

60  Pr("o)(1  "'V 


V 

V 


c 

v 


(2v+  1 ) 
v(v+  1  ) 


1 

[ r(r+ 1 )  -  v(v+  1 )  I 


_ 1 _  3v 

Mv(Poa^  Jv^O*3^  ^0 

again,  for  large  v,  one  gets 

-  s  i  n  ( v  0  ) 

_ _ u _ ^ _ 

v  fn(0Q,  0  j )  sin  sin(0o  +  0  j ) 1 


(2-61) 


(2-62a) 


Mv((30a)  *  (P0b)  ^  2vTT  <!)V  (2'62b) 

Therefore,  the  upper  bound  of  the  remainder  of  D  ^  is 
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ID 


(r. 


v  ) 
m 


4 _ 1 _  jb^Vm  1 

77  fn(0Q,0j)  a  [sin  0Q-  sin(0o  +  0  j  )]^  2 


1  v  1  v 

_  1  i  in  '  ,  _  ,  ,  m  . 

r ( 2 r+  1 )  V  +  r  (2r+l)(r+l)  v  - 2 r-  l1 
m  m 

(2-63) 


r  2  ' 


Similarly  for  D 
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R2D(r’  Vm) 


4a  m 

Jl  a 


fn(0o,  ^  i)[sin  0q-  sin(0Q  +  01)] 
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1  1 
-  +  - 


(r+l)v  (r4  1)  v 

m  m 


v-r-1  1  v  4  r 

i  /  m  ,  t  ,  m  , 

2  Inf - )  -  —  ln( - ) 


m 


(2-64) 


The  errors  caused  by  the  double  truncation  which  occurs  in  the 
computation  of  the  fields  fi.  e.  ,  anc*  (E^)^]  and  the  input  admittance' 

(i.  e.  ,  Yir^)  can  be  estimated  in  the  following  way.  A  comparison  of  the 
real  part  of  the  complex  power  radiated  as  seen  from  the  driving  point 
with  the  far-field  radiated  power  computed  from  (1-89)  gives: 
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0 


240  S' 
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(2k+  l)(k+l)  k 


1 

|Rk(P0a>l2 


(2-65) 


The  right-hand  side  converges  extremely  fast.  It  was  found  that 

5 

less  than  ten  terms  were  needed  to  achieve  accuracy  of  one  out  of  10  . 

It  is  reassuring  to  note  that  in  the  calculation,  the  left-hand  side  of  (2-65) 

3 

agrees  within  the  range  of  relative  error  one  part  out  of  10  with  the 
right-hand  side. 

It  is  interesting  to  note  that  the  computation  on  the  right-hand  side  of 
(2-65)  involves  only  the  radiation  field  coefficients  (i.  e.  ,  b^’  s)  which  is 
the  solution  of  an  infinite  set  of  linear  equations  while  the  left-hand  side 
is  a  result  involving  all  the  fields  inside  antenna  region.  The  excellent 
agreement  between  these  two  demonstrates  the  internal  consistency  and 
accuracy  of  the  computation. 


(A)  Input  Admittance 

The  input  admittances  of  modified  conical  antennas  with  half  cone 
angle  0 ^  1.1°  and  driving  gap  1.5°  and  antenna  heights  ranging  from 
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p  h  0.6  to  4.0  and  spherical  radius  ranging  from  p^b  0.  1  5  to  1.5  have* 
been  obtained  (numerical  data  are  listed  in  Appendix  A). 

In  Fig.  2-8,  are  shown  graphs  of  the  input  admittance  of  a  modified 
conical  antenna  as  a  function  of  p^h  with  the  sphere  radius  as  a  parameter. 
This  is  convenient  theoretically  and  also  for  a  later  experimental  check 
in  conjunction  with  a  partially  fixed  installation  for  which  the  frequency 
is  varied  (i.  e.  ,  in  the  experimental  set-up,  the  physical  size  of  the 
central  sphere,  center  conductor  and  hence  the  cone  angle  are  fixed  while 
the  antenna  length  varies  in  discrete  steps). 

Because  the  input  admittance  of  the  modified  conical  antenna  is 
defined  as  that  seen  at  one  driving  terminal,  the  admittance  curves  in 
Fig.  2-8  have  to  be  rescaled  by  a  factor  1  / 2  in  order  to  compare  them 
with  that  for  a  corresponding  biconical  antenna. 

In  general,  the  admittance  curves  are  broadened  and  flattened  as 
the  center  sphere  is  enlarged.  More  specifically,  the  magnitude  of  the 

maximum  values  of  ,G  i  and  B  j  and  the  minimum  B  . 

1  max1  1  max'  1  min1 

decrease  as  the  radius  of  the  centra]  sphere  increases  until  approximately 
at  p^b  1  where  the  circumference  of  the  sphere  equals  one  wavelength, 
and  the  magnitudes  of  iG  L  [B  L  and  IB  .  I  start  growing  again 
with  increases  in  the  size  of  the  sphere. 


The  behavior  of  the  susceptance  near  resonance  and  anti  -  res onance 

is  also  of  interest  with  a  change  in  the  radius  b.  In  Fig.  2-9,  are  shown 

the  resonant  and  anti  -  resonant  lengths  as  a  function  of  p^b.  It  is  seen 

that  (p„h)  . .  first  decreases  with  increases  in  P~b  until  P^b 

r0  antires.  r0  ru 


approximately  equals  one,  when  (pnh) 

On  the  other  hand,  (p„h) 

r0  res. 


Au/  A.  begins  to  increase  again. 

0  antires.  H  H 

remains  at  much  the  same  value  for 


p  b  <  0.  5,  after  which  (Pnh) 


0  res. 


increases  rapidly.  The  increasing  trend 


of  ( P qH )  r e s  1S  stopped  and  reserved  also  at  approximately  p^b  1.  The 

relative  constancy  of  (p^h)  for  p„b  <  0.  5  might  be  the  result  of 

simultaneous  capacity  loading  due  to  a  central  modification  which  tends 

to  shift  (P^h)  g  toward  the  right  and  a  thickening  on  the  average  of  the 

protruding  portion  which  tends  to  pull  (P_h) 

(J  res . 

values. 


back  toward  smaller 
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ln  the  neighborhood  of  (Pgh)res  the  slope  of  the  susceptance  curves 
decreases  with  increases  in  p^b.  The  flattening  of  the  curve  is  accompanied 
by  a  shift  of  the  point  of  zero  susceptance  (p^b  >  0.  5)  from  smaller  values 
toward  tt/2.  The  fact  that  (p Qh) r es  becomes  even  larger  than  77 / Z  as 
(p^b)  becomes  greater  than  0.  8  suggests  that  the  central  modification  it¬ 
self  becomes  such  a  substantial  part  of  the  whole  radiating  structure  that 
the  performance  of  the  antenna  protruding  out  from  the  sphere  is  changed 
beyond  any  resemblance  with  the  corresponding  dipole. 

In  the  neighborhood  of  p^b  1,  the  susceptance  stays  fairly  constant 
over  the  range  of  p^h  from  77 / Z  to  7 7  even  though  the  associated  conduct¬ 
ances  differ  significantly  (e.  g.  ,  at  p^h  0.  91,  the  susceptance  varies 
within  0.  7  millimhos  of  a  mean  value  near  0.  2  millimhos,  whereas  the 
conductance  ranges  between  2  and  8  millimhos).  This  observation  con¬ 
firms  the  experimental  results  for  the  modified  cylindrical  antenna 
reported  by  Iizuka  [8], 

Consider  specifically  the  effect  on  the  input  characteristics  of  a 

sudden  change  from  a  point  generator  in  the  biconical  case  to  a  radial 

ring  driving  a  modified  conical  antenna.  In  Fig.  2-10,  the  input  impedance 

is  shown  near  (B^h)  _  for  both  the  biconical  and  the  modified 

r0  antires. 

conical  antenna  with  a  small  central  modification  (p^b  0.  15).  Although 
the  real  parts  are  still  comparable  with  each  other,  one  remarkable 
change  in  the  imaginary  parts  is  the  ratio  (X  .  /  X  which  is 

larger  than  one  in  the  biconical  case  but  is  now  less  than  one.  This  means 
that  with  a  small  central  modification  which  provides  the  modified 
conical  antenna  with  a  driving  condition  quite  similar  to  that  of  the 
cylindrical  antenna,  the  input  characteristics  of  the  two  become  more 
comparable. 

In  Fig.  2-11,  is  shown  the  variation  of  antenna  admittance  with  a 
fixed  length  of  protrusion  as  a  function  of  the  radius  of  the  central  sphere. 
Again,  the  driving  point  admittance  of  the  antenna  is  seen  to  change 
significantly  with  the  size  of  the  sphere.  Generally  speaking,  for 
P^b  <  1,  the  change  in  susceptance  with  respect  to  p^b  is  larger  than  that 
of  the  conductance.  Specifically,  for  p^b  <  0.  2,  there  are  practically  no 


FIG.  2-10  INPUT  IMPEDANCE  NEAR  ANTI  RESONANCE  OF 

CONICAL  ANTENNA  WITH  SMALL  CENTRAL  SPHERE 


Mi  llimhos  Miiiimhos 


FIG.  2-11  VARIATION  OF  INPUT  ADMITTANCE  FOR  MODIFIED 
CONICAL  ANTENNA  WITH  FIXED  LENGTH 
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changes  in  the  conductance  but  relatively  large  changes  in  the  susceptance. 
It  is  also  significant  to  note  that  all  admittance  curves  can  be  extra¬ 
polated  toward  smaller  (3^b  and  that  they  coincide  with  the  values  of  the 
corresponding  biconical  antenna. 

In  Fig.  2-12,  are  shown  the  changes  in  antenna  admittance  with 
fixed  total  length  (i.  e.  ,  (3^a  is  fixed)  of  the  radiating  structure  with  the 
radius  of  the  central  sphere  (3^b  as  the  parameter.  It  is  seen  that  drastic 
changes  in  performance  are  introduced  by  the  central  sphere.  Again,  as 
was  observed  in  Fig.  2-11,  the  modified  dipole  with  a  small  central 
sphere  behaves  essentially  as  the  corresponding  dipole  with  the  same 
length. 

With  its  driving  point  situated  symmetrically  at  a  distance  from  the 
center  and  an  enlarged  central  volume,  the  modified  dipole  bears  some 
resemblances  to  the  sleeve  dipole.  One  of  the  important  properties  of 
the  sleeve  dipole  is  its  broad-band  behavior.  One  would  expect  the  same 
behavior  from  the  modified  dipole.  This  subject  will  not  be  pursued 
extensively  here,  since  it  is  not  involved  in  a  fundamental  check  of  the 
validity  of  the  theoretical  formulae.  Furthermore,  the  broad-band 
properties  cannot  be  investigated  in  general  without  extensive  numerical 
computations  of  the  input  impedance  for  a  range  of  antenna  lengths  and 
sphere  radii.  However,  the  fact  that  the  modified  conical  dipole  does 
possess  broad-band  properties  can  easily  be  seen  by  using  the  available 
tabulation  of  theoretical  input  impedance  (Appendix  A).  Suppose  the 
frequency  for  which  (3^h  1.  6,  |3^(h42b)  3.42  and  p^b  0.  91  is  fQ  and 

the  associated  angular  frequency  is  cj^.  For  a  variation  in  frequency  from 
Uq  -  AgJq  to  Uq+AgJq,  where  cj^/4,  the  electrical  lengths  have  the 

following  ranges: 

1.2  ^  (30h  ^  2.0 

luo€uo$luo  2-6  ^  (30(h+2b)  ^4-2 

'  0.  7  ^  (30b  <  1.  1 


Mill'imhos  Millimhos 


FIG.  2-12  VARIATION  OF  INPUT  ADMITTANCE  FOR  MODIFIED 
CONICAL  ANTENNA  WITH  FIXED  TOTAL  LENGTH. 
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Let  the  frequency  response  of  the  modified  conical  dipole  be 
compared  with  corresponding  biconical  dipole  with  half  lengths  P^h^  Pt)h 
or  P(jh^  P^(h+2b).  Let  Z^  and  denote  the  input  impedance  for 
antennas  with  half  lengths  p^h^  and  p^h^  respectively.  Tables  (A-l) 
through  (A-12)  show  the  input  impedance  of  the  modified  dipole  and  Fig. 

2-3  shows  the  input  impedance  of  the  biconical  antenna.  It  is  seen  that 
the  reactance  of  the  modified  conical  dipole  is  fairly  constant  near  a  mean 
value  of  -4  ohms  over  the  range  (i.  e.  ,  from  25.  3  ohms  at  P(^h  2.  0  to 
-38.  82  ohms  at  P^h  1.  6),  whereas  the  corresponding  reactance  of  the 
center-driven  conical  dipole  (Fig.  2-3)  varies  widely  over  positive  and 
negative  values  (20  <  <  480,  -675  <  <  650).  The  resistance  for  the 

modified  dipole  varies  within  30  per  cent  of  a  mean  value  near  111  ohms 
(77  <  Rq  <  145)  whereas  ranges  between  50  ohms  and  475  ohms  and  R^ 
varies  over  a  range  between  175  ohms  and  1550  ohms.  Clearly  the 
frequency  response  of  the  modified  conical  dipole  is  very  much  better 
from  the  point  of  view  of  broad-band  operation  than  the  response  of  an 
unmodified  dipole. 

(B)  Current  Distribution 

Due  to  the  continuous  change  in  the  cross  section  of  the  conical 
antenna,  the  current  density  along  the  conical  surface  is  proportional  to 
the  total  current  crossing  the  corresponding  circular  cross  section  with  a 
constant  whose  value  changes  with  the  antenna  length,  i.  e.  ,  I  =  2ttR  sin0x 
1^  where  b  ^  R  ^  a.  Correspondingly,  for  the  central  sphere  one  can 
also  define  I^g(0)  27fb  sin0  1^  (0)  where  4  0 ^  ^9  ^  77  -  0^  -  9 The 
current  density  can  be  compared  directly  with  experimental  data  while 
the  total  current  along  the  antenna  has  the  advantage  of  being  readily 
compared  with  its  corresponding  cylindrical  dipole  as  long  as  the  cone 
angle  is  small. 

Numerical  results  for  the  current  distribution  of  antennas  with  the 
half  cone  angle  9 ^  1.  1°,  antenna  lengths  p^h  1.  6  and  3.  1  and  central 

sphere  with  p^b  0,  0.  31,  0.  51,  0.  91  and  1.  31  have  been  obtained  and 
shown  in  Fig.  2-13  and  Fig.  2-14. 
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la)  pQh  1.6 

Both  total  current  distributions  and  the  corresponding  current 
density  for  p^h  1.  6  with  various  central  spheres  are  shown  in  Fig.  2-13. 
The  curves  show  components  1^/ Vq,  I^/Vq  as  a  function  of  [  p^(R-b)] 

and  1  q  ,  aS  a  ^unct^on  Pgb9)  with  corresponding  angle 

variations  marked.  It  is  to  be  noted  that  the  current  densities  presented 
here  are  in  milli amperes  per  volt  per  electrical  length,  so  they  differ 
from  actual  current  density  by  a  constant  p  U  pf. 

Both  the  real  and  imaginary  part  of  the  total  current  I  remain 
finite  at  the  end  and  change  their  direction  as  the  current  passes  around 
the  edge  of  the  conical  antenna  to  the  cap.  The  real  part  of  the  total 
current  seems  less  affected  by  the  size  of  the  central  sphere  and  remains 
much  the  same  in  shape  except  for  a  continuous  increase  in  magnitude 
with  an  increase  in  the  size  of  the  central  sphere.  The  real  current 
behaves  very  much  like  a  combination  of  terms  made  up  of  cosines  and 
shifted  cosines  with  half-angle  arguments.  The  imaginary  part  of  the 
current  1^  is  capacitive  along  the  antenna  for  p^b  <  0.  91  but  the  magnitude 
continues  to  decrease  with  an  increase  in  the  size  of  central  sphere.  The 
range  over  which  the  IJ  is  capacitive  is  reduced  and  inductive  currents 
begin  to  appear  near  the  driving  point.  Although  the  total  current  on  the 
surface  of  the  sphere  remains  quite  comparable  in  magnitude  with  that  on 
the  antenna,  the  magnitude  of  the  corresponding  current  density  decreases 
sharply  writh  distance  from  the  driving  point.  This  is  due  to  the  fact  that 
on  the  hemisphere,  the  surface  of  an  incremental  belt  (277b  sin@)  A  9 
increases  with  an  increase  in  6  and  current  spreads  out  over  a  widening 
area  with  distance  away  from  the  driving  point. 

The  real  current  IM  ,  wrhich  has  the  same  value  as  IM  at  the  driving 

ts  9  tc  6 

point,  remains  essentially  constant  for  a  small  sphere  (i.  e.  ,  P^b  <  0.  31). 

As  the  sphere  becomes  larger,  it  can  be  approximated  by  trigonometric 
functions.  The  imaginary  part  I|  like  its  counterpart  IJ  changes 
significantly  with  variations  in  sphere  radius.  V  becomes  larger  in 
magnitude  and  more  capacitive  with  an  increase  in  sphere  radius. 


FIG.  2-13  CURRENT  DISTRIBUTIONS  OF  MODIFIED  CONICAL  ANTENNA. 
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lb)  (3Qh  3.  1 

The  distributions  of  total  current  for  the  antenna  with  (3^h  =  3.  1 

are  shown  in  Fig.  2-14.  The  general  shapes  of  the  current  distribution 

along  the  conical  antenna  seem  less  affected  by  the  presence  of  the  central 

sphere  than  when  (3^h  1.  6.  In  general,  the  real  current  1M  can  still  be 

^  r0  ’  tc 

app r oxi  mated  by  a  shifted  cosine  and  a  shifted  cosine  with  half-angle 
argument,  while  the  imaginary  current  1^  can  still  be  approximated  by 
trigonometric  functions.  Near  the  end  of  the  antenna,  the  real  part  of  the 
current  is  negative  for  small  spheres.  As  b  increases,  this  region  with 
a  small  negative  value  decreases  and  the  real  current  becomes  more  and 
more  positive.  The  imaginary  currently  remains  essentially  capacitive 
over  most  of  its  range  except  near  the  driving  point.  The  range  over 
which  the  current  is  capacitive  remains  fairly  constant  except  for  a 
continuous  increase  in  magnitude  as  the  central  sphere  is  made  larger. 

The  current  distribution  on  the  sphere  seems  to  have  undergone  a 

drastic  change  when  compared  with  that  for  p^h  1.  6.  Again,  the 

behavior  of  I!1  bears  great  resemblance  to  that  with  S^h  1.  6  but  the  rate 
ts  6 

of  change  with  respect  to  the  sphere  radius  seems  to  be  approximately 
doubled  (e.  g.  ,  V'  begins  to  have  a  negative  part  for  p^b  >  0.  91  at 
P^h  1.  6,  whereas  F^  has  a  negative  part  for  a  much  smaller  central 
sphere).  The  most  affected  part  is  the  imaginary  current  F^  which  has  a 
small  range  of  inductive  current  near  the  driving  point  and  becomes 
capacitive  rapidly  over  the  rest  of  its  range  with  an  increase  in  the  size 
of  the  central  sphere.  At  Pf)b  1  ,  begins  to  become  less  capacitive 
near  the  equatorial  region  of  the  sphere  (i.  e.  ,  9  7 t/2).  As  can  be  seen 

in  Fig.  2-  14c,  ll  becomes  inductive  again  near  the  9  tt/2. 

In  general,  the  distribution  of  current  along  the  thin  conical  antenna 
has  much  the  same  shape  as  along  a  simple  dipole  with  enlarging  central 
sphere.  Yet,  the  current  on  the  sphere  is  quite  sensitive  to  the  antenna 
lengths  especially  for  larger  spheres. 

(C)  Far-Field  Pattern 

Far-field  patterns  with  a  fixed  central  sphere  and  various  antenna 
length  are  shown  in  Fig.  2-15.  The  size  of  the  central  sphere  ranges  from 
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FIG.  2-14  CURRENT  DISTRIBUTIONS  MODIFIED  CONICAL  ANTENNA 


FIG.  2-15  FAR- FIELD  PATTERN  OF  MODIFIED  CONICAL  ANTENNA 
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(3^b  0.  15  to  1.51  and  the  antenna  lengths  range  from  (3^h  =  1.  6  to  3.  9. 

Because  of  the  properties  of  symmetry  of  this  radiating  structure 
(i.  e.  ,  rotational  symmetry  and  symmetry  with  respect  to  Q  -  tj  /  2  plane), 
it  is  necessary  only  to  know  the  field  patterns  in  the  range  0°  ^  0  ^90°. 

In  order  to  facilitate  the  observation  of  the  effect  of  the  central  sphere, 
all  curves  are  normalized  to  have  a  maximum  value  one. 

In  general,  the  far-field  pattern  of  the  modified  conical  antenna 
behaves  quite  like  a  simple  dipole  for  small  central  spheres.  The 
presence  of  the  central  sphere  effectively  lengthens  the  antenna  as 
compared  to  a  simple  dipole  of  length  equal  to  the  portion  that  protrudes 
from  the  sphere. 

Note  that  the  far-field  patterns  are  more  sensitive  to  the  total 
length  of  the  radiating  structure  than  to  the  size  of  the  central  sphere. 

In  Fig.  2- 15a,  are  drawn  the  field  patterns  of  a  quarter  wave 
antenna  with  various  central  spheres.  The  effect  of  the  central  sphere 
can  readily  be  seen  in  the  continuous  flattening  of  the  field  patterns  which 
correspond  in  effect  to  those  of  a  dipole  with  half  length  longer  than 
quarter  wavelength. 
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APPENDIX  A 

TABULATION  OF  THEORETICAL  INPUT  IMPEDANCES 

AND  ADMITTANCES 

The  driving  point  characteristics  (i.  e.  y  input  impedances  and 
admittances)  have  been  computed  for  modified  conical  antennas  with 
half  cone  angle  0 ^  1.  1°  and  driving  gap  0^  1.  5°.  It  covers  the 

following  ranges: 

(3Qb  0.  15  to  1.  51 

and 

(3^h  0.  6  to  3.  9 
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^able  A-l  'f'heoreti  cal  data  for  the  inout  adnittance  of 
modified  dirole; 


poh 

P0a 

Yo 

Z0 

0.60 

0.75 

0.126 

+.i 

3.850 

8.48 

-1259-33 

0.70 

0.85 

0.254 

+.1 

4.833 

10.36 

-.1206.34 

0.60 

0.95 

0.574 

+1 

6.313 

14.28 

-ji57.00 

0 .  °0 

1.05 

1.303 

+1 

8.394 

18. 05 

-1116.00 

1.00 

1.15 

3.298 

+Jii  .637 

22.54 

-j  79.54 

1.10 

1.25 

9.806 

+115-930 

27.90 

-j  45.50 

1.20 

1.35 

25.410 

+1 

0.732 

34.30 

-j  13.14 

1.30 

1 .45 

20.490 

-  1 

8.294 

41.93 

+1  16.97 

1 .40 

1.55 

10.210 

_  a 

t J 

9.713 

51 .40 

+  1  48.92 

1.50 

I.65 

6.103 

-1 

7.740 

62.82 

+  1  79.67 

1.60 

1.75 

4. 228 

-j 

6 . 097 

76.80 

+1110.75 

1.70 

1.85 

3.228 

-j 

4.884 

94.19 

+  1142 .50 

1.80 

1.95 

2.630 

-j 

3.971 

115.90 

+  .1175-00 

1.90 

2.05 

2.270 

-j 

3.277 

142.85 

+i206.23 

2.00 

2.15 

I.076 

-1 

2.675 

178.70 

+1241.80 

2.10 

2.25 

1 .805 

-.1 

2.205 

222.32 

+1271.60 

2.20 

2.35 

1 . 644 

-.1 

1.764 

282.70 

+J303.UO 

2.30 

2.45 

1.553 

-1 

1 .407 

353.73 

+1320.30 

2.40 

2.55 

1.457 

-j 

1.050 

452.20 

+J325.30 

2.50 

2.65 

1 ,40P 

-.1 

0.750 

553.19 

+j294.38 

2.60 

2.75 

1.350 

-  j 

0.431 

672.10 

+1214.60 

2.70 

2.85 

1 .318 

-  j 

0.141 

750.00 

+1  80.10 

2.80 

2.95 

1.298 

+  1 

0.145 

760.80 

-j  85.20 

2 .  °0 

3. 05 

1 .2Q0 

+.1 

0.433 

696.80 

-.1233.70 

3.00 

3.15 

I.204 

+.1 

0.727 

587.50 

-1329.70 

3.10 

3.25 

1.313 

+1 

1.033 

470.40 

-1370.00 

3.20 

3-35 

1.350 

+1 

1.359 

367. °o 

-1370.40 

3.30 

3.45 

1 .410 

+  .i 

1.713 

266.50 

-1348.00 

3.40 

3-55 

1.500 

+i 

2.061 

230.81 

-1317.24 

3.50 

3.65 

1.642 

+.1 

2 . 540 

173.60 

-1277.30 

3.60 

3.75 

1.837 

+  i 

u 

2.999 

148.56 

-1242.48 

3.70 

3.85 

2.186 

+1 

3.660 

120.30 

-1201 .40 

3.80 

3.95 

2.658 

+,i 

4.274 

104.92 

-1168.74 

3.90 

4.05 

3.605 

+.1 

5.159 

C1  .03 

-1130.25 
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TableA-2  Theoretical  data  for  the  input  admittance  of 
modified  dipole; 


Poh 

P0a 

Y0 

z 

0 

0.60 

0.80 

0.148 

+.i 

3.972 

9.36 

—  j  251 .40 

C.70 

0.90 

0.310 

+.1 

5,010 

12.31 

-il98.85 

0.0 

1.00 

0.705 

+3 

6.555 

16.20 

-jl50.80 

0.°0 

1  .10 

1.577 

+3 

8.671 

20.30 

-jlll .60 

1  .00 

1  .20 

3.920 

+jii.850 

25.15 

-j  76.00 

1  .10 

1  .30 

11.056 

+  .U5.340 

30.92 

-j  42.90 

1  .20 

1  .40 

24.280 

+3 

7.262 

37.81 

-j  11.31 

1  .30 

1 .50 

18.440 

-•3 

7.750 

46.09 

+j  19.38 

1 .40 

1.60 

9.990 

**  J 

8.849 

56.92 

+3  49.68 

l  .50 

1  .70 

6.16° 

-•.i 

7.232 

68.27 

+  3  80. 03 

1.60 

1 .80 

4.336 

-•j 

5.773 

83-19 

+JH0.75 

1 .70 

1 .00 

3.331 

*•3 

4.653 

101 .70 

+ jl42 . 10 

1.80 

2.00 

2.719 

--j 

3.792 

124.90 

+jl74.l6 

1  .?0 

2.10 

2.356 

-j 

3.H7 

154.30 

+ j204. 19 

2.00 

2.20 

2.039 

-j 

2.546 

191.60 

+j239.30 

2.10 

2.30 

1.867 

-.1 

2.082 

238.73 

+ j266 .19 

2.20 

2.40 

1 .687 

-j 

1.657 

301 .70 

+J296.30 

2.30 

2.50 

1.598 

-j 

1 .302 

376.19 

+3306.33 

2.40 

2.60 

1.485 

-j 

0.951 

477.60 

+3305.60 

2.50 

2.70 

1 .440 

-  j 

0.653 

575.85 

+ j26l .34 

2.60 

2.80 

1.367 

-3 

0.338 

689.70 

+3170.40 

2.70 

2.90 

1.329 

-«i 

0.048 

751 .40 

+3  26.96 

2.P0 

3.00 

1.304 

+J 

0.239 

741 .00 

-3135.85 

2.Q0 

3.10 

l  .292 

+J 

0.527 

663.60 

-3270.95 

3.00 

3.20 

1 .2°2 

+3 

0.823 

550.50 

-3350.80 

3.10 

3.30 

1.307 

+3 

1.133 

436.80 

-3378.70 

3.20 

3.40 

1.340 

+.1 

1.463 

340.34 

-3371.60 

3.30 

3.50 

1.393 

+3 

l  .823 

264.85 

-3345.40 

3.40 

3.60 

1 .403 

+7 

2.170 

215.16 

-3312.76 

3.50 

3.70 

1 .631 

+3 

2.676 

166.03 

-3272.50 

3-70 

3.90 

2.1Q5 

+.1 

3-815 

113.30 

-3197.00 

3.90 

4.10 

3.694 

+.1 

5.332 

87.7° 

-3126.70 

-54- 


I'aole  A-3  Theoretical  data  for  the  input  admittance  of 
modified  dipole; 


(3  h 

0 

Poa 

Yo 

Z 

0 

0.60 

0.85 

0.130 

+,i  4.115 

IO.58 

-.1242.50 

0.80 

1.05 

0.856 

+,i  6.784 

18. 30 

-.1145.10 

0.90 

1 .15 

1.884 

+  1  3.012 

22.70 

-.1107.40 

1 .00 

1 .25 

4.573 

+  ;11 .950 

27.92 

-.1  72.98 

1  .10 

1.35 

12.870 

+  .114.430 

34.10 

-1  40.73 

1  .20 

1 .45 

22.820 

+  1  5.433 

41 .48 

-1  9.66 

1.30 

1.56 

17.360 

-1  6.465 

50.58 

+  1  18.83 

1  .40 

1.65 

9 . 802 

-1  8.037 

61.00 

+  1  50.00 

l  .50 

1.75 

6.240 

-1  6.738 

74.00 

+  1  79.90 

1 .60 

1.85 

4.447 

-8  5.450 

39.88 

+,1110.15 

l  .70 

1  ,n5 

3.437 

-.1  4.421 

109.60 

+1140.93 

1 .80 

2.0  5 

2.817 

-.1  3.6H 

134.20 

+  .1172.40 

l  .90 

2.15 

2.440 

-.1  2.953 

166.37 

+  .1200.64 

2.00 

2.25 

2.105 

-1  2.416 

205.00 

+1235-30 

2.10 

2.35 

1  .036 

-.i  1.956 

255.55 

+  .1258.20 

2.20 

2.45 

1.733 

-.1  1.548 

320.00 

+1236. 70 

2.30 

2.55 

1 .64° 

-1  1.195 

307.69 

+ j2do . 08 

2  .40 

2 .65 

1.515 

-.1  0.852 

501 .30 

+  .1281 .90 

2.50 

2.75 

1 .476 

-1  0.556 

5°3*16 

+  .1223.50 

2.60 

2.85 

1.384 

-1  0.243 

700.90 

+jl23.00 

2.70 

2  0  2 

1.341 

+8  0.047 

744.70 

-J258.50 

2.80 

3.05 

1.3H 

+1  0.334 

716.30 

-,1182.20 

2.90 

3.15 

1 .294 

+1  0.623 

627.20 

-1302 .20 

3.00 

3.25 

1 .200 

+1  0.922 

513.20 

404.40 

-.1366.70 

3.10 

3*35 

1 .301 

+i  1.235 

-.1333.70 

3.20 

3  M 

1.332 

+  1  1.570 

314.30 

-.1370.40 

3.30 

3.40 

3.55 

3.65 

1.387 

+  1  1.736 

244.60 

-1341.30 

3.50 

3.60 

3*75 

3.85 

1 .623 

+ 1  2.807 

154.40 

-.1267.00  i 

1 

3.70 

3.°5 

2.210 

+  1  3.971 

107.00 

-1192.30 

3.80 

4 .06 

2.732 

+  1  4.573 

96.28 

-1161.17 

q  on 
^  ^ 

4 . 1  5 

3  •  7?6 

+1  5.407 

35.06 

-.1123.00 

-55- 


,’able  A-4  theoretical  data  for  the  input  admittance  of 
modified  dipole; 


Poh 

poa 

*0 

z 

0 

0 .60 

C.91 

0.221 

+1 

12.06 

-.1233 .46 

0.70 

1.01 

0.452 

+  .i 

5.352 

15.68 

-il85.80 

0.80 

1 .11 

0.972 

+  i 

6.873 

20.14 

-1 142. 50 

O.00 

1  .21 

2  .  l4c> 

+  1 

8.970 

25.25 

-1105.40 

1.00 

1.31 

5.365 

+ jll. 910 

31.46 

-1  69.82 

1  .10 

1 .41 

12 .°40 

+  U3.100 

38.17 

-1  38.63 

1  .20 

1.51 

20.030 

+  1 

3.925 

46.15 

-.1  8.65 

1.30 

1.61 

15.730 

-.1 

5.855 

55.70 

+1  20.66 

1  .40 

1.71 

9.613 

-.1 

7.116 

6  7.20 

+.1  49.70 

1.50 

1.81 

6.331 

-1 

6.158 

81 .17 

+1  78.94 

1  .60 

!  1  .ni 

4.586 

-1 

5.065 

98.24 

+1108.48 

1.70 

2.01 

3.570 

-1 

4.142 

119.40 

+J138.50 

1.80 

2.11 

2.928 

-1 

3.333 

145.79 

+1168 .90 

1  .?0 

2.21 

2.00 

2.31 

2.18° 

-1 

2.258 

221 .30 

+  1228 . 30 

2.10 

2.41 

2.20 

2.51 

1 .791 

-.1 

1 .417 

343.40 

+1271 .65 

2.30 

2.61 

2.40 

2.71 

1.554 

-.1 

0.733 

526.48 

+1248.16 

2.50 

2 . 81 

1.523 

-.1 

0.438 

604.80 

+1173-58 
+1  64.01 

2.60 

2. *91 

1.407 

tJ 

C.123 

704.89 

2.70 

3.01 

1.357 

+  1 

0.161 

726.47 

-1  86.37 

2.80 

3.11 

1.321 

+1 

0.449 

678.50 

-.1230.67 

2.90 

3.21 

1 .2-8 

+  1 

0.741 

531.06 

-,1331.63 

3.00 

3.31 

1 .28" 

+  1 

1 . 042 

469.09 

-1379.29 

3.10 

3.41 

1.296 

+  i 

1.360 

367.32 

-1385.36 

3.20 

3.51 

1.323 

+  3 

1.700 

285.01 

-1366.34 

3.30 

3.61 

1.377 

+  i 

2.074 

222.17 

-1334.72 

3.40 

3.71 

1 .463 

+  .1 

2.491 

175.56 

-1297.95 

3.50 

3.81 

1 ,62Ci 

+  .1 

2.865 

£150  #  04 

-1263.79 

3.60 

3.91 

1.849 

+  1 

3.393 

123. 56 

-1227.10 

3.70 

4.01 

2.237 

+1 

4.161 

100.22 

-1136.45 

3.30 

4.11 

2.790 

+1 

4.736 

92.33 

-1156.75 

3.90 

4.21 

3.934 

+  i 

5.681 

82.38 

-1118.97 

Table  a- 5 


Theoretical  data  for  the  input  admittance  of 
modified  dicole; 


P0h 

v 

Yo 

zo 

0.6c 

1.01 

0.299 

+  .i 

4.507 

14.64 

-1220.90 

0.70 

1.11 

0.603 

+1 

5.602 

18.9° 

-1176.45 

0.60 

1 .21 

1.274 

+1 

7.122 

24.34 

-1136.50 

0.90 

1.31 

2.730 

+1 

^ .  073 

30.41 

-1101 .10 

1.00 

l  .41 

6 . 088 

+  .111.192 

37.50 

-1  68.95 

1.10 

1.51 

12.710 

+110.770 

45.81 

-1  38.80 

1 .20 

l .  6l 

17.430 

+  1 

3.H7 

55.60 

-1  9.94 

1 .30 

1.71 

1^.070 

-•1 

4.232 

65.14 

+1  19.61 

1  .40 

1.S1 

9.375 

~1 

5.675 

78.06 

+1  47.25 

1.50 

1 .91 

6 . 50S 

-j 

5.209 

93.65 

+1  74.97 

1.60 

2.01 

4.840 

-j 

4.421 

112.64 

+1102.90 

1.70 

2.11 

3.814 

-j 

3.671 

136.IO 

+1131 .00 

1  .so 

2.21 

3.143 

-j 

3.024 

165.24 

+1158.98 

1.90 

2.31 

2.745 

-j 

2.524 

196.26 

+  1181 .60 

2.00 

2.41 

2.343 

-j 

1.992 

247.74 

+1210.61 

2.10 

2.51 

2.217 

-j 

1.535 

304.96 

+1211.13 

2.20 

2.61 

2.019 

-.1 

1.173 

370.39 

+1215.18 

2.30 

2.71 

1.747 

-j 

0.850 

462.75 

+1225.26 

2  .40 

2.81 

1.737 

-j 

0.533 

526.13 

+1161.51 

2.50 

2.91 

1.528 

-j 

0.226 

640.36 

+1  94.48 

2.60 

3.01 

1 .450 

+  .1 

0.069 

688.10 

-1  32.50 

2.70 

3.11 

1.388 

+  1 

0.358 

675.60 

-1174.20 

2.60 

3.21 

1.340 

+  1 

0 . 648 

604.84 

-1292.29 

2.90 

3.31 

1.307 

+  1 

0.944 

502.97 

-1363.10 

3.00 

3.41 

l  .239 

+  1 

l  .251 

399.48 

-1337.73 

3.10 

3.51 

1.288 

+1 

1.577 

310.70 

-1380.28 

3.20 

3.61 

1.310 

+  .1 

l  .029 

240.05 

-1354.80 

3.30 

3.71 

1.361 

+  1 

2.316 

18 8.56 

-1320.92 

3.40 

3.81 

1.506 

+  1 

2.610 

165.84 

-1287.48 

3.50 

3  .°l 

1.613 

+  1 

3.248 

122.66 

-1246.98 

3.60 

4.01 

1.878 

+  1 

3.629 

112.46 

-1217.34 

3.70 

4.11 

2.263 

+  1 

4.263 

97.16 

-1183.02 

3. 80 

4.21 

2.905 

+1 

4.976 

87.50 

-1149.89 

-57- 


Table  a-6  Theoretical  data  for  the  input  admittance  of 
modified  dipole; 


Poh 

fV 

Yo 

z 

’0 

0.60 

1.11 

0.385 

*•-3 

4.698 

17.32 

-J211.43 

0.70 

1.21 

0.764 

+3 

5.784 

22.43 

-J169.94 

0.80 

1.31 

1.575 

+3 

7.235 

28.72 

-3131.97 

0.90 

1.41 

3.225 

+3 

8.926 

35.81 

-3  99.10 

1.00 

1.51 

6.581 

+110.306 

44.01 

-3  68.93 

1.10 

1.61 

11.845 

+3 

8.992 

53.56 

-3  40.66 

1.20 

1.71 

14.788 

+3 

3.139 

64.71 

-3  13.74 

1.30 

1.81 

12.550 

-3 

1.970 

77.79 

+3  12.22 

1.40 

1.91 

9.238 

-3 

3.709 

93.22 

+3  37.42 

1.50 

2.01 

7.060 

-3 

3. 960 

107.74 

+3  60.43 

1.60 

2.11 

5.131 

-3 

3.769 

126.58 

+3  92.97 

1.70 

2.21 

4.095 

-3 

3.192 

151.89 

+3118.42 

1.80 

2.31 

3.390 

-3 

2.649 

183.14 

+3143.12 

1.90 

2.41 

2.940 

-3 

2.182 

219.33 

+3162.78 

2.00 

2.51 

2.522 

-3 

1.723 

270.30 

+3184.70 

2.10 

2.61 

2.273 

-3 

1.320 

328.99 

+3191.06 

2.20 

2.71 

2.023 

-3 

0.971 

401.83 

+3192.85 

2.30 

2.81 

1.847 

_1 

V 

0.617 

487.00 

+3162.79 

2.40 

2.91 

1.707 

-3 

0.324 

565.37 

+3107.14 

2.50 

3.01 

1.601 

-3 

0.045 

624.02 

+3  17.46 

2.60 

3.H 

1.498 

+  3 

0.271 

646.41 

-3116.94 

2.70 

3.21 

1.421 

+3 

0.562 

608. 32 

-3240.63 

2.80 

3.31 

1.361 

+3 

0.855 

526.76 

-3331.00 

2.90 

3.41 

1.316 

+3 

1.157 

428.46 

-3376.74 

3.00 

3.51 

1.288 

+  3 

1.474 

336.22 

-3384.63 

3.10 

3.61 

1.280 

+3 

1.811 

260.28 

-3368.27 

3.20 

3.71 

1.296 

+3 

2.177 

201.93 

-3339.14 

3.30 

3.81 

1.346 

+3 

2.582 

158.80 

-3304.55 

3.40 

3.91 

1.532 

+3 

3.120 

126.82 

-3258.23 

3.50 

4.01 

1.618 

+3 

3.557 

105.98 

-3232.94 

3.60 

4.11 

1.920 

4.144 

92.05 

-3198.67 

3.70 

4.21 

2.390 

+j 

4.831 

82.28 

-3166.29 

3.80 

4.31 

3.310 

+j 

5.499 

80.35 

-3133.^9 

3.90 

4.41 

4.470 

6.180 

76.84 

-3106.24 

-58- 


Table  A_7  Theoretical  data  for  the  input  admittance  of 
modified  dipole; 


Poh 

poa 

Yo 

z 

0 

0.60 

l.3l 

0.534 

+  3 

4.946 

21.59 

-J199.85 

0.70 

1.41 

1.006 

+  3 

5.923 

27.88 

- jl64.10 

0.80 

1.51 

1.924 

+  3 

7.099 

35.56 

-,jl31.23 

0.90 

1.61 

3.498 

+  .1 

8.199 

44.02 

-jl03.18 

1.00 

1.71 

6.005 

+3 

8.710 

53*66 

-j  77.82 

1.10 

1.81 

9.044 

+  3 

7.630 

64.59 

-j  54.50 

1.20 

1.91 

10.986 

+  .1 

4.681 

77.04 

-j  32.82 

1.30 

2.01 

10.760 

+  3 

1.485 

91.22 

-j  12.59 

1.40 

2.11 

9.278 

-3 

0.542 

107.41 

+  j  6 .28 

1.50 

2.21 

7.669 

-3 

1.446 

125.91 

+j  23.73 

1.60 

2.31 

6.341 

-3 

1.706 

147.06 

+j  39.56 

1.70 

2 .41 

5.324 

-3 

1.655 

171.27 

+j  53.25 

1.80 

2.51 

4.553 

-3 

1.470 

198.89 

+j  64.22 

1.90 

2.61 

3.962 

-3 

1.230 

230.20 

+j  71.43 

2.00 

2.71 

3.500 

0.970 

265.30 

+j  73.50 

2.10 

2.81 

3.106 

0.726 

305.26 

+j  71.40 

2.20 

2.91 

2.782 

-3 

0.445 

350.47 

+j  56.11 

2.30 

3.01 

2.360 

-.1 

0.178 

493.38 

+  ,)  43.77 

2.40 

3.H 

1.940 

+  3 

0.089 

514.39 

-j  23.52 

2.50 

3.21 

1.742 

+  3 

0.395 

545.98 

-jl23.80 

2.60 

3-31 

1.591 

+  3 

0.700 

526.68 

-j231.77 

2.70 

3.41 

1.476 

+  3 

1.010 

461.37 

-J315.77 

2.80 

3.51 

1.382 

+  3 

1.330 

375.80 

—  j  36 1 .47 

2.90 

3.61 

1.309 

+  3 

1.663 

292.20 

-J371.30 

3.00 

3.71 

1.258 

+  i 

2.018 

222.49 

-J356.91 

3.10 

3.81 

1.233 

+  .1 

2.401 

169.30 

-3329.60 

3.20 

3.91 

1.242 

+  3 

2.821 

130.76 

-3296.94 

3.30 

4.01 

1.298 

+  3 

3.289 

103.92 

-3 263.23 

3.40 

4.11 

1.469 

+  3 

3.790 

88.91 

-j229.4o 

3.50 

4.21 

1.642 

+  3 

4.410 

74.13 

-3199.15 

3.60 

4.31 

2.129 

+  3 

5.110 

69.47 

-3166.75 

3.70 

4.41 

2.617 

+  3 

5.821 

64.25 

- jl42 .90 

3.80 

4.51 

3.690 

+  3 

6.710 

62.93 

- J114.43 

3.90 

4.61 

4.871 

+  3 

7.392 

62 . 16 

-3  94.32 

-59- 


,'able  A-8  Theoretical  data  for  the  input  adrittance  of 
modified  dirole; 


Poh 

poa 

Yo 

zo 

0 . 6 

i  .51 

0.520 

+  i 

5.053 

20.507 

-.1195.74 

0.7 

1.61 

0.925 

+3 

5.977 

26.840 

-3 162 .86 

0.8 

1 .71 

1.713 

+1 

6.922 

33.690 

-1136.12 

0.7 

1.31 

2.861 

+.i 

7.314 

41.310 

-1112.35 

1.0 

1 .91 

4.537 

+ ,i 

8.402 

49.760 

-1  92.lt 

1.1 

2.01 

6 . 64  0 

+  i 

3.258 

5°. 137 

-1  73.54 

1  .2 

2.11 

8.851 

+3 

7.046 

69 .490 

-1  56.60 

1-3 

2.21 

9.820 

+1 

4.938 

80. 950 

-.1  41.12 

1 .4 

2.31 

9.860 

+1 

2.839 

93.660 

-1  ?6.97 

l  -5 

2.41 

9.126 

+1 

1.197 

107.726 

-1  14.13 

1 .6 

2.51 

8.103 

+1 

0.173 

123.360 

- 1  2.64 

1.7 

2 . 6l 

7.087 

-1 

0.366 

140.720 

+1  7.27 

1.3 

2.71 

6.100 

-1 

0.597 

160.050 

+1  15.^2 

1  •? 

2.81 

5 .484 

-  1 

0.638 

179.920 

+1  20.94 

2.0 

2.91 

4.303 

-1 

0.564 

205.370 

+1  24.13 

2.1 

3.01 

4.314 

-1 

0.440 

229.449 

+1  23.39 

2.2 

3.11 

3.332 

-1 

0.250 

259.860 

+1  16.92 

2.3 

3.21 

3.484 

-  i 

0.068 

286.960 

+1  5.58 

2.4 

3.31 

3.128 

+  i 

fj 

Q.172 

318.740 

-1  17.56 

2.5 

3-41 

2.863 

+1 

0.386 

342.462 

-3  46.02 

2.6 

3.51 

2.596 

+1 

0.657 

362.030 

-1  91.70 

2.7 

3-61 

2.377 

+  1 

0.919 

365. °50 

-jl4l .49 

2.3 

3-71 

2.185 

+  1 

1.202 

351.300 

-1193.29 

2.o 

3.31 

2.019 

+1 

1.502 

318.760 

-1237.26 

3.0 

3.91 

1.880 

+  1 

1.830 

272 .°00 

-3265.94 

3-1 

4.01 

1.780 

+1 

2.148 

228.750 

-1275-84 

3.2 

4 . 1  1 

1.702 

+  1 

2.587 

177.460 

-1269.80 

3.3 

4.21 

1.6^3 

+.1 

2 . 966 

144.699 

-1255.10 

3.4 

4.31 

1.737 

*1 

3. 508 

113.400 

-1228.90 

3-5 

4.41 

1.939 

+  1 

4.005 

97.930 

-1202.28 

3.6 

4.51 

2.142 

+  i 

4.496 

66 . 356 

-1181.28 

3-7 

4.61 

2.731 

+1 

5.041 

83. 080 

-1153.36 

3.3 

4.71 

3-273 

+1 

5.576 

78.290 

-1133.37 

3.9 

4.81 

4.050 

+1 

6.271 

72.670 

-1122.50 

-60- 


Vable  a-9  Theoretical  data  for  the  inout  admittance  of 
modified  dipole; 


Yo 

z 

0 

0.6 

l  .61 

0.477 

+,i 

5.138 

17.93 

- jl°2 .98 

0.7 

1.71 

0.853 

+0 

6.006 

23.19 

-jl63.22 

0.8 

l.Sl 

1.457 

+.1 

6.033 

29.04 

-J138.14 

O.o 

1  .41 

2.307 

+  0 

7.657 

35.53 

-j ll 6. 44 

1.0 

2 .01 

3-785 

+  .1 

8.626 

42.65 

-j  97.22 

1  .1 

2.11 

5.650 

+j 

8 .  o4l 

50.51 

-j  79.93 

1.2 

2.21 

7.764 

+  .1 

8.429 

59.12 

- j  64.18 

1.3 

2.31 

9  •  55° 

+  0 

6.928 

68 . 59 

-j  49.71 

1  .4 

2.41 

10.442 

+  J 

4.804 

79.04 

-j  36.36 

1.5 

2 .51 

10.315 

+  .i 

2.733 

90.57 

-j  24.04 

l  .6 

2.61 

9.530 

+  .1 

1.172 

103.37 

-j  12.71 

1 .7 

2.71 

8.500 

+  0 

0.176 

117.65 

-j  24.41 

1.8 

2.81 

7.470 

-0 

0.370 

133.60 

+ j  6.6 l 

1 .9 

2.91 

6.540 

-  0 

0.612 

151.55 

+J  14.19 

2.0 

3.01 

5 . 744 

-0 

0.665 

171.79 

+j  19.90 

2.1 

3-11 

5-067 

- 0 

0.596 

l 94 . 64 

+J  22.91 

2.2 

3.21 

4 . 404 

-0 

0.454 

220.27 

+J  22.25 

2.3 

3.31 

4.004 

-0 

0.264 

248.70 

+j  16. 4l 

2.4 

3.41 

3.501 

-0 

0.039 

27°. 18 

+3  3.01 

2.5 

3.51 

3.215 

+  .1 

0.213 

309.65 

-j  20.47 

2.6 

3.61 

2 . 8°7 

+  0 

C  .480 

335-95 

-j  55.68 

2.7 

3.71 

2.618 

+,1 

0.773 

351-38 

-OIO3.76 

2.8 

3.81 

2.375 

+0 

1.083 

348.62 

-0158.93 

2.9 

3.91 

2.166 

+5 

1 .41 6 

323.42 

- j211 .44 

3-0 

4.01 

I.092 

+j 

1.777 

27°. 60 

-0249.31 

3.1 

4.11 

1.857 

+.1 

2.16° 

227.84 

-.1265.99 

3.2 

4.21 

1.769 

+3 

2.592 

17°. 66 

-0263.25 

3-3 

4.31 

1.738 

+0 

3.042 

141.57 

-0247.83 

3.4 

4 .41 

1.784 

+.i 

3-528 

114.12 

-0225.73 

7.3 

4.51 

l  .n0l 

+  0 

4.093 

93.34 

-j200.°7 

3.6 

4.6i 

2 .104 

+  .1 

4.657 

82.79 

-0175.73 

3.7 

4.71 

2.624 

+.1 

5.243 

76.21 

-jl52.44 

3-3 

4 . 8l 

3.240 

+3 

5.723 

74.91 

-0132.32 

3.9 

4.91 

4 . 04  0 

+.1 

6.035 

76.60 

- jll4 .41 

-61- 


Table  A-10  Theoretical  data  for  the  input  admittance  of 
modified,  dipole; 


V 

IV 

Yo 

zo 

0.6 

1.71 

0.411 

+0 

5.248 

14.84 

-j 189. 38 

0.7 

1 .01 

0.729 

+  .1 

6.122 

19.18 

-j l6l .07 

0.3 

1 .91 

1 .237 

+  i 

t 

7.071 

24.00 

-J137.22 

0.9 

2.01 

2.031 

+0 

8.063 

29.34 

-J120.20 

1.0 

2.11 

3.234 

+  0 

9 . 024 

35.19 

-J104.32 

l  .1 

2.21 

4.958 

+  0 

9.721 

41 .64 

-j  81.64 

1.2 

2.31 

7.180 

+0 

9.791 

48.71 

-j  66.42 

1.3 

2 .41 

9.534 

+  0 

8.831 

56.51 

-j  52.29 

1 .4 

2.51 

II.298 

+  0 

6.772 

65.12 

-j  39.03 

1.5 

2.61 

11.889 

+  0 

4.223 

74.69 

-j  26.53 

1.6 

2.71 

11.374 

+0 

1.953 

85.40 

-j  14.67 

1.7 

2.81 

10.249 

+j 

0.364 

97.45 

-j  3.46 

1 . 8 

2.91 

8.960 

-0 

0.571 

111  .16 

+j  7.09 

1.9 

3.01 

7.746 

-0 

1.024 

126.88 

+j  16.77 

2.0 

3.U 

6.688 

-0 

1.166 

145.10 

+j  25.29 

2.1 

3.21 

5.795 

-0 

1.120 

166.36 

+j  32.15 

2.2 

3.31 

5.046 

-0 

0.956 

191 .32 

+j  36.25 

2.3 

3.41 

4.416 

-0 

0.724 

220.49 

+j  36.16 

2.4 

3.51 

3.886 

-.1 

0.446 

253.97 

+j  29.12 

2.5 

3.61 

3.435 

-0 

O.I38 

290.65 

+j  11.67 

2.6 

3.71 

3.049 

+  0 

0.194 

326.62 

-j  20.81 

2.7 

3.81 

2.720 

+  0 

0.545 

353.40 

-j  70.81 

2.8 

3.91 

2.441 

+  0 

0.917 

359.05 

-jl34.82 

2.0 

4.01 

2.207 

+0 

1.309 

335.26 

-jl98.80 

3.0 

4.11 

2.01O 

+  0 

1.722 

286.79 

-j 244. 50 

3.1 

4.21 

1.880 

+0 

2.157 

229.69 

—  J263 . 51 

3.2 

4.31 

1.795 

+.1 

2.613 

178.60 

-J260.00 

3-3 

4.41 

1.780 

+  0 

3.091 

139.91 

- j242 .95 

3-4 

4.51 

1.829 

+0 

3.625 

110.95 

-0219.91 

3-5 

4.61 

1.973 

+0 

4.263 

89-59 

-J193.02 

3.6 

4.71 

2.245 

+0 

4.349 

78.62 

-jl69. 84 

3-7 

4.81 

2.651 

+0 

5-399 

73.28 

-J149.25 

3.8 

4.91 

3.216 

+  i 

5.889 

71.43 

-J130.81 

3.9 

5.01 

3.944 

+j 

6.253 

72.08 

-0114.38 

-62- 


Table  a-11  .heoretical  data  for  the  input  admittance  of 
modified  dipole; 


Poh 

(30a 

Yo 

zo 

0.6 

1.91 

0.2-1 

+j 

5.534 

0.49 

-jl8C .20 

0.7 

2.01 

0.518 

+  .i 

6 .46° 

12.31 

-3154.08 

0.8 

2.11 

0.889 

+  3 

7.525 

15.48 

-313I.O6 

0.9 

2.21 

1 .4°0 

+.1 

8.730 

19.01 

-3IH.31 

1.0 

2.31 

2.471 

+  3 10 . 082 

22.94 

-3  93.56 

1.1 

2.41 

4.06° 

+3: 

11.507 

27.31 

-3  77.24 

1 .2 

2.51 

6.607 

+312.712 

32.19 

-3  61.93 

1.3 

2.61 

10 .2°2 

+J12.932 

37.63 

-3  47.34 

1 .4 

2.71 

14.4^6 

+J10.967 

43.88 

-3  33.19 

1.5 

2.81 

17.159 

+  .1 

6.400 

50.97 

-3  19.30 

1.6 

2.91 

16.761 

+  3 

1.554 

59.15 

-3  5.49 

1.7 

3.01 

14.384 

-.1 

1.753 

68 . 74 

+3  8.4l 

l  .8 

3.H 

11.569 

-3 

3.246 

80. 13 

+3  22.48 

l  .9 

3.21 

9.233 

_  i 

3.611 

°3.94 

+3  36.74 

2.0 

3.31 

7.434 

-.1 

3.427 

110.95 

+  3  51.14 

2.1 

3.41 

6.073 

-,i 

3.000 

132.29 

+3  65.31 

2.2 

3.51 

5.050 

—  i 

2.473 

159.59 

+3  78.32 

2.3 

3.61 

4.260 

-3 

1.929 

1 94 . 82 

+3  88.21 

2.4 

3.71 

3.642 

~,i 

1.382 

240.01 

+3  91.06 

2.5 

3. 81 

3.156 

0.349 

295.51 

+3  79.53 

2.6 

3.91 

2.766 

-3 

n.329 

356.55 

+3  42.41 

2.7 

4.01 

2.455 

+3 

0.209 

404.41 

-3  34.46 

2.8 

4.11 

2.206 

+  i 

0.661 

415.89 

-3124.69 

2.9 

4.21 

2.049 

+  i 

O' 

1.160 

369.60 

-3209.21 

3.0 

4.31 

1.910 

+.1 

1.622 

304.00 

-3258.19 

3-1 

4.41 

1.823 

+,i 

2.242 

218.26 

-3268.55 

3.2 

4.51 

1 .7C2 

+3 

2.824 

160.21 

-3252.48 

3-3 

4.61 

1.820 

+3 

3-391 

122.85 

-3228.96 

8.4 

4.71 

1 .914 

+3 

3.957 

99.09 

-3204.80 

3.5 

4.81 

2.086 

+3 

4.533 

87.79 

-3182.04 

3.6 

4.91 

2.351 

+j 

5.H7 

74.14 

-3161.36 

3.7 

5.01 

2.727 

+  j 

5.697 

6c8.35 

-3142.80 

T  P 

5.11 

3.236 

+.1 

6.254 

65.26 

-3126.12 

3.9 

5.21 

3.901 

+j 

6.754 

64.13 

-3111.02 

-63- 


Table  a- 12  theoretical  data  for  the  input  admittance  of 

modified  dipole; 


poh 

poa 

Yo 

*0 

0.6 

2.11 

0.217 

+  3  5.854 

6.319 

-3170.58 

0.7 

2.21 

0.393 

+  3  6.881 

8.268 

-3 144. 86 

0.8 

2.31 

0.691 

+  j  8 . 084 

10.492 

-3122. 80 

0.9 

2.41 

1 .202 

+3  9.530 

13.032 

-3103.29 

1.0 

2.51 

2.107 

+jll .308 

15.926 

-3  85.46 

1.1 

2.61 

3.771 

r 313 .488 

19.230 

-3  68.78 

1.2 

2.71 

6.939 

+315.911 

23.030 

-3  52.81 

1.3 

2.81 

12.828 

+317.406 

27.440 

-3  37.23 

1 .4 

2.91 

21 .226 

+314.161 

32.600 

-3  21.75 

1.5 

3.01 

25.196 

+3  3.062 

38.730 

-3  6.09 

1.6 

3.11 

20.705 

-3  4.484 

46.134 

+3  9.99 

1.7 

3.21 

14.673 

-3  7.108 

55.200 

x3  26.74 

1.8 

3.31 

10.401 

-3  6.938 

66.540 

+3  44.38 

1.9 

3.41 

7.682 

-3  5- 981 

81.047 

+3  63.10 

2.0 

3.51 

5.926 

-3  4.961 

99.960 

+3  82 . 92 

2.1 

3.61 

4.745 

-3  3.922 

125.220 

+3103.50 

2.2 

3.71 

3.062 

-3  3.029 

159.768 

+3123.53 

2.3 

3. 81 

3.315 

-3  2.239 

207.175 

+3139.93 

2.4 

3.91 

2.864 

-3  1.497 

274.258 

+3143.32 

2.5 

4.01 

2.515 

-3  0.948 

348.190 

+3131.19 

2.6 

4.11 

2.347 

-3  0.35 

416.81 

+3  62.16 

2.7 

4.21 

2.105 

+3  0.248 

465.453 

-3  55.25 

2.8 

4.31 

1.958 

+3  0.870 

426.473 

-3189.59 

2.9 

4.41 

1.858 

+3  1.461 

332.539 

-3261.61 

3.0 

4.51 

1.799 

+3  2.015 

246.564 

-3276.09 

3.1 

4.61 

1.782 

+3  2.554 

183.727 

-3263.33 

3.2 

4.71 

1 .807 

+3  3.090 

l4l .020 

-3241.18 

3-3 

4.81 

1.877 

+3  3.634 

112.202 

-3217.20 

3.4 

4.91 

2.001 

+3  4.191 

92.781 

-3194.31 

3.5 

5.01 

2.189 

+3  4.767 

79.555 

-3173.25 

3.6 

5.11 

2.456 

+3  5.361 

70.621 

-3154.17 

3.7 

5.21 

2.823 

+3  5.973 

64.688 

-3136.85 

3.8 

5.31 

3.320 

+3  6.594 

60.014 

-3120.98 

3.9 

5.41 

3.984 

+3  7.205 

58.772 

-3106. 2Q 

-64- 


APPENDIX  B.  COMPUTER  PR#OGRAMS  FOR  MODIFIED  DIPOLES  (1st  set) 


FQRTr.JN  IV 


L  :v  *  l  1  L} 


MAIN 


00  31 
0002 
0003 


F  3R7R  ' 


IV  PRCG~ 


CATE  =  71110 
•CO  1  f  I  EC  CONICAL  Af 


15/ lb/39 


C 

c 


r 

Tr 


c 

c 

c 

c 

c 

c 


AL 

Cl) 

TH 

i  P 

IS 

L* 


HIS  PROG  I<  AM  J0LV5S  • 
IhE  fvCO'Al  ICN  PFGICU  or 
1  CHAR  ACT r*T5T IC  FUUA.  iCf 
ic  ANGLE  ANC  »  H r- ! .  >•  :  . 

■jEP.RAJC  r:\’u  v  7  J  'U  5  .  "r«r  « 

krent  oisi  irun:r<s 
E  IriPUT  C  H  .*  p  A  C  i  £  R  I :»  7  I  C  :»  , 
P*>.  ATT  K.'.N'  L  ‘  O'rl'  v  >  C.  nm 
V  A  L  i  D  P  C i  i  A  L  r  A  N  G 1.  i  L 
TO  5  AND  Cco-kAL  SPr^f 


-he  COEFFICIENTS  D,-:*  OF  EL  ECT  L  CKAGNET  F  C  FIELD 
;i-c.  modifies  con  ic a l  •'MEnv*  o.y  cjp.St  sclv1.:;.: 

L  6.  )  -•  •**  vs  Ml  Crt  GIVES  A  SET  C*  \  S  FO*  A  ritFC 
V.  V-FLr  T-UVC.UfNG  7t»r  !  W  T  1  r«  »  <  2  T  Of  Lt-irAi 

i  ec  » p -3 v* a .i‘, l r i :  fields  r<  T  knmuna  afci:;n, 

V  5 '  •  SURFACES  Cr  TmE  P  AT.  I  A  'f  I.'JC  ST'uCT’JSC  *.  M; 

P Z  EXPRESSES  IN  TEAMS  Or  t/S*  TCGlTVEP  KITH 
• :  FUNCTION  A  NO  VESSEL  $  **  U  N  U  7  1  C  W  S  •  !  H  I  S  >f>C-Gi 

35  1>  VU  IwJ  CEC-U'f:  S  AN;)  ANTENNA  iLICH*  c!  FT*  * 
el1  a-'  op  rc 


EXTt  ftNAl  fCT 
r*N (  }<¥.  *i  >=KKA(KK  <1  1  —  Z  v  <  i 
COMPLEX  v  (  50  •  FO  I  ,.At50#S; 


1) 


1 20)  f  CUflCI  .301  fU-TiSt’O  *  i CO 


>  H  {  If 

:i  3-;  i 


0  1  ME  N  S I  :  tN  FM7SI  *f  S  i:0,5  3) ,  F  I  (  50  <  SO  I  *p 

1  , f :iH  T  (  /*>  )  ,  ASr.  (  3(0  ,  GriOHfM  .u  ))  ,OKXJP  (6  0  )  «G 

2  PNM  <110/,  Cr'N  0(1 2  I  t  P A L  <  )),  Cf  NL  <60  )  , 

30)  •  B  L  T  A  <  60?  ,  GY1  $U  f  6  0  ,P?>  ,0Y  ‘  <  77  (  60  ,  ?9 
4,r>K  <  bO)  ,F<  VC  )  ,C US  u:  J  ,Cl?9  <60)  .C7JO  Jp<  60  > 


1  JO  )  . 

DTBl c0 ) 

i  0  S  0  S 

(  2  »:■■’  i  • 

r  CO 

MC 

( 

3)  , 

,  a  \ 

MG» ADMl 

» At; M2 

*  i\  C  M  3  1 

,  A  C-M  A 

? 

\  L/:' 

iS? 

» CUN S3 , 

CUN$<* 

t  ( ; 

K: 

/  f op  1  ( : 

10  '  ,  I 

t::t  (- 

;  0  • 

i>C 

) 

?M> 

■  i  i 

51), PNP 

<  1  i  0  ) 

,CPN?< 11 

r*  i 

» 

:t  1 

L  ( 

60  i  ,OYx 

.jp' 

'  nSN 

1>.J 

1  < 

r 

1 1^. 

>  M  ? 

OMIbO.J 

0  5  ,  OM 

1GKP  (J 

:  0  ♦ 

20  ) 

c 

c 

c 


IE  AO  l KrCT  OATA 

H  A  v-  M A  y  T  ,VI  \Mt  (J f  OF  K 

/ax:,;o.*  cr 

TA-r.M:-  CCNC  /Nil.  E  In  CEGPcLS 
TA-OH  i  VINO  {»ftF  i  N»  DEGREES 


00  )5 

1  <'  FAC  <  >,?)  WAX,  NF*AX,TAf  BETA 

CCC6 

2  F C A f •  A T <  15  f  15  »2  JP!  )  ) 

0  30  7 

CALL  ERRSFTS20  7.  :rfltliO*JMPr.  R  ? 

C  •?  C  5 

CALL  ERRStTf 2! 0. 2,1,1, PUMPER) 

CCCs 

MAXX  =  r MAX  +  1 ) /2 

0  010 

Ml -HAM*  1 

0011 

M2  - M  +1 

001  2 

MAI  - MAX* 1 

0013 

0=1 rf- 07 

CC1  A 

PIC- 3.1  a; 59  2*5 

0015 

TAR^FIOTA/lcO 

-  C 1  6 

T  AF  B  =  P I C  * ( T  A  *  S  E  T  A )/lQ0 

0 Cl  7 

TAB  C--TARB/2 

CO  I  P 

T  A  2  =  T  A  P  /  2 

CO  19 

IT^=lf I X  < TA)*1 

C  0  2  0 

CF r 1 TA-1 A 
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00  2  5 

KIT  £9  =  ALCC-«CCT;iN;TA2  )/  Cl  ?  AN  <  T 

002/: 

m  T.\  =  need*  irji  a:  >  > 

0  y2  7 

t  \  “I  G  3 •  ' F> i  1  A 

(  CO 2 

w  >  - 1:  C  S  (  T  A  N  ) 

' )  i 
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OPTfiAN 

17  .  Lt/n 

1  ►:  Mi  I’i  r,4TE  -  71110 

002  ^ 

r 

c 

i  x  =  i  I  ,J  ( T  '  ; 

c 

SCl\(  ii-iCTc'*!''!  .  CCotTIOfi  l  <  9.  >  =  0 

c 

R7=  :  Or?.  <:r  i  ,<  i-m  *  <> 

r 

CM-  : 

c 

0030 

CALL  PTM  ( ,1  M  ,M  w  ,5  r'<r  f  fCT, -?.2  ,  !  7  *  If  -  ,  I  ,  P  T  .  -\;T  ) 

0  03  i 

1='-  -?  i  /> 

00  5" 

i.»C  301  J<  ►'  SV'.', 

:  r  ?  ? 

C T‘ b ■  c  i  l  » 

003- 

CALI  IF  rr  f  T  T  '  O  ,r.  rj>  ,'irn,  •.  .  ^\p  ,  .^M*-  ,s»N  V  ,0*-  N*  ,  P  M  ,  0  PM.  » 

003  3 

<  L  T  A  {  h  ' 4  w-  f:  t  <  >i ) 

C  C  3  6 

301 

r 

s» 

com  i.v;: 

r 

V 

r  c  A  U  3  N  P  0  I  C  .A  T  A 

o 

h7  -  4NU,V.A  hi  I  AM  7  *  :**•  VI  r CT ? t CAL  LENGTH 

c 

rS-’  CEMhAC  SPUES1  k£,  fjs  IN  ElKTRICAl  L  SHUT  D’ 

r 

KON-  NWOcB  CL  rU’C.'L  f.  FOMENTS  /  «.OWG  7  r  »  &tfT£.*'NA 

u 

c 

c 

FOiJ  whlCr  CCrMONl  L'iIOTRI  3LTiON7  C CPPuTfc 0 « 

0C3  7 

3 

R  EAG(  6*4  ) H7  , 

0C3  6 

HIM  MRS 

OC34 

■rtRilElfc.4JHT»«S.K«N 

0C40 

4 

FORMAT  1  ;n?.  4  •  15  > 

C  "41 

CT*-2«M  AM  HI  -RSl/ftLTA 

f.|)4  > 

f  A  l  L  '  1  •  •  .i  T  «  SH  ?  L/S r»  »  .i  .#*>,>  « 

004  3 

00  3  0  0  K  L  -  1  »  A f'  A 

0044 

0  L  C  *  R  T  <  K  L  ) 

0045 

CALL  LrGf  (  IT  A,  OP  ,  C  L  0  .•  ~  f  PM  .  FiPfcP,  t  DPM.^  *PNL  ,  0PM  > 

0046 

CALI  GT  Jf.fM  LI  C  ,HV  ,  Vi,  ,  KI-SN  ,  CMC.-VP  ,  o*XJP,  GW-iX  J ,  C*X  J?\l  ,G?PXJ  1,07 
2ST1  ,GY].$T2  ,Gf'l  CJ*  tr,PICPPfGJl  LJ,!Dc*  ,  I  » 

CCA  7 

3  CO 
C 

CCt»J I NUC 

c 

CENERA1F.S  FATRIX  £Lc«tMS 

CC4  6 

c 

00  1  C  3  ,MX,2 

CCA  <5 

IS-<  KM  1  )  /.? 

0050 

ASM 

0G5 1 

DO  102  KK  •  r'r  ,  ^  A X  ,  2 

00  52 

KK2-KKf4 

005  3 

|K-?KK*1 I/O 

0  05  4 

I  Tt  sm*,  IK^O 

0055 

r S\  ip :f  i  y )  r0* 

0066 

i  f.  x  -  0 

0  7  5  7 

f‘0  1J  I  N  -  1  f  N  ^  A 

CCS  b 

1  T  ~  S  7 1  IP,  I  K  )  -  I  T  f  S  T  ( I  h  » I  K )  ■•  1 

005^3 

C^C - R T% I  A ) 

C  06  C 

*M  U«  J  2*-0?  f  ♦  1  •  )/t'-P£f  (  A'  C  <  1  .  )  ) 

0061 

PM  I  N  )  -  P  (  JM*H7  A(  !  K  »  *INT  (  l\  )  /  (  SX  ^5  X  M'M  <  IN  5  ) 

006  2 

F  5P  -  P  (  f  N  ?  *  Mi  i\  }/  (  g:  /;•  P  (  IN  »  .<0,  (  KK  ♦  r  PO)  ) 

"  ^  .6  " 

,lK)=r  M  I  «*' ,  T  K  )  ♦  P  .SP 

flMf  =  AOcj(F  SF/PSCIP*  IK)  I 

Tr6  5‘ 

1F<  HR  .Mr  .KK>GC  T C  17 

I  CX-  i 

CC67 

ASPh-?»;<  :m/:.  hM^.cr^v.Vivjpj  (  in;  i 

1  «/U/39 


The*  S  P  W£K 


C  J  <  1  V7  Y  1 
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FORTRAN 

C06  6 
oot>5 
CC7C 
0071 
0072 
0  C  7  3 
CC  74 
0075 
0076 
0077 
007? 
0C79 
00P0 
CC61 
0032 

C  C  c3 
CC  S 4 
0  C  6  5 
C  C  86 

cce7 
0038 
C  Ce  9 
CO^O 
0091 
0092 

tC.  'O 

0C94 
CCS  5 
CCS6 
0097 
0098 
CC99 
0100 
0101 
0102 
0103 
0  104 
0105 
01  06 
0107 
01  C  8 
0109 
0110 
011  1 
0112 
0113 
C 1 1  4 
0115 
0116 
0117 

out 
Cll  9 
0120 
0121 
012? 


W  G  LEVEL  13 


MAIN 


CATE  =  71 1 1 C  19/16/3* 


ASi\<  m=A$B<  1  R  J  *  A  S  P  B 
R ATA  =  AB$( ASPB/ASk ( I R ) ) 

17  I F ( CR  C—  K  K  2  1 1 0  » I  0  » 2 L 
21  I F  (  IGX.LT.l »G0  TC  25 

IPiRATF.LT  .1  .t-7'f  ANC.BMA.L  T.1.E-04>C0  TC  19 
GG  70  1C 

25  IF(FATF.LT.l.r.  .J4IGC  TO  19 
10  CONTINUE 

19  F  I  (  »R  »  m«-2  1  >«PI  (KK«  t  l*F3  <  IR»  IK  I 

F  I  (  IK  r!  ft  I  *F  I  (  iw  ,  I  Kl 
Iff  h  R  . :  C  *  P  K  )  G  C.  I'”  51 

V(  IR  ,  I  Kl  =C  T*T  I  (KB  +  l  )  *PI  (  Kr;  +  )  )  /  (  K  m*KK  *  ( KF  «  1 1  *  (  K  K  ♦  1  )  ) 

V(  IK  ,  I R  >  -  V  (  I»,  Ifv) 

GO  TO  101 

51  V ( l P , I« ) =CT*PI <KK*1 >**2/ ( KPXKR  +  l ) I **2-2. M l ./FT ♦ G 3 0 S ( K R ♦ 1 1  I / ( ( 2 • * 
1KR-*  1  )^KR’*  (  KB  ♦  l )  ) 

101  A [ i F  , i K)=f I ( I ?, IK ) +V< IB , I K 1 
A  (  I  K  «  I  R  !  -  A  (  I  R  »  I  K  ) 

102  CONTINUE 

DT  181  IR  )*P  I  <  KR*1  l*<  -?/  I  2C*CO$  (HT-RS  I  *KRM  KR  *1 1  I  OSo!  IR  I  /  <  &  0*RLTAB  I 
1  ) 

DTR  {  IP  )=CKHX0.  ,-OTI E< IR>  > 

A ( IR  fFl  ) =DTP (  Ift  ) 

103  CCN1 INUE 

CALL  CS  l  FTC ( A , 60 , FA\X , 1  ) 

HPnE(6,7JhT,R$ 

7  FORK  AT( /////46H  **'*  ANSWER  PC R  COEFFICIENTS  Cf:  FIELD  ECR  HT=F10«4 

lflX'/MAND  R$=  F10  •  **/5  H  I#10X,5H  REAL,?0X,5H  I  MG  I 

do  r  :  :  r  ?  a  x  >: 

WR|Tb(6,9)lfAU,Ml> 

<5  FCRK  AT  (  15,2(1  OX  «  F.  15*3«) 

8  CONI  I NUC 
RIVC=(0.,  1  •  ) 

A0K1*BIKG*TAN(HT-PS> /ZC 

acm; mc.  ,o. ) 

CO  7 CC  Kl =1  , M  A  X  X 

ADM  =  A U  F 2  +  P I  <2*f  1  I  *  A  (  Kl  ,  Ml  )  /  (  2*K1M  2*K1-1  )  ) 

700  CCNT INUE 

ADM  =  120*ACM2/(  2C*CGS(HT-RS  )  ) 

ADM}- (0. ,C,  ) 

DO  9  01  K.2  =  1  ,  N  VA 

AD3-  PK(  K2  MGY1ST1  (K?  ,  1 ) 

ACMr=ADR2*A03 
9  AT  - CA3S ( AC  V AOM3 ) 

I F\ RA T. I T . 1 .E-C3 )GG  TC  901 

701  CONTINUE 

901  AC?'3  -F  IMGOCM/ (120.  +  RLTALM 
A0M<  =  (0,  ,C.  ) 

DC  703  K3-1.KAXX 
K  :  v  K  3 
A  0  4 - ( C. ,C* ) 

00  7  0**  K 4 "  1 , N M A 
CKC-RT (K4 } 

AC  3  r  P  {  K  4  )  *  G  M  1  C  f  ‘  (  K  n ,  l  >*D.NT  (K4  >  /  fK.NI  (  K  32-1  >  ,OXC  ) 

A  L*.  -  A L ♦  A 0 3 

IT  i  (  f  Kf,  •  I  •  L  T  •  K 2  I  GO  TO  7  04 
ft  AT  *  *•*'*><  A03/ADA  ) 

IF  f  f*  a  f  •  L  T  •  1  •  E  *  0  3  I  GO  TC  902 
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FORTRAN 

IV  G  LEVEL 

16  MAIN  DATE  =  711 1 C 

012? 

704 

COM  IMl 

0  1 

902 

*CK4=An*A-MK3tPl  )  "PI  (2*K3)  *AC4 

0 1  2  b 

RATI  -CA  1  (  *  C4';A  (  K3,M  M  PI  (’*K3)  /AC^4 

Cl  2  6 

IF  l  RAT1.  *17*1  .  p-02  ):a  TO  903 

012  ? 

703 

C  0  N  T  I  ;  i  I  j  E 

012  8 

903 

ADM*  2*  (  A  ry :  ♦  AC**;:  ♦  A  '3  *■  A  CM  4  > 

012  9 

r  im=i  ./  a: 

C  1  30 

AAC.V*C  AGS  I  AlM  ) 

C 12  1 

AR  IM  =  1/  A/.CV 

0132 

TA0M-A1  ;•:>?(  AlMAr,  (/.»::/ If  prai.  (ACM)  »  *160 /  PIT 

012  3 

TRIM  -  ta;.> 

0134 

r.0  IT  9(6,16) 

01  3  5 

16 

FORMAT  (/////  4 1  H  *  INPUT  A  CM  I  TT ANC9  AND  IMHECANCE 

1TANC  E  ,50X  ,  UM  IMPtDANCE  ) 

0  1  3( 

WRI  T  c ( 6 , 1 ? ) 

0137 

W*  1  T&(  6,  OC-  )  AD?- ,  A.VJP,  TACM  »P  IM-tAR  IP, TRIM 

0136 

hR  IT£  (  6, 27)  ACM  1  ,AO?  ,ADM3  ,ASP4 

013  9 

27 

FORMAT(jX#E15«Cf?/*»F15»?/> 

0140 

1  £ 

FORMAT (2(2* »6H  PtAL  ,10X,6H  1  MAG  »20X»11H  MAGNITUDE 

0141 

904 

FGRMAT(2(?X,cl5*e)f2>tfrir*6»2X,F6.:»?(2X,E15#e),0XfF 

0142 

SR  IT*:  (6,2  0) 

C 

C  jM  PLIES  CLFRENT  OISTR  IBLT  I  C  N  S  ALONG  ANTFKf.A 

01  h  3 

20 

FORMAT ( ////4  9  H  **  CURRENT  DISTRIBUTION  ALONG  CONICAL 

0144 

D I V  = (HT~R$)/KMN 

0145 

KMN1-KMN4 1 

Cl  46 

CUNC (  1  ) =0  *  5* ACM 

0U7 

C  UNO ( 1 )  CUNC ( 1 ) /(2”P!C*RS*SX) 

oi^e 

ACIJNC  =C  AES  (  CUNC  (  1  )  ) 

o\ 4? 

AC.  N  P  ►  r*  l  >.  O  u*  i  k  j  < 

0150 

T  CUNC* AT AN2 ( AIM AG (CUNC (1 >  ) , PEAL (CUNC ( 1 M ) *180/ PIC 

0151 

TCUN3=TCUNC 

Cl  52 

SR  I  TE  (  6 , 9  04  )CUNC  (  1  )  ,  ACUhC  ,  T  CUNC ,  CUNC  ID,  A  CUN  C ,  TC'JNC 

01  53 

00  405  K5'2,KMM1 

C 1  5  4 

HTC=R$«C!VM K5-1 ) 

01  55 

CUNC 1 *-RIPG* ( $IN(HTC-RS)-CGS(hTC-RS ) M AN (HI -RS >  >/ 2C 

0156 

CUNC  2  =  (  c. ) 

Cl  57 

DC)  3 OC  Kc  - 1  ,  MAXX 

0158 

CUNC  2  =  CUf«'C2  ♦Pi  (2*Kc)*MK6,Ml)/(2AK6*  i2*K6-l  )  ) 

C  1  59 

800 

continue 

OlfcC 

CUNC2-CUNC?*CGS(HTC-RS)/ (CCSU  T~RS)*RLTA) 

016! 

C  UNC  3  -  (  C  ,  ,  0  •  > 

C!  12 

I  F  (  K5.rC.KMM  )GC  TC  666 

0153 

DO  6fl  K7*1,NMA 

01  64 

AD3*PK(  K7  ^GYISU  (K7,K5) 

Olt  5 

CUNC  ?*CUNC 2  +  AC3 

0166 

*  AT«C  AB  $ ( AC3 /CUNC? ) 

01  67 

IF(RAT.tT.i#F-03 )GC  TC  1^01 

C 1  6  6 

901 

CCNTiN-jr 

016  9 

1  101 

CUNC3 9  I  MU*  CUN'C  2/  (  KCNLTA3  ) 

01  7  C 

66c 

CLNC4 “ ( 0  «  ,0. ) 

0171 

DO  3C3  KC-1 , 6AXX 

0172 

KC2=2-'M 

01  73 

CL.\4r(D.  ,c.) 

0174 

(;U  PCs  K  4=1,  NR  A 

Cl  75 

CKC  =  M(  K9) 

01  76 

A  H 3  -P  (  K -V  1  -G.V  1  C'M  M.Kl  )  ^CNT  (K9  )/PK?i(  (K  82-1  )•  CKO 

C  1  7  7 

C  UN*.  -  t  UN 4  ♦  AO 3 

01  7  t 

lr((CKC-V«).LT«KsZ)GD  TO  8C4 

****/ 1 1 


1 6X i 6M 

1  5  •  8  •  2  X 


»  SUk p sC 


18/13/?? 


H  A CM  I  T 

PHASE  )  ) 
t  F  6  •  1 ) 

r  #*  ) 
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FORTRAN 

IV  0  level 

16  MAIN  GATE  =  71110 

01  79 

P.  AT--C2.P,S<All3/CUN<t  1 

01  ac 

IF ( RAT.LT. i . 9-33  1GQ  TO  1002 

0 1  5  1 

COM  I  NO  6 

016? 

1002 

cunC4=cunc<-  -mks.mi  )  * r i  f  ?-:k/  >  cu;a 

Oi  c 3 

All  ^CA65(  r-8  ,  r  I  1  *PJ  i?»kp  )/  CUN';^  I 

01  P  4 

I  T  (  R  A 1  1  •  L  T  *  1  •  E -  C  3  •  ’.1 0  TO  100  3 

Cl  £5 

633 

CONI  INIJE 

01  6fc 

1  303 

Ci  NCI  K5  >  =C  UN'C  1  +C IJNC  2  ■*  CuM"  3 4  CUiiC 4 

01  67 

CUUL  I  K5  J  =CUr  C  lr‘5  )  /  i  ?  l  f.«M  C  '$*  t 

Oise 

ACONC -  C  AES 1 CH\C I  K  5  > ) 

Bits 

A,a  \c  =  cAtft / a1/  r  ( K5  > ) 

0140 

rCUNC*AT  AN2  (  M  »AGfCi!Nr  (  K3  H  'MtAL  (OUNCI  Ki:  )  1)*  180/ PIC 

0 1 Q  I 

TC  UNH=  rCVMC 

c:*2 

V  P  I  T  r  (  b  1 1  6  ) 

01-3 

V-«R  (T£ (6 , SC*  J  CUNC  (  K5  >  ,  *  CU^C  ,  7  CC  N'C  ?  CL!  N  0  C  K3  1  *  ACUhC ,  TCUU  C 

o:?- 

WR  J  7  El  ft  ,22  > KB  .CiJNC  !  ,:tNC2  ,CtMC  3  fCUNC  V 

Cl  9  5 

2  2 

rC-F.MC  TC2Xtl>/S2X,£lf.e.?X*bl3.3/M 

C  i  S  "j 

v  00 

CONTINUE 

r 

COMMUTES  CURRENT  DI5TR  fLTICNS  ON  THF  SPHERE 

G1S7 

UK( T£<6,?3) 

Rise 

2  3 

T  Ci  R  M  A  T  <////'  1 1  Cl  fFf  M  C.  1  S  TR  I *  u  f  1 0  A  A  L  f  NO  SPFEFICAl 

019  9 

C  L.i'  S  1  =p  UC  '-TANIHT-R  S)/2C 

0200 

CUNS2-C C. tC. > 

0  ?  3  i 

no  501  KM  1-  1  ,MAXX 

0  202 

C  LMS2=CUNS£«  n  l 2*KML)-Ai KMl >Mli/ i2»KMl*i2 *KMl-i ) ) 

0  ?  0  3 

501 

COM  INUt 

0  2  0*1 

CONS?  =  Cl.  Vl7  /l  COS  iHT-RS  1‘ftOAI 

C  2  3  5 

00  53  2  K  K  2  =  i  ,  N  •'•’  j 

o:  06 

CU3CKM2)*  RK(  KM2)  ♦C"UT?  \  \H2  A)/i  SX'  SX  Of  KL(  KH2  5  3 

02  0/ 

CUM  ^2  )  *P<  NM2  )  ♦CMiDl*  (KK-M  >  »!)Ni  (  1/  i£.X~S*  ^OPNL  ( K M2  )> 

0206 

5  0? 

CUNT  1  f-UE 

02  C  9 

DO  502  KM2-1  » 10 

02 1  C 

TF  (KMJ. EG- 1 0) GO  TO  504 

021 1 

I  TABS-  |TaCj*1CM  KM3-1) 

0212 

TAPti$*FMC*(  1  T  A  P  $  -  Gf  B  W 1  8  3 

C  2  1 3 

0FB$*fJF8 

0214 

GU  TO  5C5 

C215 

5C4 

ITABS--65 

0216 

DFeS^OFB 

021  7 

T  ARBS -P I C* ( 1 T ARS-CF BS ) / l 30 

0216 

50  6 

CO  506  KM4 - t , NMA 

021  v 

CKD*RTUM4) 

C22C 

CALL  LE CF  (  ITABS,DF6S*CKD,C,f  \NP ,  0  fNP  »  PNM f  DPM  ,  PNL  ♦  DPNL  ) 

0221 

50ft 

CONTINUE 

0222 

C  UW$  3  =  I  C  •  ,0. ) 

0223 

DO  507  K i ^ 5  -  1  t'VMA 

r>  22  4 

A  C2  -C  U3  (  KM5  )  *S  1  M  TAR  PS )  A;>2-  Cf’NU  KM5  I 

0  22  f 

CUIIS3  “C  ON  S3  <  A  03 

022  e 

h;-T*Cabscal:3/c  uns?) 

02 2  7 

I F f  R  A  T *  L T • 1 • t - 0 3 ) G G  TC  5C2 

C22  6 

537 

CONI  I NUE 

02  2- 

•5C  J 

ri  N  S  3  a  -  V  I  MG  *  C  U  N  S  3  /  (  ICC x;  R  L  T  A  H  > 

02  3  C 

CUNS4-I  C.  ,0.  ) 

">23i 

DC  5C9  K  ‘A6 ^1  •  MA X X 

02  3  2 

Jcr*t?^2*KM6 

023  3 

CIS-'-  i  j  •  ,  0  •  ) 

0^34 

r;0  510  KM7-1  ,  N ’Ni 

0  22  6 

CRD- R  F f  hM?  ) 

16/18/2; 


r'liCF  «*  ) 
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OPTRA,; 

IV  G  Lcv'rL 

It  H!N  CMC  =  71110 

02  2  6 

AD  3 '-CUM  Kfc7»  T  A ?  OS)*  *21  -C^IU  KM  7  3 /  f  K  \'l  I  KMi.2-11  .  GPO) 

07  3  7 

CUS4=CU<*«A!J3 

e  22  P 

IF  I  (fif-C-6)  .LT  .12  -K  r-V.  )  3  GO  TO  5 1 0 

0239 

PM- AHSi  A L' 3 / C 1  •  *•  «  ) 

G?4C 

I  F  i  K  A  *t  i  17*1  *  F  •  !  CO  *  C  511 

0  241 

510 

COM  !MI 

n  24  2 

SI! 

Cc*M  M  MjOS  Ai  M  6  Mi!  0  i  < 

32^3 

k  a  r !  =\,\£z  t  c  j£  •  -  .m  ••  ,  v  \ )  -*-v>  i  :o  •  <  ve.  >/  :  jn  34  > 

ir  <r>ui  .n .  j  .  >vj  u.  nz 

0  ?  4  5 

~C  c 

COWMMifc 

0 ::  H  6 

M2 

Cl’NSS  K  li  i  *c  UN5.i  C  JivSM  <’uMS4 

■:  ?  4  7 

C  UN  OS  <Ky*3  i  -  CUN  Si'  K  1  )  '(6  4:f  O  PS*S;?i<  r  !>S  )  ) 

3?- 8 
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3147 
C  14  8 
Z  L*-c 
015° 
3151 


>154 
315  5 
*'  13  6 
315  7 
0158 
Ci.  59 

0160 


0161 
31c? 
0162 
01  6h 
0165 

0166 
01  67 
3163 


0169 
0 1  7  C 
0171 
C  1 7  ? 
017  2 
0174 
.3175 
0176 
01  ?  7 
01  73 
01  7  9 
Cl  60 
0131 
Cl  C2 
313  3 
0134 
013  5 
0 1  H  6 
016  7 
153 
01  c  0 
01  90 
0191 
Cio? 

0 1  8  it 


lc(  TLSn-AUUlHtUS.lIS 

114  CALL  eESJ(*Sfc«GK£«*J.DttJ9 l.f-C/,  IFR > 

C  ALL  5 r  £  Y  t  ^  S p  *  W H  •  0  Y  •  T *5 Y  *  I  •  T  - 0 7  ) 

GT.’M  i  K  )  s-Ce^STBu1- '  H3  fK?3  ) /PN4  S  I  Kb) 

GTjl  iKIsSTPUMEJlKP.  J 
GTMCKi*STHLK-0Y«K3) 

dgtnl  (  k  )  -=s re*  m o .5-e\  < k ? )  *rtSP< nev (kb  ;  > 

!  r  (  I  TFT ,WC . 1 ) G C  TC  31 

CALL  BESYThT  fCwp  tf-V  ftfcVf  i  .  8  -  0  7  > 

Call  BrS J<HT ,lP j «  BJ  »CRJ • l .c-0  7*  7£*  i 
CTJ<  \  )  -  5 THU  *  3  J <  K !.*  J 

DGT  J  (  K  >  -  ST H!>«  f  C  •  5  *  B J  ( K 8  MH1  *rh J  OCP  )  « 

GTN( K>=STHU*3Yl KE ) 

0  C  T  N  (  K  )  -  .  TH C  *  DFY  {  KB  )  ) 

GTMf K)=~Cb *0 T  k >/BMOTN<  k ) 

OGTM(K» s-LB  OGT J i K) / EN+ OGTNl K 5 
31  G^IO  MUC  tK TG  )  =  GTMl  (K  )/STHO<) 

GJ10 j(K  fKTG)=GTJl (K  l/GTJ (K I 
GMlOHP(KtKTG) =G  THl ( K I /DGTrttK) 

GY  IS  T  It  K  ,  KTG  )  -  GTNl  iK  )  *6TJ  <  k  ) )  *  <  l-~  f  GTM<  X  )  /  G  TjW  I  (  K  >  )  *  (G  T  J  l<  K  )  /CTO  ( 

If  )  J  )  /GriX  JP1  c  K  ; 

GY1  S  J<:i  K  ,  KTG  )  -- (C>  TNI tK>  “OGT  J<  K)  >  *  (  U  -  (  DOT  M  K  >  /G  7M  (  K  >  )  -  (  G  T  J 1  (  K  )  /DO 
IT J( K ) ) } /GPMX Jl( KJ 
GO  TO  103 

115  CXRR  2=CXRR*CXRR 

UTRR  (1)  =CXRR*  "(XRKZ  )/?4 

UTRR  (2)  *CXRR2* t  41.  *CXKR2*  <— G?***  85*CXPP2  ;  1/  1 1  52 

UTKR  (3)  -CXRK'-CXBR?''  ,O0*37  5*CXKfc2*  (-36060:  *CXRR  2  M  765  746 -*2  5*CXfcft 
i 2 : n /41472C 

YTPfULI  =C*S**  < -9  7*CXKR2  1/24 

VTkP ( ?) =CXRR2 *< -12 5fCXRR2* ( S94-4  55*CXRR2 >  i/1 1 52 

VTPK ( 3 ) -CXk  ?*CaRR  21* { *  4 2 5425* CXAR 2  *  1*5 1 737 +OXRR Z * (-«6 ili 7 5 *4754 73*C X 
1RR2  )  )  1/414720 
DE  J*  R  =}  . 

DC  JPRR- 1  • 

D6YRK*1  , 

DEYPRR* 1 « 

J  T“  1 

DO  119  I  P P= 1  *  3 
JTt.-  JT 
Cb=ORu**  IPP 

D6JRR-0E JRP4UTRK ( IPP)/0° 

OCYRR^DEYRR ♦ JT  *OTKR ( IPP  >/CB 

cCjprr  or jfrp<vtrk( tpp)/oe 

0fc‘YPRR=rPYFRP4jT  tvTTR  UPP1/0S 
1 19  CCNT1 NUE 

oRMO*  6-RS*S  INHt  At  PHR)*C'£vpr  /CFYp  ) 

CXPNT - ( ALPHRP-ALPI  H )- (T  ANh( tLPHPP  >-T ANH ( ALPHH)  1 
EXP NT -c XPNT"OPD 

EXPfITB*  (  A  L  P HP -AL  PI  IPP  )-  (7  ANM  (  A  L  Phft  ) -T  ACH  {  Al  PHRR)  ) 

6  XPN  T  t  X  »>NT  R*ORG 
V2r-SGRT  (RS-5/HT  ) 

Y  5-  SGK7  (  T  A .NH  (  ALPHHJ  /TANH  (  At  phPR  )  ) 

/.]  OF  JKR/DE  JH-PEJP  ♦DEYRR*f:XP(-2i  Kc  XPMti )  /  (  DC  JH*DE  Yk  +  TBR  )  > 

GMlOMv  K  r  KTG)=Y?"  v:-  rfcxPl-rXPNT  )  !A 
G  JIG  JiK  >KTGl*GMlCP.I  K«KTG> 

G.1lOMP(  KtKTG  }  *GMlGPCfttKTG)4rHT*DEJH./A 

GYlSTi  (  K,KTO  )  =SCKT(  PSP*RS)-SCRT  f  TANHi  At  PHKT/7  ANH(  Al  P HkR  1  ) y* D t V R R  - E X 
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FORTRAN 


01SA 
C  1 5  5 
cise 

CIS? 

ci  s  e 


FORTRAN 
OOO  1 


CC02 
0003 
OOOA 
C  005 
0006 
0  00  7 
C  j  C  R 
0000 
OOLO 

con 
0012 
CC1  3 

oou 

0  01  5 
C01  t 
001  7 
CC1  s 
0015 
0020 
0021 
0022 
002  3 
C  0  2  4 
C  02  5 
CC2  5 

cc:t 

M2  E 
C-C  2  5 
0  0*'  0 
003  1 
C  CO  2 


IV  C  LEVEL  18  GUNN1  CATE  =  71110  10/13/30 

IP  (-EM  N  II  )/( 19.  M  ) 

G V  i  i-7 2  (  K  ,  K  T  G  )  =  C  V  1  ST  1  t  X  ,  K  T£  ) 

M3  CONTINUE 

GR’O-CPC'-O,  5 
RE  TOOK 
END 


IV  G  LEVEL  16 


CSIMSC 


CATE  -  71110 


15/ 16/06 


C 

c 

c 

c 


SUBROUTINE  C  5  I  MFC  (A,^,W,Nn 


THIS  SlioRCUTINE  SCIVES  SET  Or  N  LINEAR  ALGEBRAIC  cCUATluWS. 


C  CMP  LEX  Z  ♦  X  >  M  ^  t  3  i 
NPN!  =  N+  N1 
CO  5  I  =  1 y  N 

zno. 

cn  2  J  -  I  ,  N 

XI *CAKS< A( J.I  )) 

IK  (71- Xi*  1 9  2  « 2 

1  11-  XI 
U  -  J 

2  CONTINUE 

IF  (ZD  3  1 11 1  2 

3  Z  -  A  (  I  !  »  I  I 

Ai  i 1  ,  D  *  A<  I  ,  I  ) 

1  P  1  =  I  ♦  1 

DD  L  DP!  9  NpM 
X  =  AC  1 1 1 L  ) 

A  (  1 1  ,  l  )  ^  A  (  I  9  L  I 

4  A  (  l  ,  U  =  X  /  Z 
DO  8  J  =  1 » N 

IF  ( J-I  >  5.  8.  5 

5  IK  :REAL(A(J,ID  .e;‘.  0.  •ANC.  MMAG(  A(  J,  l)  >«ec.  O.)  CO  to  0 

2  ^  -A ( J  j  I  > 

Du  7  L  ~  I  FI.  MPN1 
7  A  (  J  f  L)  ^  A  (  J  « L  )  *  l  *  MI,L) 

Q  CONTINUE 
s  CONTINUE 

10  RETURN 

11  WR  J7 t  ( 6 , 1 2  ) 

12  FORMAT  (  31  H  COEFFICIENT  MATRIX  1C  .SINGULAR  ) 

RETURN 

ENTi 


FORTRAN 

9001 


CC02 
0003 
C  004 

ro*'  5 

C0G6 
0  C07 
00  OB 
0009 
0010 
0011 
003  2 
0G1  3 
0  C 1 4 
001  « 
CC»  6 
001  7 
0018 
r oi  9 
0020 
0  021 
00?  2 
0023 
00?4 
CC2  5 
00 2  6 
0027 
002* 
C02  9 
TCrC 
0  0  3  I 
00  32 
OC33 
GC34 
CC3  5 
C  0  3  C 
CC?  i 

*  r  rs 


FOP.TR-V! 

0001 
0002 
CC03 
0004 
00*5 
caws 
000  7 
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\v  »;  level  18 


MINI 


CATE  *  71 1 1 C  16/16/3^ 


C 

c 

c 

c 

c 

c 

c 


SJ3R0UT  I  NE  f; T  M  (  X  ,  T  A ,  f.'MAK  #  F  t  DERF  ,  FC  T  ,  >  ST  ,  tP  S  .  IEND»  tc«t  I  » R  T ,  DNT  I 


THIS  Sir.'RDJ’MNF  SOLVES  THE  CHARACTERISTIC  ECL'ATICN  Lv<9,)=0  FQP  4  FIX- 
EC  CONE  ANGLE*  -T  P*9'/Ii"P.>  ROTH  FOOTS  OF  L  vi  p.)  AND  CEP  IVATIVE  OF  MiH 

s  Fircc t  t:  cus*5 


DIMENSION  CT(FO)  f  CM  (  50  ) 
T A 2-14/2 


TClf-10C**F.F5. 

00  1  0  M  *  1  •  N  .v  A  >  ,2 

RTNC3179 

X  =  X  S  I 

RTNC3h40 

701  *  X 

CALL  FcrtTA  ,  rci fNl.FfOFFF  #FFf  n 
pn  t  JJatrlENr 

RT  h'j3«50 

lMFiM.l 

RT  NO  392  0 

l 

l Ff DERF  )?,?,? 

RTN03950 

? 

DX-F/OFR* 

WTH03980 

X - X— D  X 

KTN035S9 

\  01- X 

CALL  FCT(TA.TOL.Nl.F*i,f  RT #FP,!  1 

R7NC400P 

T  OL -E  PS 

RTNC4040 

A- A  3  S ( X  ) 

KTNC4f»*u 

IF<  A-i.  14  ,4.3 

RT 004060 

3 

tclmcl  +  a 

RTWCsu  70 

4 

IF  (  ABSiDxI-TOUE  ,5  ,6 

in  nosoum 

5 

!RABS(  F)-TeiF>?tT*o 

RING'- 090 

6 

CONTTNOF 

*.Trr 

7 

Kc(Nlil>/2 

RT{  K >  «X  ♦M 

DNTi K  )  *~2 *C OS (T *2  I *C0$ ( T A  ll  * FF/CERF 

IF ( ( 2 *-X) « IT. C*0C1 ) Or  TO  L? 

IF< X.CT.l .E ICC  1C  11 

xst  =  ifi  x<  xnccoo. )  / ia.ro. 

GO  TO  ICf 

1? 

XS  T  =  1 *ocqoOO 

10 

CONTINUE 
i  1  R  -  N  M  A  X 

RETURN 

1? 

IEP  *  VI 
i?  F  T  Ui>  N 

H 

ICk-rU 
,:Cs  URN 
ir  N  O 

IV  G  lcVEw  1.8  FCN  OATC-  71110  16/16/39 

FUNCTION  FCNU  .XN) 

T~ «  51*  X 

>  a-l.  570 79 6 3 7*  XN 

KH-.fi.it61r/76«-  <P5IG(XM1«  )  ♦ ALCG ( S IN ( T  I  I «  •  5772154649  >  -  T  A  N  { Y ) 

FCN- -FCN 
R 1  T  URN 
UNO 


-76 


FORTRAN  IV  G  LEVEL  13 


L  ECr 


DATE  =  71110 


0001 


C 

c 

c 

c 

c 

c 

c 

c 

/"■ 

u 

c 

r 

C 

c 


$ 10  ;i'  F  LET *  t  • ,C  ,PNF,OFNP ,FNP, FNL,  C?NL  > 


THIS  StbRCUTfNE  GENERATES  VALUES  FOP  FRACTIONAL  ORDER 
T IONS  * 

ANGIE-  IF  IX  f  &.  ♦*  J-CF 
VNr.y-CLR  rr  LcGENCPE  FINC7  !Cft 
PNP  =  P  v  (/-’> 
liFNP-  Fv?^  J  , 

PNM*PV  (-J‘>  e 
DPNM-Pvl-1*  I  - 
PNL=FNF-?N  P 
DPM  -DPNP-  DFA  A. 


0  00  2 

DI^FNSICN  TEFKl  (50)  ,TEP H2(50)  ,SUH1(5 
]  DP/iP  (  100  >  ,UPNM  ICO)  ,PNi.  (50)  ,  CFNL  i  50) 

0003 

PJC**3*1<1592E5 

CCOA 

X  =f  HA-DF  >  *  P  I C  / 1  5  0 

CCC5 

C  X  P  ®  C  G  S  (X) 

cccs 

CXK--—CXP 

CCC  7 

K  =  <  IF  IX'XNUl  )/2 

ccce 

C=5«  E-05 

0009 

N= IF  IX ( X N 4 C ) 

001  c 

FNaXN-FLCATCM 

0011 

NNl-N+1 

C Cl  2 

XO-F  N 4  1  * 

cc:  3 

XI  =  X0*  u 

CCI 

2 

CALL  TERHSiXvFA/»DfTLkHl«  iLRH2f2» 

OOL  5 

PNP ( 2) ~C. 

001  c 

PNPI 1 )  -  C  • 

001  7 

ID  ;.C  11  =  U  I 

oci  a 

PMrmaFNFCZMTERMl  (III 

C Cl  9 

PNP(  3) -PNP43) 4TL  RM2 ill' 

0020 

30 

CONTINUE 

0021 

6  C 

CALI  SUPS l F  N , I iS  UM1  ?$U!*2  } 

0022 

4 

SMO~  SIN  5  XC*P l C  > 

002  3 

SN1-  SIAM  XI  *P  IC  ) 

0  02  4 

SU1  =^C. 

002  5 

SU 2  - 0* 

00?  A 

L»L  SO  1=2,1 

002  7 

SU3  =  S!)l4  SNCMTERPH  U-SUHl  <L>  l/PJC 

0C2  5 

$U2~SU2*SF  1  v {  T E R V 2 (  L  MSUM2CL  )  )  /PIC 

C  C2  9 

5C 

CCNT INUE 

C  03  0 

PWM(  2>  =  (  i  .-SfiC^-Crit  X,  >01  >*PNP< 2 >*SU1 

0031 

K-.M  3  )  -  ( 1  *  - SN 1  *  F  C  N (  X  ,  1  )  )  *  P IJ P  C  3  )  ♦  SU2 

C03  2 

PNP  U  )  =  (  i  2.  xcn.  l^Cxp^PM*  (2  )-  C  /O*  l . 

C  03  ? 

PNH ( 1 )  =  ( ( 2» -  XC4 j  -  ) *C ( 2) •( XOvl • 

COS  4 

C  X  P  2  *  C  X  P  *  C  X  P  - 1 

3C35 

CXM£rCXF^CXK-l 

0  03  6 

DPNPl  I )^XC -  I PNF(2)"CXP^PN0( 1 ; ) / C  X - 2 

003  7 

PPNMO  )  =  XC  M  PNNC  2  Wlxf^PN*  ,*  i  )  J/CXS2 

C03  S 

Dr-  MP  i  2  )  -  *(>(C  XP*  PNpt  2  )~*>NP  (  1  )  >  /CXF2 

cc3  9 

DFNN(  2  )  -XC*  I  CX^-FNfi  v2  >~P’,n  (  i  )  J/CXM? 

C040 

DPNP  (?) “X 1 •  i C  XF*  FNP ( 3 )-°N»  i 2 1) /CXP* 

00  41 

OPNfM  3)  •  xi  -  (  C  xr~  -  PNM  (  ?  ;-PNr  (  2  >  i  /CX‘-2 

0042 

IF  ; f;l\ J  - 4  11  0 0  *  5  ,5 

004  3 

5 

DO  4  C  LK='*tNM 

0)YSUft’(50)«PNP(100li 


)  *PNP(  5  )  )/X0 
) +PNMI3 > )/X0 


ie/ie/19 


UGENCRE.  FUNC- 


PNMtJlOO  > , 
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FORTPAN  IV  C  LfcVEL  1 8  LEGF  DATE  =  71110  18/18/39 


0C44 

KKl-KK-1 

0045 

KK2=  KK-2 

CC46 

K  K3  =  KK- 3 

004  7 

XX=i(0«FL0AT«KK?) 

c:4  e 

XXD*X0»KK2 

0C4S 

PNPIKXI  =  H2.*XX ♦!.) *CXP*PNP I KK1 »-XX*PNP(KK2) I /tXX+l) 

0050 

PNMI KK)= ( ( 2.«XX«1* 1  4CXX*PNMIKK1 > - > X* PNN <  K K2 )  )/( XX  +  1 ) 

CC  5 1 

DPMP  I  KK  1  ■* > XDMCX  P*  PNP (K  K  I-PNPIKK1  )  I/CXP2 

0052 

DFMMtKK )  =  XXD"ICXM*PNMiKK)-PNMim ))/CXM2 

005? 

40 

CONTINUE 

0  C  54 

ICO 

PNL ( KI-PNP I NNl ) -PNMfNNl ) 

005  5 

DPNL  (K)  -DPNP  (  NN’l  J+CPMIMNNl  1 

0056 

PNL ( K)-C.5*PNL(K) 

CC57 

DPNL I K) =0.5»DPNLIK) 

005  8 

RETURN 

CC5  5 

END 

FORTRAN  IV  G  LEVEL  18  SPYNX  CATE  =  7111C  18/18/39 


C001 

C 

FUNCTION  SPYNX IN. X) 

c 

c 

THIS  SU8R0UTINE  GENERATES  VALUES  FOR  SPHERICAL  FUNCTIONS  OF  THE 

c 

r 

SECOND  KIND. 

L 

c 

0002 

IF  IX  •  G 7«  0.  )  GO  TO  1 

0003 

$PYNX=0* 

C004 

WRITE  (6,2) 

o:o5 

2 

FORMAT  1 78H  SIR*  SOME  THING  WRONG  WITH  YOU  ...  YOU  K  NOK  THAT  THE  A 

1 RGUMCN  T  SHOULO  NOT  BE  0) 

C  CC  6 

RETURN 

0007 

1 

IF  IN  *  E  0«  1)  GO  TO  3 

0  CC8 

A  =-C  C  S  (  X  )/X 

CCC  5 

IF-  IN  •NE.  C)  GO  TO  3 

0  310 

SPYNX* A 

0011 

RETURN 

0312 

3 

B«-(COS( X)+X*SIN!X) l/X/X 

0013 

IF  ( N  • NF •  1 )  GO  TO  4 

0014 

S  P  Y  N  X  -  •  B 

CC1  5 

RETURN 

001  6 

4 

U=  1  • 

CO)  7 

00  5  1=2, N 

0018 

U=U  +  2 • 

CC1 9 

C  =  U*3 /X- a 

0020 

A-B 

0021 

5 

B=C 

0022 

SPY  NX  =  B 

0023 

RETURN 

0024 

END 

FG RTR-VJ  i  /  o  LEVEL 


BIS'* 


DATE 


71110 


1S/1S/39 


0001 


0002 

0CC3 
G0?4 
OC35 
0  OC  6 
0007 
CCC8 
OC04 
0010 
0011 
0  2.12 
C  01  3 
0014 
0015 
001  6 
001  7 
0012 
001  <? 
0020 
0C21 

0023 
C02  4 
002  *3 
002  6 
002  7 
0  12  8 
CG2  <5 
0030 
0031 
C03  2 
0023 
0034 
OC35 
OC?  6 
00*3  7 
CO?  8 

CC3S 
C  C  4  c 

0041 
C  0  4  2 
1043 


SUiM.GUT  IKE  PT  SV(X  **  :  »  n\  .CbV  ,1 ) 

c 

C 

C  T  H  j  f«  Mjfc3CLT»Nr  Ot>.r&ME$  F**  4  T  IZfJAL  6  iSSLl.  S  rUNC  7  10%  Or  THE  SECOND  K 
C  rOk  45-  GU*lcM  Lt’.S  THAN  S 

C 

C 

DIPT  NS!  CN  L3  Y  C  *:  0  )  ,CB  t  (  ?C  )  *H  VT1  (  50  )  ,BYT2  l;0)  ,  3  JT  1  (  3  1  >  •  FJT2  (SO  > , E JP  (  1 

i  o )ti  Jv(  i c  i  • «  v  i  i  s c  > » ?. vr?  < 3 o ) 

PI  =  *  -  l4;5S2t~ 

c-i  :  -  o : 

( 1  !  X ( X  N*C ) 

N  1*1.4-  1 

f N“)N-FICAT(N> 

I H  X • G E  * 5 • 0 ) G l  TC  ICO 

C4U  £Y1  FfcP  l  >  ,rW  #0,3*71  »  3Y7?  ,VM  ) 

CALL  C J T  E k  V ( X  » f  N  f  0 1  2  J T 1 V3JT/ »PNI 
D  JP ( 1  )  =  C. 

B JP< 2  )  =  0- 
B JK( 1  )  -C. 

3JP<21*C. 

XTH  =  f  X / 2 •  )  vS  * F N 
XTFM  =  l  X  /  2  ♦  I*  *  Q  .4FN  ) 

00  40C  1 K  =  1  ,f'N 
3  JP  I  1  J  *  3  J  P  <1  i  <  >J  T.l  (  !  K  > 
bJPt 2)  =  BJP( 21 *  J  T  2  t  l K  J 

400  CON' INUF 

BJPll)  *XTFN*fcJ  IT  1  ) 

>  J  1  i  1  .1  \i  •'  i  .*  # 

DO  4C1  I H  =  1 ,  P  P 

EjrJ  1  HEYTIUF) 

e jrt mpyt2( 1 h > 

401  cc:r.  IUJE 

(1J-4M  l«KJH  1  !/Xl  T N 
3Jtf(2  2I//.TFNT 

by  (*  imejpci  kccsif.\*pi  i-«jNm  )/s  im  fn*p  n 
6 y(  ;  j  - 1 2 jp ( 2 1 4crs ( (  i.*fn>*p  d-cjm(2)  )/su;(  ( i  .♦fnmph 

CPV(l)a-8V(  2  >+PN*BY(l  l/X 

i- 0Yi  2  )  =E  Y  ( 1  )  -*  (1  ,  +  FNI*m?  )/X 

!F(r.'l-LE*2l  GO  TC  101 

Dr;  200  P=3 , N 1 

OP.U!  ^  ^LCAT  1  V-l  MEN 

0*1 2=CkCl-l . 

3Y<  i-U  =2  .  ▼  CR02*i5Y  (  F-l )  /X-CYl^-2) 

DriYt  M  )  -  *1 Y  <  H-  1  )-CPDl*EY(.M)  /  X 

200  C  C  NT  f.VJE 

10} 

I  CO  C  /  L 1  BE  $  Y C ( X , XK «  £  Y • £  S Y  »  3  »  I  E  P  ) 

FT  MyF-r* 
tfiD 


TORTPAM 

0031 


0002 
0003 
0004. 
OC  05 
0006 
cro? 
ccca 

0039 

OCIO 

bon 
0012 
0013 
001  A 
C  Cl  5 
C  01  6 
0317 

ooie 

001  5 
C02  0 
0021 
0022 
C C?  3 
0024 
0025 
0026 
002  7 
0028 
CC25 
00  3C 
C  3?  1 
0032 
003  3 
032  H 
0C3  5 
0036 
0C2  7 
0038 


FORTRAN 

OCCl 

0002 

0003 


0004 

CCC5 


FORT  RAN 

0001 
C  CO  2 
C  003 


0004 
OCC  6 


-79- 


IV  G  LEV  11  1?  fc  f  5  V  G  DATE  *  71113  1  8/1  8/ 29 

SIB ROUT INF  Bf ^Yn(X,XKiPY»E&Y»C« It*) 

c 

c 

t  THIS  SLBRCUT  INF  Cf.Nt' RATES  VALUES  FOR  FRaCTICNAl  ORDER  BESSEL  S  FUNtT- 

C  TICNS  OF  THE  SfCCNC  MNP  WITH  ARGUMENT  GREATER  THAN  5. 

C 

C 

DIMENSION  BY i 50  J  » CR Y ( 50 ' r  PY  c 1 ( 50 i #  BYE  2 ( 50  ) 

PI*  2.  141  59.765 
C=1 .E-05 
N* I r I  X ( X  N*C ) 

Nl=rj+l 

F  N -  A  N  -  F  L  CAT (N> 

CALL  BYTE  KG { X .PM w BY  El  tfiYE? .MA) 

XEI*l./f  2'-X) 

X  E  1 2  =  XE  l  *  *  2 
6  Y i  1*1* 

BY)  2=  B YE  1  12  XFI 
B Y2 1  * l • 

B  Y22  =  BYE  2 ( 2  )  XE I 
JT  =  i 

DO  400  IK-I.MA 
JT  =- JT 
I K2  =  2*1 K 

BYl 1-BYll *  JT  "BYE  I i  l K?*l I  * XE  12*  *  I K 
BY12*BY12  + JT  -BYF1 ( 1 K2>2 ) *Xfc  I *XE  I2«*JK 
BY2  1  =  BY  2 1  ♦  JT*6YE2(  IK2*l J*X£  12**IK 
BY22  =  BY2  2^  JT*CYF2( TK2*2 ) *XE I*XFI7  >*IS. 

400  CONTINUE 

AG  1  * X-C.  25M 2-FN^l) APJ 
AG2-X-0*25M2 v(l*FN)< 1)*PI 

8YU  )  =  <S1N(  AGU*  BYl  i*COS(  AG1MBY121  *SCRT  <2./<PJ*X  >1 
8Y(2)~(5TNfAG2)*r  BY2  1  *COS  ( AGO  )  *BY22  )  'SGRT  I  Z-  /  l  f> I  *  X  > ) 

DBY ( i>=— BY(2  )  ♦  F  N  *  BY {  i  )  /  X 
DBY(Z  )  =G V ( 1 )-(l**FN)»BY(2)/X 
IF  C  Ml •  LE • 2  )  GC  TC  101 
DO  200  «M  =  2  f  N 1 
OP.Dl*FLOAT(M-l)+FN 
0RD2=0RD1-1. 

£Y(H>*2  »  * ORD?  *  BY ( K- 1 >/X-BY(M-Z) 

D8YC  M)=BY(H-1 ) -ORDi * BY  I  M) /X 
200  CONTINUE 
101  RETURN 
END 


IV  G  LEVEL  18  CPSi  DATE  -  71110  18/18/19 

FUNCTION  DFSI  <Z) 

X- Z- 1 • 5 

DPS  I  *0.9  34  90200*  X* t-C*  E2  8  7  9  30  «■  X"  i  O  *  704  502  94*  X  57956094*  XM 9  *4 

1  642 634 5 *X*(-G. 35206372*  X* (0.26484 460  +  X* (-0.27630538*0.21 785379*X)  > 

imn 

RETURN 

F.NQ 


IV  G  LEVEL  18  PSi  DATE  =  71110  18/18/39 

FUNC 1  /UN  PS  I  i 2 > 

V  -  Z  -  1 «  5 

RSI  =  0  *  3  6  4  F 5  5 1  41—01  *  X  *  ( L*  *  93*  BC09?  *X-  (-C.4  )  439B37*  X*  .'0. 23.4*22  7  7*  X*  L~ 
10.  144  79SC4*  >>  i  C.  911  ?>  -  77 1  -  01  tx  { -0 . 594Z2363E.-01  ♦  X*  (0 . 49350996E ~0 i- 
j  0-  32fcfrl  5S2L-C1  X  )  )  l  I  )  ) 

RETURN 

END 


80- 


FCRTRAN 

0001 


C002 
0002 
000* 
C  CCS 
00^6 
CCC7 
GCCc 
0C09 

cc;  c 
OCil 
0012 
00^  3 
GC1  A 
CC1  5 
0C1  £ 
001  7 

cci  e 

CC1S 
002  0 
0021 
0022 


FORTRAN 

OOOI 

ooo: 

0003 
0004 
0  0  05 
0006 
CC  07 
0000 
0009 
0010 
001 1 


FORTRAN 

0091 
OC  02 
0  0  0  3 
COCA 
CO-OS 
cr  06 

C  CC  7 
OCCR 

coc« 
con 
ooi : 


IV  G  LEVEL  1«  eCSJG  GATE  =  711 1 C  16/18/39 

SUBROUTINE  BESJGO,XNfejT) 


c  this  Subroutine  generates  values  for  fp.act  icnal  order  cessel's  funct 

C  T1CNS  OF  TFE  FIRST  KINO  vs!TH  ARGUMENT  GREATER  THAN  5. 


DIMENSION  &JC1C50 > 

Pi-3*  141  55265 

oi«  =-rs 

\-  I  -  !  X  (  X  N  4  c  > 

FN-X \-F LCAT ( N ) 

CALL  B  J  T  E  R  G  C  A  »  r  N  |  F  J  E 1  ,  N  A  ) 

XEIM*/(2*X) 

XEI2-XF 1**2 
B  J1 1- 1* 

BJ12-BJE1  (2  )  *  X  E I 

JT*1 

DO  *0C 

I  K  2  -  2  *  I  K 

3J1  l«BJmJT*5JEl  l  UiU  ;*>E!2**IK 
B  J12»BJ12+  JT*eJEl  i  IK2*2  )*XEI*XEI2*  *  IK 
4C0  CONTINUE 

AG=X**C«  55M  2**e\+l>  >  F I 

e  JT=  I  CuS  (  AG  D9  J1  l-S  IM  AG  I  *r3J12)*$QRT  {2./(f  I*X  I  ) 
RETUf  N 


IV  G  LEVEL  16  PSIG  DATE  -  71110  18/18/39 

FUNCTION  PSIG(X) 

IF  (X-2.J  10,11,12 

10  PSIG-PSI ( X) 

RETURN 

11  PSIG'- 3* 4227843351 
RETURN 

12  l  =  !F  IX( X )-l 
Z  =  X-F  LC A  T ( I ) 

PS  I G- PS  KZMSUMPC  Z-l.tl  > 

R  ET  URN 
END 


IV  G  LEVEL  1?  CPSIG  D AT c  =  71110  13/18/29 

FUNCT ICN  CFSIGIX ) 

IFix-2.)  13,11,12 

10  UP5IG-0PSI (X ) 

R  ETLR  N 

11  D  P  5 1  G  -  0 »  6  4  4  9  3  4  C6  fc  8 
RETURN 

12  !*£FIX(XI-1 

7 *X~F LCATI l ) 

OPSIG-OPSK  Z  )  -SUPER  (Z-l.  ,  I) 

RETURN 

EDO 


-81 


FORTRAN  IV  0  LEVEL  16  t'VT-RG  DATE  =  71110 


0001 

S  1.8  P  OCT  IKE  b>K 

•j r, :  •<  ,•  n  ,8v ei  .evri  ,ra) 

0002 

D I  V*S  vSIC  N  rfYfl  ( 

0003 

PA*l  *13/ALCG(X  5 

CCC4 

% e I  1 1  t  =  i. 

0  C  0  5 

BYE?,'  1)  =1  . 

OCO  6 

OCC  7 

CO  1  :c  I  3  2  #  v.*'  2 

OCO  6 

11*1-1 

CC0  9 

I  ?-i  -2 

001  0 

xu^ v'ati  n  > 

CO  L  1 

XI  2=  LCMU?! 

0012 

8  YF1 < II *6 Y£1 { 1 1 

)  *  (  30  ♦FN-XI  O  «(  J  .O+Ff  •:  XI?  >  'X  1 1 

001  3 

6  Y  E  2 1  I  )  -  B  Y  E  2  t  ll 

)Ml.f>*KN-XJ  1 )  *1  1.5+FFMXI?  >/XU 

0014 

100  CONTIMJc 

CC1  5 

RETURN 

001  6 

E  NO 

FORTRAN  IV  G  LEVEL  IR  t?JI  6RR  GATE  =  71110 


0001 

SU3RGLT  I  At  6  jT  E6  R  !  X  ,  FN  ,  V  ,  61  ,  R,\  ) 

0002 

DIRE  NS! CN  67(60 

0003 

I  IX ( 1 -6»>*1C.  ) 

000  4 

CALL  GRFR.AI  J  •  «ri\,GX,  ISR) 

CC05 

Xf  2  =  I X/2. )**2 

0006 

JT=1 

ooot 

ptti  »  *  .  •rr. 

uCl  o 

u  u  vis  -*  y  r  yj 

OCO? 

JT=-  ,J1 

001 J 

PNi*=rN-l 

0011 

XMN  =  PlCAT  <*M 

0012 

X*\l  r  LCAT  (  M N  1  > 

0013 

6T{MfO-J7*f3BS(BT(PNl)>*XE2/(X>M*(  XMNMFN  )  ) 

0014 

IF(Aii$(tiT<MN)).L£.0)  GC  TC  10 

0015 

20 

CCNT I KUF 

0016 

10 

P.ETUf  A 

CC1  7 

EKO 

FORTRAN  IV  C  LEVEL  16  eJTERG  DATE  =  7111C 

0001  SLER0LT1NE  P JT! RC I X  ,  FN, L«JE I  ,M A ) 

0002  OIFENSJCN  3  J  E  1 ( 5 C  > 


0003 

F*A=-1  «:0/ALCG(X> 

CO  3 4 

bj:;  <  n-i - 

1PPU 

?♦? 

C  0  0  6 

DO  100  I  =2,  M2 

0007 

13-1—1 

002* 

I  ?  ~  I  -  2 

OCO? 

x i:  - n t ( ii  > 

0010 

X  ]?.*?  LC.  n  (I?) 

0311 

BJrli  I  >  =  2  J  e 1 C 11 

0012 

100  C  C?N T  IMT 

0013 

-  E  T  l.V  .\ 

0014 

fc  NO 

* (0. S+FN-X! 1 >*<  0.5+FN4XI2  »/Xl 1 


15/13/35 


18/16/3° 


1 8/16/30 


-82- 


FORTRAN 

IV  G  LEVEL 

13  PACK 

0001 

Ftirjt.T  !C>  FACKIKI 

000  2 

FACIEI. 

CCO  3 

00  12  J  =  : f K 

000  A 

xj- j 

000  9 

13 

f  AC  K  =•  f-  A  C  ft  *  <  J 

CCC  6 

RETURN 

CCC  7 

fi  NO 

FORTF  AN 

IV  G  LEVEL 

IS  PACK/ 

0301 

FUNCTICN  FACKP<X,K> 

CCO  2 

<  AL ►  P  ”  a  e 

COO  3 

1  F  {  K  *  L  £  •  1  »  <>C  TC  12 

0304 

K  1  =  K  -  i 

CCC5 

DC  11  I  s 1  , F  i 

o:c6 

XI  »I 

CCC  7 

p AC KP  =  F  * C K F * i X*  *2 -X : * *L  i 

ocoe 

11 

(  0>‘T  I  N'JE 

3C0  c 

12 

*-ACKP«FACKF*X*iX*KI 

COIC 

t*.  v.  T (j  ft  a 

0011 

c  NO 

FORT  ft  A’'! 

I\  G  LEVEL 

18  P  CTO 

0001 

FUNCT  IC<\  F  C  T  C  { T  A  t  X  #  M.  1 

0007 

y»'  v  4f  *  * 

W  b  V<  i 

r-x.  5  f  J  /v* 3  J*  O’  xl  0  1 

GOO  <* 

4-T  a  /2 

0005 

FC  TO  =T  AN  <  Y 1- C • 63  97  6*  <  ALCG  <  $  IN  t  A  II 

0006 

rtf  TU7N 

CCC7 

ENG 

FORTRAN 

IV  G  LEVEL 

18  CFCTO 

0001 

FUNCTION  CFCTC(TAfXtM) 

OOC  2 

X  N  -  X  ♦  M 

C  303 

Y  - 1  <*  57079  633*  { X*  1«  ) 

0304 

A  =T  A/ 2 

CCG5 

Dr  CTO-O.  M  6  61  9  7  6*DP$  I  CIXN+l.  1  ♦1.570' 

CCC6 

R  t TURN 

0007 

FNO 

FORTRAN 

I V  C,  LEVEL 

U  SUMP 

ocoi 

FLNGT1CN  SUFHXtKI 

000? 

SU.,#P*C. 

0003 

on  13  J -1  f!< 

CC04 

X  !  -  J 

0005 

ii 

$U*‘P-S*J>P<-  I  •  /  (  X  ♦  X 1 ) 

OCCt 

F  F  T  0  P  N 

OOC  7 

END 

date  =  71110 


'CAT 


S/**  *  - 


71110 


DATE 


7111 C 


DATE 


71110 


DATE  =  71110 


18/18/ T9 


18/16/39 


18/18/39 


\  ) 


13/18/39 


13/18/29 


“83 


FORTRAN 

OCOl 
0002 
CC03 
CC04 
0005 
C  006 
0007 
0008 
000  5 
0010 
0011 
0012 
001  3 
C Cl  A 
0C15 
0016 


FORTRAN 

000  1 
0002 
C0D3 
0:04 
OCG  5 
C0C6 
CCC7 

occe 

i  no  *- 


FORTRAN 

C001 
CCC2 
OC  33 
0004 
0005 
0  Cn  6 
0307 
CCC8 
0009 
001  0 

oo :  i 

0712 
CO  1  2 


FORTRAN 

0001 
C  002 
CCG  3 
0004 
COO  5 
0CC6 
noc7 


IV  G  LEVEL  1*  SUMS  DATE  =  71110  10/18/39 

SUNfOUl  INE  StfS  i  Pr..  I  *SUH1,SU*«2I 
Cl:- ENSIGN  SUV i I  5 7 > . SIP? «  SCI 
X0---1  « 

X  3  =XD>K 
S U  M 1 (  11 =c« 

su>°( 1  >=•:. 

ON  2  0  "  2  •  1 

01  = 0-1 
J-2 

FJ.1  =  :LCAT(  j») 
r  J2--31C-AT  ( jt?> 

SL.  iiJuSwHUli  - ...  /f  Jl  «  1.  /(  XO  +  ^Ji  i- :  •  /(  v>FJ2  ) 

SL”  2  »  J>  aSIKt  1  v.D-2.  /FJl  +  l./<  Xl»FJl  1-1  #/(Xl-pJ2) 

20  CCM  I  MJ  6 
RE  TO  N 
EM' 


IV  G  LEVEL  13  DATE  =  71110  15/18/39 

FUNC  r  ICN  SUPPCX  f  K  ) 

SuKM«0# 

K1--K+  1 

00  1  I  J - 1 »  K 1 
X  J 1  =  J  1 

13  S  U P ! 5UF  P  ♦  1 «  /  (X-X  J1  i 
REXI'NN 
t 


IV  C  LEVEL  1$  SUMNS  DATE  =  71110  18/lo/25 

SOUR  OUT  INE  SUi'NStTAfXiNlt  I  »  SUMN  ,  5UMN  f; ) 

DIMENSION  $UPM5C)tSUPNC(50) 

PIC- )• 1*1592*5 
XN-X  *.\1 

rc -f*dC(TAf  XtNl  I 
CFC-CFCTC( TA  ,  X , M  ) 

00  2  0  H  = 1 1 1 

P1=P-1 

SUMN;  P|  -{FC*PIC«  2*5UP(M1-S0PP(XN,P)*SUPK(  XN,  PI  )  ) 

SUKMWP  )  r(SUPP(XN?Ml  >  ♦  SU  M ,°  {  X N  t M )  )  *SUPN  CM )  +  C  0  FC*P  l  C  +  S»JP  FP  (  XN  ,  M  ) -SUP 
1PP<  >:?.  f  f  1 1  ) 

20  CONTINUE 
RETURN 

tfv  J 


IV  G  LEVEL  I*  SLP  DATE  *  71110  18/19/3S 

FUNCTION  SU>'CK> 

$U>  -o. 

DO  13  J*I,K 
x  J  -  J 

13  SUP ~  SL w  *  I . /X J 
"  n  t  ■  k  n 
t  K!> 


84- 


FORTRAN 

COOl 
COO  2 
0003 
0C04 

C  CO  5 
0036 
0007 

occe 

c  ccc 
0010 
0011 


FORTRAN 

0001 
CC02 
C  CO  3 
0004 

000  5 
CC06 
C  C  0  7 


0  CO  9 
0C10 
0011 


FORTRAN 

0001 
COO  2 
CCC3 
0004 
0005 
000  6 
CCG7 
CCC8 
0C09 


FORTRAN 

0001 
0002 
CCO  3 
0004 
C  C  05 
0006 
0007 


IV  G  LEVEL  18  EFO  DATE  =  71110  16/18/39 

fumct i c; n  erc(  xi 

if  (  x-3 •  i  l ,1  ,2 

1  Y«X«X/9. 

ft  10  =  l.*Y*<-2.2  499  997*  YM  1.  26  S620  84  Y*  f  -  .  2  1  63fl  66*Y*  (  . 0 4444 79* Y*  < - . 0 0 
13344 W.0CC21*YI) I  II 
RETURN 

2  Z=*./X 

r  0*  •  79 76  £4 5 o ♦  Z * (  -•  7CC0C  3  774-  Z*  < 005 5274*Z * (  -  •  COC 0951  2  ♦  /*  (  .001  3  72  37 
00 C726C5*.CC  .14476  1  1  )  1  ) 

TO-X-.  725  2P?i  £  */*■  C4  U6  357* /*  (- •OOOC?.5:>4  «  Z*  I  .0026 2  S  7  3 «  Z«  i  - .  000 5  a 

ii?o»  :•  (  -  «cc  : ; ; « .ccc:  25cs*z  >  >  >  > : 

DFC- -0*C  CS  <  T3 I /S CRT  (/ ) 

«pTU7A 

END 


IV  G  LEVEL  18  8F1  DATE  *  71110  13/18/39 

FUNCTION  BFl<*> 

IF  (X-3. )  1,1,? 

1  Y=X*X/9. 

6F 1-  X*(  •  5  +  Y*  (  -•  5  6249  985  OM  .  21C93573*  Y*  C 03954209 *Y*  (.  C044 3  3i 9+ V* 
l  (-.  0  3031  761  +  Y+.0CC011C9I  >111  ) 

RETURN 

2  7--3./X 

f  1^.  797664564-  Z*  (.0300C1 56*2*  (  • 016 596&7+Z*  ( .0C01  7 1 05  ♦  Z*  (  -.0024  9511  ♦ 

1  /  *001  ivs^.r.:r;ro?*/)  >  > >  > 

*  >  •  f  «  ^  c  *  c  **  t  T4  ^  t 

1  8  ♦  Z  *  ( . DC C79B24-.0CC 29166)  ) ) ) J 
SF1=-1*CC3(T1>/SCRT(X> 

RETURN 

END 


IV  G  LEVEL  IP  SLMMM  DATE  =  71110  13/18/3° 

FUNCTION  SURFR(X,k) 

sur 

K 1  -  K  ♦  1 

DO  13  J«lfKl 
Jl= J-l 
>J!  -J1 

13  su*?Vi.-v-:i*F,H4  l.  /<x-XJl>**2 
R  1  1  ijo  N 
C  NO 


IV  G  LLVGl  H  SUMPP  DATE  *  71110  18/10/39 

FUNCTION  SL»V  r  F  (  X  ,  K  > 

SUN'1>p  =  C 
CO  '  0  J  *  1  •  J< 
r.l'J 

i:  surff*s ./<  x4xim*2 

F  C  1  U*  N 
END 


0001 

CC02 

0003 

000  A 

0205 

0006 

0CC7 

C  006 

CC09 

C01  C 

0011 

OC  i  2 

001  2 

001  4 

OCX  5 

001  6 

CC1  7 

COl  6 

OCl  9 

CC2  0 

0021 

002? 

0023 

On2  4 

002  5 

OC2  6 

0027 

oc2e 

00"  7 

0030 

:ORTR 

0001 

0002 

0003 

0004 

0CC5 

?CC6 

CCC7 

0006 

0009 

C010 

0011 

OCl  2 

001  2 

0014 

C  C  1  5 

0016 


-85- 


IV  C  LEVEL 


4 


6 


10 

15 


50 

60 

62 

64 

70 

ec 


no 


7  20 

120 


18  CMP  A  CATE  *  71110  18/18/39 

SL3P0LTIM  Gf  P.'*A  (  XX  .  CX  ,  I  ER  > 

!  F  (  X  <-r 7  .  )6  ,t>  ,4 
IER  =  *> 

G  X-l  •  E  ♦  7  5 

RETURN 

X^XX 

ERR -1*  CE-06 
IER^O 
G  X=  I  . 

I  F  (  X-2.015C,  if.lii 

ire  x  -2. c  mr-,iiOfi5 

GX=GX*X 
GC  TO  10 

if- ex-i.o)  e.c,  nemo 
IFCX-ERf.  162,62,80 
/-FLOAT ( 1 M  C  x ) i -X 
I  M  Ai31(  Y  I-Hf  i  ;  *i  C  T  1  00,-4 
I  f-  (  1 • 0- Y- E R  R  >  J  :  0 , 1 0  C  ,  )C 

Irrx-.U0)3C, 80,110 

G X-G  X/X 
X  —  X  1 « 

GC  TO  VC 
Y'X-IO 

GY  =  1  .0<  Y*(  -0.5  7  7*.  Cl  f  t+Y^  IC.9E  58  5  3  99  4  y  *  ( -0 .6  7  642 1  8  7  iYM0.6328  212<Y  + 

1  (  -  0.  5684  7  29«-Y*(O./54£2C4  9*Y*(-0.0514c»93)  )>>>>> 

G  X  -  G  X  *  G  Y 
RETURN 
l  ER*  1 
F  •'  I‘uRk 
FNO 


IV  G  LEVEL  18  FCT  CATE  =  71110  18/18/29 

SU3P  JLiTt*F  FCT  ( T  A ,  X  ,  M  ,  F  f  CE  RF ,  FP ,  I) 

D  I  HE  \  S !  C  N  TERFN ( 50 )  ,  SLIMN  (  50  ) ,  SUMN  0  (  5  C  ) 

P 1=3.141 39 265 
FP^O. 

CALL  TFRFNSf T A , X , M , T ERFN  ,  I) 

CALL  SUVNS ( 1  A ,X «  M ,  I , SUWN , SUMNC ) 
f  =f*  C  T  C 1  7  A  ,  X  ,M  ) 

DLPr =GFCTC(  TA,X,Mi 
DO  4  1  K*l,  I 

P=MTE.?  PM  y  I  *SUriMM)/PI 
DL>F=  CEPF  *T  £RWN(  F  )  *SUFND (  F» I  /P  I 
r  P  =  rP-(TERM«(P)*  C  P*S'JPN  «  F.  I-  1  )  ) 

40  CONTINUE 

FP-(  1  .♦FPl/PI 
50  RETURN 
END 


FORTRAN  IV  G  LEVEL 


IS 


TERMS 


DATE  *  71110 


18/18/39 


0001 

SUE  ROUT  I  N  E  TERM'S  (X,  PN  TERM  ,TERR2t 

I  ) 

0002 

DIMENSION  TERF1C5CI ,TER«2<50> 

0003 

T  =■  •  5  *  X 

O.CO  4 

S«SIMT)*s2 

0005 

X0  =  1.  if N 

0006 

X  1- X  0  ♦  1* 

CC  C  7 

TERMl  (11=1, 

OCO  8 

TERW2(I>=1« 

0009 

JT=1 

0010 

DO  10  < -  2  f  5C 

001  I 

JT  =  - JT 

0012 

I  1  =  1-1 

0013 

TER  Ml ( I  ) =JT*S* M  1*P'CKP< XOt  II  >/FACK( 

1  !  >**2 

001 A 

TER M2 ( I >=jT*S**I !*FACKP(X1# 11 !/r ACKI 

T1  J**2 

0015 

IF  (AeS*TERM2(  I)  i-OStlCflC 

C  Cl  6 

10 

CONTINUE 

CC  1  7 

1=50 

CC1  8 

RETURN 

0019 

9 

1=1 

0020 

RETURN 

0021 

END 

FORTRAN 

IV  G  LEVEL 

18  TERMNS  OATE 

0001 

SUBROUTINE  TERMNS(TA,XtM  ,TEPNN.I ) 

0002 

P!7*';  NGI ..  .*■  -.’om 

Guu  3 

XN=X*N1 

000 A 

T  =0 • 5  *  T  A 

00  C  5 

S*SIN(.T  >**2 

0006 

JT=1 

C007 

DQ  10  1=1,50 

C  009 

JT=-JT 

0009 

T  ERMN(  I  >=JT+S**I *FACKP(XNf I )/FACK( I >**2 

0010 

I F ( A35(  T EPMN ( I) ) -1 • E-07 ) 9 , 10 , 10 

coil 

10 

CONTINUE 

CC1  2 

I  =50 

0013 

RETURN 

0014 

9 

I*! 

C  Cl  5 

RETURN 

CClfc 

END 

FORTRAN  IV  G  LEVEL  iS 


FCN 


DATE  =  71110 


0001 
CC  02 
C003 
COOA 

OCO  5 
OC  06 
OCC  7 


FUNCTION  ECN(X,XN> 

T  =  . 5  *  X 

Y  =  1  •  57C79652*>N 

FC N  =  *  63  661 97  5* ( PS IG ( XN  +  1 • M  ALOGC  S INI T  >>  +  •  5  772  1  5664  9  >  -T  AM  Y  > 

F  CN=-FCN 

RETURN 

END 


16/18/29 


18/18/29 
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F JRTRA N 
0001 


IV  G  LEVEL  16  BESJ 

SUP  f  C  LT  I  4  F.  I’ESJU  <  X  .\  ,  3  J  ,  C  P  J ,  f.  ,  !EF<I 


DATE  =  71110 


C 

c 

c 

c 

c 

c 


n  IS  SI  3.-CU1  V:c  C^PUhS  VALUES  fop  fractional  order 
t  ions  c  r  ti  e  f  its:  :  m. 


i e/13/39 


eFSSEL  S  FUNC7- 


COO? 

D  I  vr  rsic.ix  &J  (1031  ,  VJ( 

0003 

P !  14  1  f  L/C  5 

CCO<* 

c--r>.  E-r.4 

00?5 

N- 1 '  I M />. t  C> 

0006 

N 1  --N  <  1 

OOC  7 

W2-M  2 

CCC  6 

FN-XN-FLG at ( M 

CCC9 

20 

I  r  (  >: )  "0,30,24 

0010 

3  0 

IE«-2 

0011 

R  C  T  u  R  v 

0012 

29 

If( X-15. >22,3? ,?• 

0  013 

2? 

U7E  S T  =4  C  ♦  ♦  1 C t  *:<->-*  2/ 

0014 

GC  TC  *e 

0015 

34 

NU  ST  =9  0.  +  X/2. 

C  Cl  6 

3$ 

IF(  M -MESD4C.3F.  ,38 

0017 

38 

I  ER-4 

ccia 

RETURN 

CC1  s 

40 

IF(x«3t:*6«C  )  GC  71.  4 C 5 

0020 

CALL  B  JT  EPF(  X  ,  Fr.;  ,  D  ♦  3  1  ,i 

0021 

C  JT  -  0  . 

0022 
*1  ?  T 

DO  40  )  IK-1  ,M\ 

GC2  4 

400 

CCNT  I  s,l«E 

0023 

</?. 

00?  6 

GC  TO  3  5 

C02  7 

4  C  5 

CAl  »  X  ♦  FK,EJT  ) 

CC?  8 

35 

9  PR L V-0* 

OC2<5 

7  C 

MPA  X  - ME  ST 

0030 

PJtFON*  10 

0031 

100 

DJ  190  P**ZSBCiPr/»Xt* 

C032 

r  PI *1  . OF-79 

GC33 

F  M  -  •  ') 

0034 

130 

CC35 

Mli *-) 

0036 

2  10 

DO  lfc<  K---I.P2 

OC  3  7 

M  K“M- K 

0035 

MK’.-  w  j/  _  j 

003  S 

XPK-r  ICATI  vkh-FK 

0040 

B  j;  P*  i-2.  *>:PK*FP1/X-FM 

0041 

FP=r  n 

004  2 

160 

K  > 

B  JC 1  i  ^2  •"  (1  -  +FM  -F  N'l /X 

0044 

4  L  f  M  4  T  =  9  J  1 1  X  /  3  J  T 

0045 

*»0  J  .>  1  U  i  **  i  ,  2 

C04t 

bJ(  J;  1  -  f;  J  (  j?  )  /  A  L  r  h  A  T 

004  7 

3  01 

CCf<1  I  N i 

0C4  r 

IF  (  4  5  >  c  »? .» {  K  1  ;  -  ?  ?  c  2  V  1  -  A ! 

004  S 

10  2 

IM>-  Fy  =  “J(  M  1 

C050 

1  Erl  -  :1 

0010 

G»i  1  *J  20  j 

.*>  a  c  7 

201 

I  L'L  ---  " 

BESC  3'>70 
BESC04SC 

3ESC0500 

BCSC0520 

D  E  $  C  0  5  4  0 
PESC055- 


BESC  37?C 
BF  SC  36  30 


6ESC0850 


BE  SC0S70 


-88- 


FGRTRAN  IV  G  LEV  EL  19  St  $  J  PATE  *  71110 


00  5  3 

2*0 

090(1  )  r  -  E  J  (  2  1*F.\-*J(1)/X 

C  05  A 

IFIH1-2  )5CC,-»Clt^01 

0055 

AO  l 

DC  A  2C  L L r 2  »  f> 2 

0.C5  6 

ii:=Li“i 

0057 

CSDJ“FLCAT«LH><FN 

00  5  8 

•OBJ(  LI  3  -  3  J  <  L  L  1 3  ~  C  R !  j  1  *  E  J  (  L  L  5  /  X 

005  5 

A  20 

CCMIMiE 

0060 

£  CO 

R?  RRN 

0061 

END 

FORTRAN  XV  G  LEVEL 


TOG  UM  GATE  =  71110 


000  1 


c 

c 

c 

c 

c 

c 


SUOR  CLT  i  .\E  TC$UW  t  I  \  A  ,  CF  t  b  ,  7CSU ) 

THIS  SUBROUTINE  GENERATES  VALUES  FC®  C  r 
HALF  CCNE  AML‘  =  IFIX(£i  +  1>~0F 


CC02 

CI^ENSICN  PI  (100),  OP  HI  00) 

0003 

COUPLE  PRECISION  TER (2HQ * » PS ( 200 ) ,TAf TA2f SX2 -CSUiPIC 

0904 

P  10  0*2  141  5S2C3D+  01 

CC05 

7  A  - (  H A-DF)*FIC/> 60 

0  C06 

T  A2  =  T  A  /  2 

0007 

Ml*Ml 

c  c  c.  *> 

CC09 

y  i  j  r  i  **■  i 

M2-K*2M 

0010 

$X2*DSIMTA2M*2 

COlt 

TER ( 1  )=C  .1C  +  C1 

0012 

J  1  =  1 

0013 

00  20  02  ,K1 

0G1A 

JT~~  JT 

001  5 

TER(  1  )  =  JT*0ABS(7  FR(  1-1)  )MK-l*2l*(M«1-ll*SX2/l(  I-  1  >  *  (  1-1)  ) 

C  01  6 

20 

CCNT I MJE 

0C1  7 

DO  4  0  I  ORI  1  ,P2 

ooi  e 

TEM  I  I  )-C*0  +  01 

CCl  5 

40 

CCN7  IMJE 

0C2  0 

DSC  1  >  -0*  1Q*01 

0021 

DC  30  II-2t«2 

002  2 

os(  n )- c.of ci 

0023 

DO  11  12*1,11 

0024 

DSC  I1)*LS!  It  MTERU2  1  *TER  <  1 1- 12*  U/ 1 1 

0025 

11 

CCNT  )  MJE 

0026 

10 

CCNT  it U i. 

CC2  7 

o so*  o.-ho: 

002c 

DC  3C  KK*ltP? 

00?  P 

t>SU*9S'J*DSI  KK> 

0030 

30 

C  Ci^T  I  NR 

0331 

CALL  LFGrM  ITAiU-.Si  tPJ  fTP!) 

003? 

-  L5  I N  ( T  A  )  *  '2  1  ?  I  i  )  *Di;  I<  1  3*2  •  (M*  1 )  *SX  2*0  SU  ) 

0033 

P  E  ’♦  UR  N 

0C34 

t  NO 

19/18/39 


18/18/39 
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F0kTPA>4  Jv  C  S.EVF.L  13  SFJ\X  CAT  5  -  71110  IQ/18/3'5 


0001 

SUHJ  DC:'  IN:  S  V  JN X  i  L  ,  >. ,  SO  > 

c 

c 

/ 

c 

'HIS  Sv^FOUT  :.\£  GfNcRATTS  SPHERICAL  eZSSEL  3  FUNCTIONS  OF  THE  FIRST 

c 

K  J 

\T). 

r 

C 

C002 

Oi  v.-NSiCN  59  U  ) 

0003 

L»U  TO  NK»J.  r  L 

0004 

r  o ^ 

l  O  Y  :  0 

0006 

T'I  C  *3*  1'  1  i5.?t  5 

C  OCT 

Ter.  i * x *- f i c 

CCOtt 

tfs  '2  =  x-:*r:c 

occo 

30 

If  (  <.M  .  C  •  >Gl  TO  1 

001  c 

IF  tf,%  „EC .  *»>  o']  TO  2 

0011 

;u(  sr: j  "C • 

001  2 

00  0  20 

CO  l  3 

*1 

SB  (  N N  )  =  I  . 

OC 1  4 

GO  TO  20 

001  5 

1 

If  (  N  .GF.  :>  X  tT*  1  * )  GO  TO  1C 

0016 

IF  N  • F C -  11  GO  TO  3 

0017 

A-STNUI  /> 

0C1P 

IF-  f  N  •  NF  •  0  )  GC  T  ^  4 

CC19 

.SB(NN)«A 

002  0 

GO  • 0  20 

0021 

4 

IF  iX  ,CE.  FLOAT  .M  )  CC  fC  3 

C022 

IF  UBSI  Tf  ST1  J  .  LE  .  1  .  L-C6.  uK.  ABS(  TtST*  KLF.  1.  L--05)GC  TC  *5 

*7  "* 

r «  t  <*•  /  ' 

C  u*  4 

-  i  1  (HAr*A’UMAli/A/^ 

002  5 

IPX=1 

002  6 

46 

Nl  =  3  4M0.5  +  EXP{-FLC/n  <N-2>/<*>  >*X 

002  7 

F  *N*.  M 

0C2a 

U=? >K45 

002  0 

r!W  .S-3C 

0  03  C 

r-Ni  --0. 

0  C  3  1 

N  1  =  N  1  *  1 

0032 

00  5  r*!fM 

0033 

l. r  U*  -  2  • 

0  0?  4 

f IN  U-FN/X-FN1 

0035 

PN1=FN 

0C3  6 

\ 

FN-F IN 

CC3  7 

<>3tNN)=FN 

0039 

IF (  !PX«NE«1 ) GC  TC  55 

0039 

N=N-1 

0340 

53 

00  6  I- 1  ,N 

0  041 

L'*-U~2. 

0042 

FIN  =  *J*FN/X-  M 

0C4  3 

FN1  F\ 

004'* 

6 

F »\“F  IN 

004  5 

Sr<NN)^SB(NN)*7/FN 

004  6 

f-n  TO  2  0 

OC4-? 

2 

f*  *  C  S1MU  -x  *0.0  S  (  X)  )  /X/X 

c  o  a 

IF  If.  ,M.  n  GC  TO  7 

0045 

5  C  \  N  M  ’  2 

0C5C 

GO  TO  20 

00  6  3 

7 

U*l. 

00*? 

on  9  l  ~ 2  •  N 

-90- 


FO&TRAN  IV  C  LEvEL  13 


005  3 

U  =  LH  2. 

0054 

FN 1 =U*B / X-A 

0055 

A  =  3 

0C56 

p 

B  =F  N 1 

0057 

S B 1 J  N  i=a 

0058 

GC  -C  20 

0C59 

1C 

A*  1 «. 

CC60 

B '-I  * 

CCfcl 

X  ?=-*X*  X  /  2  • 

C06  2 

f -i-o« 

C  C  5  3 

p  ?=: *k+i 

0064 

It 

si-:  1 4i. 

0065 

02=  B2  +  2. 

CC66 

B  =  B :  1 X  2  /  B 1  /  3  2 

006  1 

A=A*B 

0063 

IF  1 AflS ( D/A  ) 

0065 

H  =1  . 

C  07  C 

PQ  12  1  *  1  *  N 

OC  7  1 

£> — B+2  • 

0072 

12 

A«A/B 

0073 

SP( N  N  >  *  A  *  X  *  *  N 

0074 

20 

CONTINUE 

0075 

return 

0076 

END 

SPJNX  DATE  =  71110 


.  1 «  E-6 )  GC  TO  11 


FORTRAN  1 V  G  LFV£l  10 


0001 
OC  0  2 
CCC2 
C  00  4 
CCC5 
ccce 

0  CO  7 

cccs 
oco<? 
0010 
001  1 
0012 
0013 

oou 

001  5 
001  6 
cc:  7 
001  ? 
CC1  *? 
C  02  0 
0C21 


SUO!<OCTINE  fiJTEr.YtX  , FN» 0 , 8 JT 1 ,CJT 2,  TN) 

DIMENSION  0J7115C ) »ej72( 50 
J  =  l  F  IX  (  1  *c*X*l  D.  > 

CALL  G^FFA(l,4FNfGXti jRI 
Xf2  * ( X / 2 • >**? 

JT  =1 

3  JT1 (1 ) *1./GX 
3JT.M  1 

or  ro  mn=2,vj 

JT=- JT 
MrJl  -s  MN-1 
MN2*  KN-2 

XMi\~rlUAT(MN) 

X*N5 =FLC AT {*N1 ) 

XMN?=-FLCA7  (MN2I 

BJTI  (  >  =  JT*ABS(  3  jTK  FM)  >*XE2/  (  X*M  ♦  (  XMN  1  +  FM  >  > 

&J72i  >N)«JT*ABS( ?JT2 (FM  )  )*XS2/ (  XMN  1  *  ( X.v\ 4 F\i  )  » 
If(ASS(RJTl(FM)  .Le.C.AN'C.A0S(3j:2(MNM.LE.C)  GC  TC  10 
20  CCNT  IMJE 
10  RETURN 
E  NO 


lb/lo/39 


-91 


FORTRAN 

IV  0  LEVEL 

15  LEGFI  DATE  =  71110  13/i6/3 

arc  i 

SLF.I  OLTINE  LEGF  S  <  \X  ,SF ,K » F ! ,CFI 1 

000  2 

mu  f.SICN  SHU  ,P  I  (  30  1  ,CP  U  50  1 

000  3 

PIT-*  3.1  4155265 

oco^ 

g=p; C/1  50 

0005 

II)/  C 

CCC6 

IFIf.X.EC.lflO.CP.N/.EC.GIGC  TO  305 

0007 

IK  it  X-  002,2,1 

core 

1 

ID/  1 

C  CCS 

f*  X  ”  1  6  w  -  N  X 

0C10 

x* ( i  x+m  i*g 

c:m 

00  *  c  3 

0012 

2 

xm  ;  x*  j-  >  *  o 

001  3 

IF(,  E$( X  )  .Lc-  1.F-C2  'CC  TC  3C5 

0C14 

-> 

cx^i  r.s:  x ) 

0  015 

X2“  >:<••  X 

0016 

p  1 1  M  - 1  . 

CC1  7 

CPI ( 1  )  -C  . 

0013 

r  I  (  )  =C  X 

001  a 

D?J<2  >*L  . 

C  02  0 

CO  30  L  =  3,6 

00?  1 

l  1  =  L-  1 

CC2? 

L2-1 -2 

00?  3 

F  Ll - F  LOAT ( LI ) 

002  4 

r  L2- F  L04T ( L2i 

0025 

PI  i  L>  =  *<2.  *-f  L2*l  .  ?»CX*PT  (Li  )-FL2*PI  <L2  H  /  <FL  2*  1.  ) 

002  6 

OP  If  L  >=L1  *<CX*P  I  (  L  >-rjHLl  i > /<CX**2-1  .  1 

CO?  7 

30 

COM  I  N'J i; 

002  b 

1*1  K*.L!:.6»GC  TG  103 

0  029 

CALL  SPJhX(4fX,S3> 

i .  -  t 

p»*  <  -  VM  '  t  /  t  4  -  S  «•'  i  '•  i  1 

0  0  31 

ph3- sg ( 3)  /( :4.*sen  » ) 

0  0  3? 

PH4 *  St<  ( 4  )  /  <  1  92  SE  ( 1 1  ) 

003  3 

!F< A  3  S( >1-0.02113,13,14 

0034 

13 

AO=SLV<T  U.*X2t<  1#/d#  +  ?#*X2/36C*  )  ) 

C  C  3  5 

GO  1C  13 

00  3  c 

14 

AO50rT(X/$IM X) > 

C03  7 

1  5 

A 1  =*  •  5*4C*  FH 2 

0033 

A?«*C*3.*C-.5*FH2+.  3Y5*PH2**2) 

00  J9 

A3=*0*15.*(-«5*°H**.75*PH2*PH2-*3I25*P*2«*3> 

CC4C 

00  ^0  >’  =  7  ,  N 

0041 

«!-►-! 

C  042 

M2=y-2 

C.rt* 

cF) ^  F LCA7 ( M 1 

C04  4 

FM?-F l C A  T ( F2) 

004  5 

<*v=f  I  r  a  t  t  iv-i  i 

0C4C 

Xrt*XNO.S 

0047 

X  X  ■=  X  N  *  X 

CC48 

P  i  (  p  1  *  }P?{  XX  1*  1  A0-A7/  XNv*« i-4.*A0/i  X^XN'^^4  H  ♦e.r  1 «  X  <  1  *  (  A]  /XM2.  * A2/  l 
1X*XN*  <3)  +  A2M  E./XX*  *?-?,«  )/XN**3> 

0049 

CPI  C  M=M  *«CX*PHM>-PI(  Mi  1  >/f  CX**2-i#  ) 

CC5C 

40 

COM  ?MJE 

:o6 1 

ivi 

;p<  IDX.NE.l ) G C  TC  I  03 

00*  <; 

JT  1 

CC5? 

K  f  -  -  1 

0  054 

N\-  1  f^C-tv  X 

005  5 

DO  1  ( i  f.  r  I  r  N 

CC 5  6 

J7 -- J f 

CCS  7 

«T=~K  T 
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FCRTR  AN  I  V  G  LEV  “L  1  P 


LECFI 


DATE 


71110 


005a 

?>MM  *iCT*:-i  (.<; 

CC5  9 

DHIIK  )-JT*()FI  (*> 

C  Ct-G 

102 

CONTINUE 

CC61 

r»n  to  1 02 

0062 

3C5 

no  2  C4  IL  =  1 t N 

0  C  fc  2 

p  i  ( m  - 1  • 

0064 

304 

CON1 1 NUE 

CC6  5 

DPIi 1 )=C. 

0054 

DPI*  ?  )  -  1  • 

C  06  7 

OPM  3  )*?  * 

i  n  v  7  n-/  „  k 

CC69 

OP  1  *  ICJ«(  ID-  1  )♦ « DPI  UC-D-PIC  IDI)/(  !C- 

CC7C 

3  03 

CCN*’  I  NUE 

0071 

IFCCX.NE.  130  )  G  C  TC  1C3 

OC72 

J$*l 

CC73 

KS  -•  1 

007* 

DO  : Ct  KP*1 «N 

CC  7  5 

J$  —  J$ 

f  C76 

K  $ -— K  $ 

0  r  7  7 

Pl( KP)=K$*Fl <KP) 

0076 

OP!  (  KP)*»JS«Dtl(KP) 

CC79 

304 

CCM  INUE 

003  0 

1  C  3 

RETURN 

0061 

END 

FORTRAN 

IV  C  LEVEL 

12  BY TERM  f  CATE  =  71110 

0001 

SUBROUTINE  BYTERFIX *FN ,C , BYT! f 3YT2  **M ) 

crc2 

DIMENSION  BY  T  1  ( 50  ) »  BYT2 ( 5C ) 

0C02 

HJ-1F IX (10. 4 S CRT (1. f3.*X*X)l 

CC04 

CALI  GPPPA ( 1 •  — FN  »GX  ,  IEP) 

000  5 

XF2* (X/2* )  *  *  2 

COCf 

JT-1 

0007 

BYT]  (  1)-1. /GX 

oooa 

9YT2 ( 1 ) "-FN/GX 

CC09 

DO  20  PN^2  tPJ 

OCIO 

JT=-JT 

0C1  1 

MP1*PP-1 

003  ? 

N’M2=FK-2 

0  0!  3 

XPP*FLC.*T(PMI 

0014 

XPM1*FLCATCF PI) 

0015 

?  =  T  L  C  A  T (PP2) 

C01  h 

GY  T  i 1 F M )-JT*fi0S(EYTl(PMl )  >*XE2/ ( XP M 1  * C XMM 1-F N )  ) 

OCi  7 

(jYT2(PtM)-JT*AP$(BYT2(PMl  )  )*XE2/(XPPl*A5S(XMP2“FM  ) 

OOIP 

I  F(  AH  S(  BYU  (MM  )  )  •  LE#  C.  ANC. ABS  (BYT  2(M,H)  )•  LE.  0)  GO  TO 

001  9 

20 

CONTINUE 

C02C 

10 

PbTURN 

0C2  1 

END 

16/1S/29 


16/16/39 
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fortran  : 

onoi 

OOC2 
C003 
CO  04 
GOO  5 
000 6 
OOC7 

coce. 

nco? 

COi  0 
001  1 
001? 
uCi  3 
00  14 
0015 
OCi  6 
001  7 
C  01  >3 
0  019 
002  0 
0021 
0022 
0  02  3 
GQ2  4 
00?  5 

002  6 
>  ,  >  / 

CC?  b 

0029 
003C 
0  021 
0032 
002  3 
0  03  4 
CC35 
OC36 
003  1 
CC36 
0039 

0  0“r  c 

00s  1 
004? 

0043 
00:4 
$C45 
0  C  4  6 
004  7 
C04S 
GO'  S 
C  05  0 
0051 


V  G  LEVEL  13  SFMMX  DATE  =  7111C  18/18/'9 

SU5f.  DtTfNf  t  L  *  X  »  S  H  ,  '•  £  H  t  D  £  2  S  )  > 

C  THIS  SrjRCLTINe  GENERATES  VALUES  rC?  HAfJKcL  S  FUNCTIONS. 

CC^P.CX  SH  CiC!  ,D3MtO,C:SCSi50) 

CIME\M(A  S3C50)  #5Y45C>  , V  1 3 < 5 1 , UT S < 5 ) 
call  SPJNXil  ,X,52> 

CALL  SPYNXU  ,  X,S  V) 

00  2°  f-Ml.l 

Sh(N\)=CVrLXtSE(NM  ,-SYiNM  ) 

20  CONTINUE 

OSh  {  1  )--S)  (  2  5 
00  i  ^  1 1  -  c-  *  L 

i;SH<  LlS  -  SH!  LL-1  S -FLOAT  <U)-*SMIU/X 

30  CONTINUE 
ir-(L.LE.r7JGC  rC  31 
DC  130  KP^CCtL 
QKD*  <P-0  *5 

H  S  =  X 

SX$=H</CRD 

ALPHS  =ALCG( ( 1 •♦SCPT <  1  .-$XS**2  H/SXS  ) 

TXS'SQRTd.-SXS^SXSl 

cxs*i»/tx$ 

CX S2  *CX  S*(.XS 

EXPN-GRDM ALPLS-TANH(ALPHS) ) 

UTS(i  >  -  C  X  S  *  (  3  .~5.*CX$2>/24 

UTS  i  2)=CXS2*<  ei.+CXS2*(-4<:2*3e5*CXS2)  1/1152 

UTS  (  3 ) - C  X  S*C  X  S2*  (  10  37  5  +  CX52* (-369gC3  *CXS2*<  76 5765-42542 5*CXS2  >>>/'• 
i  1 4  7  2  0 

VTS<l)*CXS*C-9+7*CXS2)/24 

y  i  >  isi  y. 

VF$<  3)sCX$*CXS2  *(-4  2542  5*  CX$2M451737*CX  $2*1-2  83  575*475475  *CYS?)  »  ) 
I/M472C 
OL-  JS  -  1. 

CEY5  =  U 
0  L  J  P  $  =  1  • 

0 1 Y  P  :S  =  1  • 

JT  =  1 

DO  19  I P= 1 »  2 
JT  --  JT 
CU-ORC**  I  P 

OE JS  *0c  J$*LTS ( I P)/09 
UE YS*DE  YS* JT*LTS(  I  P) /CP 

oejPs=Gr jps*vts ( ip) /ce 
otr ps =Dl y fs * jt*vt  s( i P)/ce 
19  C  C  ft  T I N  U  2 

H  l  i*-0.5*EXF(-2.^6XPNi)  *(CEJPS/CEYPS*D  =  JS/0EYS  5 
DSCS(KPMCKPLX(1  .,R  IN) 

DSOS(KP)-C. 5/FS*SlNM AL  FHS  >  *CEYPS*CSCS(  KP  l/CtrYS 
OSOS  (  KP  i  =■— 0 S C S  (  KF  ) 

100  C  C N T I  My  E 

31  CO  4C  1  iC«)  tl 

D.snsc  I K. )  *CSH  (  I  KI/5H(  IK) 

40  CONTINUE 
X L  TUP N 
END 
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which  project  the  ends  of  a  thin  biconical  antenna  is  chosen  to  simulate  the  actual  sphere- 
centered  thin  dipole  The  conical  antenna  is  driven  at  its  junction  with  the  sphere  by  a 
rotationally  symmetric  electric  field  maintained  across  the  gap  by  a  biconical  transmission 
line  excited  by  the  TEM  mode.  The  attractive  features  of  this  model  include  the  fact  that 
it  has  surfaces  that  permit  a  simple  specification  of  boundary  conditions  and,  hence,  a 
rigorous  formulation  for  the  electromagnetic  fields  and  a  shape  such  that  its  properties 
should  come  reasonably  close  to  those  of  a  modified  cylindrical  antenna  as  the  cone  angle 
becomes  quite  small. 

The  measurements  of  both  input  admittances  and  current  distributions  on  modified 
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curves  for  modified  cylindrical  and  conical  antennas  involve  only  slight  shifts  suggests  that 
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An  infinite  set  of  algebraic  equations  was  solved  numerically  in  Volume  II  for  small 
cone  angles  Comparisons  were  made  between  the  modified  conical  antenna  and  its  limiting 
biconical  antenna  which  provides  both  an  extrapolatory  numerical  check  for  the  modified 
conical  antenna  with  shrinking  central  sphere  and  an  understanding  of  the  underlying 
physical  phenomena.  Theoretical  and  experimental  results  are  in  very  good  agreement 
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